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1.  Summary  of  the  significant  work  accomplished 

Bright  Incoherent  Solitons 

Self-consistent  modal  theory 

The  self-consistent  modal,  theory  was  developed  by  our  group  as  an  alternative  way  to  describe  in¬ 
coherent  solitons  in  non-instantaneous  nonlinear  media  (Phys.  Rev.  Lett.,  vol.  79,  pp.  4990,  1997). 
Unlike  the  coherent  density  approach  which  is  by  nature  more  appropriate  for  dynamical  evolution 
(under  arbitrary  initial  conditions),  the  modal  theory  is  better  suited  in  identifying  incoherent  soli¬ 
tons,  their  range  of  existence,  and  coherence  properties.  This  method  relies  on  the  idea  that  every 
incoherent  spatial  soliton  must  be  a  combination  of  the  modes  of  its  own  self-induced  waveguide.  The 
total  intensity  of  this  soliton  is  the  sum  of  all  the  modal  intensities.  The  range  of  existence  of  inco¬ 
herent  solitons  in  saturable  media  such  as  photorefractives  was  also  obtained.  For  bright  incoherent 
solitons  we  found  that  the  coherence  increases  at  the  low  intensity  tails.  Similar  results  were  obtained 
numerically  using  the  coherent  density  method  (Opt.  Lett.,  vol.  23.  pp.  418,  1998). 


Incoherent  solitons  in  Log  nonlinear  saturable  systems 

Theoretically,  incoherent  spatial  solitons  were  first  identified  in  Log  saturable  nonlinear  media  using  the 
coherent  density  approach  (Christodoulides  et  al.  Opt.  Lett.  vol.  22,  pp.  1080,  1997).  The  incoherent 
bright  solitons  allowed  in  this  system  exhibit  a  Gaussian  intensity  profile  and  they  exist  below  a  certain 
cut-off.  It  was  thus  natural  to  ask  whether  the  self-consistent  modal  theory  could  also  apply  in  this 
system.  We  have  found  that  the  modes  of  these  incoherent  Gaussian  solitons  can  be  obtained  in 
terms  of  Gauss-Hermite  polynomials  as  long  as  the  mode-occupancy  function  is  Poissonian.  Gaussian 
soliton  solutions  were  also  obtained  in  two  dimensions.  Unlike  their  coherent  counterparts,  they  could 
also  be  elliptical  depending  on  their  2-D  coherence  characteristics.  This  implies  to  some  extent  that 
incoherent  spatial  solitons  may  be  more  versatile  in  terms  of  shape  as  long  as  their  coherence  properties 
are  tailored.  Above  all,  this  result  demonstrated  that  the  two  methods  give  the  same  results  (in  the 
Log  case)  and  thus  they  are  equivalent.  Their  equivalence  was  later  formally  established  on  the  basis 
of  Mercer’s  theorem  and  the  associated  Karhunen-Loeve  expansions.  In  fact  the  modal  theory  is  also 
dynamic  (it  can  be  used  for  dynamic  evolution)  unlike  what  we  thought  in  the  very  beginning. 

Coherence  properies  of  incoherent  spatial  solitons  in  non-instantaneous  Kerr  media 

In  this  part  of  our  research,  we  have  investigated  the  coherence  properties  of  multimode  incoherent 
spatial  solitons  in  non-instantaneous  Kerr  nonlinear  media  (Physical  Review  E,  vol.  59,  pp.  1193, 
1999).  In  doing  so  we  employed  the  self-consistent  multimode  description  first  developed  by  our  group 
(Phys.  Rev.  Lett.,  vol.  79,  pp.  4990,  1997).  In  general  we  have  found  that  the  intensity  profile  of 
these  incoherent  Kerr  states  is  of  the  sech^(3:/a:o)  type.  The  modal  constituents  or  modes  of  these 
solitons  were  found  in  terms  of  associated  Legendre  polynomial  functions.  The  spatial  width  of  these 
incoherent  states  was  explicitly  obtained  as  function  of  the  total  number  of  modes  supported  and  the 
amount  of  nonlinearity  provided  by  the  system.  Even  more  importantly,  the  complex  coherence  factor 
associated  with  these  Kerr  incoherent  solitons  was  also  obtained  in  closed  form.  Our  results  were  in 
excellent  agreement  with  previous  semi-analytical  results  (Phys.  Rev.  Lett.,  vol.  79,  pp.  4990,  1997). 
One  interesting  aspect  that  emerged  from  this  study  is  the  fact  that  the  coherence  manifolds  of  such 
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multimoded  states  are  ’’topologically”  the  same  depending  on  whether  the  number  of  modes  involved 
is  odd  or  even. 

Elliptic  incoherent  solitons  in  saturable  nonlinear  media 

In  Optics  Letters,  vol.  25,  pp.  972  (2000),  we  identified  elliptic  incoherent  solitons  in  isotropic 
saturable  nonlinear  media.  We  were  able  to  do  this  by  employing  a  two-dimensional  version  of  the 
coherent  density  approach  and  by  assuming  that  the  material  nonlinearity  depends  on  the  optical 
intensity,  I,  in  a  way  similar  to  that  in  photorefr actives.  We  also  found  that,  even  in  this  case,  elliptic 
solitary  states  exist,  provided  that  their  mutual  coherence  function  is  anisotropic.  Therefore  these 
soliton  states  can  possibly  be  excited  from  anisotropic  incoherent  sources  such  as  edge-emitting  LED  s 
or  lasers  operated  below  threshold.  The  propagation  dynamics  of  this  new  class  of  solitons  were  further 
studied  by  use  of  numerical  simulations.  Collisions  between  two  such  elliptic  incoherent  solitons  were 
also  investigated.  In  particular,  we  showed  that  in  certain  collision  regimes  the  intensity  ellipse  for 
the  two  solitons  rotates,  whereas  at  the  same  time  their  centers  of  gravity  tend  to  revolve  around 
each  other.  This  interaction  behavior  is  possible  even  though  these  states  are  launched  parallel  to  the 
propagation  axis. 


Coherence  enhancement  of  spatially  incoherent  light  beams  through  soliton  interactions 

We  have  shown  that  it  is  possible  to  enhance  the  spatial  coherence  length  of  a  bright  partially  incoherent 
signal  beam  through  interactions  with  a  coherent  (or  an  incoherent)  dark  spatial  soliton.  We  have 
demonstrated  that  during  this  interaction  part  of  the  incoherent  bright  beam  is  trapped  within  the 
dark  or  gray  notch  of  the  controlling  soliton,  thus  forming  a  sharp  intensity  spike.  In  this  region, 
the  correlation  length  increases  by  at  least  2  orders  of  magnitude.  Therefore,  incoherent  light  can 
be  effectively  cooled  (its  entropy  being  reduced)  at  any  arbitrarily  chosen  point  upon  a  partially 
incoherent  wave  front  by  use  of  a  dark  spatial  soliton.  This  is  the  only  passive  system  that  we  know  of 
that  exhibits  an  increase  in  both  the  local  intensity  and  the  local  correlation  distance  simultaneously. 
These  results  were  presented  in  Optics  Letters,  vol.  25,  pp.  826  (2000). 

Incoherent  dark  solitons 

Theory  of  incoherent  dark  solitons 

The  theory  of  incoherent  dark  solitons  was  developed  soon  after  their  observation  (Phys.  Rev.  Lett., 
vol.  80,  pp.  5113,  1998).  To  do  so  we  used  the  self-consistent  modal  theory  (Phys.  Rev.  Lett.,  vol. 
79,  pp.  4990,  1997).  We  have  found  that  these  dark  incoherent  solitons  are  comprised  from  a  belt  (a 
continuum)  of  radiation  modes  as  well  as  bound  states.  Unlike  coherent  dark  solitons  which  involve 
only  an  odd  radiation  mode  at  cut-off,  the  radiation  belt  contains  both  odd  and  even  modes.  Similarly, 
the  bound  states  can  be  odd  or  even.  The  presence  of  even  radiation  modes  and  even  bound  states 
explains  why  the  incoherent  dark  solitons  are  in  fact  gray.  On  the  other  hand,  close  inspection  of  their 
modal  structure  reveals  that  in  the  dark  notch,  the  odd  modes  dominate  which  in  turn  explains  the 
’phase-memory  effect” .  The  coherence  properties  found  theoretically  were  in  agreement  with  those 
obtained  earlier  numerically  (Opt.  Lett.,  vol.  23.  pp.  418,  1998).  Another  interesting  outcome  of  this 
study  was  the  fact  that  these  solutions  are  by  no  means  unique.  Instead  they  depend  on  the  choice 
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of  the  radiation  mode  distribution  function.  This  problem  is  equivalent  to  Manakov  system  of  infinite 
dimensions. 

Incoherent  dark  solitons-phase  memory  effects 

The  dynamics  of  incoherent  solitons  in  non-instantaneous  nonlinear  media  were  investigated  using  the 
so-called  coherent  density  approach,  first  developed  by  our  group  (Phys.  Rev.  Lett.,  vol.  78,  pp.  646, 
1997).  In  this  approach,  we  assume  that  the  input  can  be  modeled  as  a  continuum  of  sources,  which  are 
mutually  incoherent  with  respect  to  each  other.  We  further  assumme,  that  at  the  origin,  the  source- 
fluctuations  are  stationary  which  is  typically  the  case  in  statistical  optics.  This  method  was  augmented 
with  a  version  of  the  Van  Cittert-Zernike  theorem  in  order  to  monitor  the  evolution  of  the  coherence 
during  propagation.  This  particular  version  is  better  suited  for  multi-Fourier  split  method  techniques, 
used  in  our  numerical  computations.  In  Opt.  Lett.,  vol.  23.  pp.  418,  1998  we  have  simulated  the 
dynamics  of  incoherent  solitons  in  biased  photorefractive  crystals.  Surprisingly  enough  we  have  found 
that  dark  incoherent  quasi-solitons  can  form  after  a  certain  distance  of  propagation.  In  addition,  these 
incoherent  dark  solitons  were  found  to  be  in  fact  gray.  What  was  even  more  interesting  was  the  fact 
that  the  nonlinear  system  could  exhibit  phase  memory  effects  in  spite  of  the  incoherency.  At  higher 
biased  voltages  (higher  nonlinearities),  an  even  field  depression  was  found  to  lead  to  doublets  whereas 
an  odd  one  to  triplets,  pretty  much  the  same  way  one  would  have  expected  from  a  coherent  system.  It 
was  thus  natural  to  ask  the  following:  (i)  why  are  the  fundamental  incoherent  dark  solitons  gray?  (ii) 
how  can  the  imposed  phase  be  remembered  in  the  midst  of  random  phase  fluctuations?  Subsequently 
this  behavior  was  observed  experimentally  in  SEN  crystals  (Science,  vol.  280,  pp.  889,  1998).  In  this 
experiment  1-D  stripe  dark  (gray)  solitons  as  well  as  2-D  gray  vortex  solitons  (of  topological  charge 
1)  were  observed. 

Incoherent  dark  solitons  Y-splitting  and  ’’phase  memory”  effects 

The  properties  of  these  incoherent  dark  Y-soliton  doublets  were  further  investigated  in  a  subsequent 
paper  of  ours  (Physical  Review  E,  vol.  59,  pp.  R4777,  1999)  where  we  have  also  provided  the  first 
experimental  observation  of  this  very  intriguing  behavior.  In  this  paper  we  have  shown  that  the 
dynamics  of  these  self-trapped  entities  are  associated  with  strong  ’’phase-memory”  effects  that  are 
otherwise  absent  in  the  linear  regime  (i.e.,  during  diffraction).  The  numerical  simulations  carried  out 
using  the  coherent  density  theory  (developed  by  our  group,  Phys.  Rev.  Lett.,  vol.  78,  pp.  646,  1997) 
revealed  that  the  Y-splitting  process  is  rather  insensitive  to  the  degree  of  spatial  coherence.  More 
specifically  we  found  that  over  a  wide  range  of  parameters,  the  Y-splitting  is  approximately  the  same 
irrespective  of  the  degree  of  spatial  coherence.  These  conclusions  were  successfully  validated  in  the 
experimental  part  of  our  study. 

Waveguides  formed  by  incoherent  dark  solitons 

We  have  demonstrated  experimentally  optical  guidance  of  coherent  light  beams  using  incoherent  light 
(Optics  Letters,  vol.  24,  pp.  1160,  1999).  This  was  made  possible  using  either  dark  incoherent  solitons 
(previously  predicted  and  verified  by  our  group,  Opt.  Lett.,  vol.  23.  pp.  418,  1998;  Science,  vol. 
280,  pp.  889,  1998)  as  well  as  Y-dark  soliton  doublets  as  described  in  the  first  section  of  this  report. 
In  all  cases,  the  waveguides  were  probed  at  higher  wavelengths  and  it  was  found  that  they  were  in 
fact  single-moded.  This  is  in  agreement  with  our  earlier  predictions  concerning  the  modal  structure  of 
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incoherent  dark  solitons  (Phys.  Rev.  Lett.,  vol.  80,  pp.  5113,  1998).  These  experiments  suggest  that 
it  is  possible  to  control  high-power  laser  beams  with  low-power  incoherent  light  sources  such  as  LEDs. 

Vector  solitons 

Multi-component  vector  solitons 

It  has  been  previously  shown  that  saturable  Manakov  systems  can  be  implemented  in  photorefractive 
crystals  (Christodoulides  et  al,  Appl.  Phys.  Lett.,  vol.  68,  1763,  1996).  In  principle  these  Manakov- 
like  systems  can  involve  an  arbitrary  number  of  components  as  long  as  they  are  mutually  incoherent. 
It  was  also  known  for  some  time  (at  least  in  the  Kerr  regime)  that  bimodal  vector  solitons  should  also 
be  possible  (multi-component,  multi-humped),  Christodoulides  and  Joseph,  Opt.  Lett.  13,  pp.  53, 
1988,  TVatnik  and  Sipe,  Phys.  Rev.  A  38,  pp.  2001,  1988.  These  soliton  states  were  observed  for  the 
first  time  in  saturable  nonlinear  crystals  (Phys.  Rev.  Lett.,  vol.  80,  pp.  4657,  1998).  Their  theory  was 
also  formulated.  Among  other  things,  it  was  found  that  double  and  triple  hump  solitons  are  possible. 
Solitons  which  can  contain  no  ground  state  were  also  identified  and  observed. 


Multi-component  t^vo- dimensional  solitons  carrying  topological  charge 

In  Optics  Letters,  vol.  25,  pp.  61  (2000)  we  have  proposed  for  the  first  time  two-dimensional  vector- 
solitons  for  which  each  component  carries  a  different  topological  charge  (or  “spin”).  We  have  found 
that  this  family  is  characterized  by  a  unique  triple-point  phase  diagram,  which  is  completely  absent 
in  the  two-component  case.  Moreover,  we  have  shown  that  these  higher-dimensional  vector  structures 
are  composed  of  lower-dimensional  building  blocks.  With  soliton  collisions  in  mind,  it  is  clear  that 
the  spin,  the  multimode  nature,  and  the  multihump  structure  offer  new  opportunities  for  interactions 
between  two-dimensional  vector  solitons.  In  this  work  we  employed  the  thresholding  nonlinearity  and 
made  use  of  the  self-consistency  principle.  The  core  ideas  and  findings  presented  are  expected  to 
be  universal.  These  composite  solitons  can  provide  insight  into  ways  to  realize  these  interactions  in 
materials  with  saturable  nonlinearities.  In  addition,  we  have  suggested  (Physical  Review  Letters,  vol. 
84,  pp.  1164,  2000)  composite  solitons  that  carry  topological  charge:  multicomponent  two-dimensional 
[2+lD]  vector  (Manakov-like)  solitons  for  which  at  least  one  component  carries  topological  charge. 
These  multimode  solitons  can  have  a  single  hump  or  exhibit  a  multihump  structure.  The  spin 
carried  by  these  multimode  composite  solitons  suggests  2-D  soliton  interactions  in  which  the  particlelike 
behavior  includes  spin,  in  addition  to  effective  mass,  linear,  and  angular  momenta. 

Observation  of  two-dimensional  multimode  solitons 

In  Optics  Letters,  vol  25,  pp.  1113  (2000),  we  presented  the  first  experimental  observation  of  (2+1)- 
dimensional  multimode  (composite)  solitons.  A  single-hump  component  and  a  double-hump  (dipole- 
type)  component  are  jointly  self-trapped  as  a  composite  soliton  in  a  biased  photorefractive  crystal. 
We  observe  that  the  two  modes  self-trap  when  they  are  launched  together,  forming  the  composite 
soliton.  However,  the  stand-alone  components  do  not  trap  on  their  own  for  the  given  value  of  the 
nonlinearity  that  supports  the  composite  soliton.  This  self-trapping  of  multimode  (2+1)  D  solitons 
opens  up  the  possibility  of  distortionless  image  transmission  through  highly  nonlinear  self-focusing 
media.  If  the  envelope  of  a  highly  multimode  soliton  is  modulated  to  contain  an  image  (superimposed 
on  the  intensity  profile  of  the  soliton),  then  this  image  can  be  transmitted  through  the  self- focusing 
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medium  and  remain  unchanged  through  propagation.  This  method  is  in  contradistinction  to  image 
transmission  through  a  multimode  fiber,  for  which  the  fiber  modes  are  coherent  with  each  other,  yet 
they  propagate  at  different  velocities  (intermodal  dispersion)  and  thus  destroy  the  image. 

Modulation  instability 

Modulation  instability  of  incoherent  beams  in  noninstantaneous  nonlinear  media 

In  Physical  Review  Letters,  vol.  84,  pp.  467  (2000),  we  demonstrated,  analytically  and  numerically, 
the  existence  of  incoherent  MI,  that  is,  modulation  instability  of  incoherent  wave  packets.  Incoherent 
modulation  instability  exists  when  the  nonlinear  index  change  exceeds  a  well-defined  threshold  that  is 
determined  by  the  amount  of  incoherence.  This  is  in  a  marked  difference  with  coherent  MI,  since  there 
does  not  exist  a  similar  threshold  for  coherent  MI.  The  nonlinear  index  change  is  determined  by  the 
spatial  degree  of  coherence  (angular  power  spectrum  of  the  source).  We  used  analytical  and  numerical 
methods  to  study  the  properties  of  incoherent  MI  in  a  general  self-focusing  noninstantaneous  medium. 
Assuming  the  input  beams  have  Lorentzian  angular  power  spectra,  we  were  able  to  find  closed  form 
solution  for  incoherent  MI  when  the  nonlinearity  is  either  of  the  Kerr  or  saturable  type.  We  confirmed 
our  results  with  numerical  simulations. 

Bright  spatial  solitons  on  a  partially  incoherent  background 

We  have  demonstrated,  both  theoretically  and  experimentally,  that  stable  partially  coherent  antidark 
solitons  can  exist  in  noninstantaneous  nonlinear  media.  We  have  analyzed  the  modal  structure  of 
these  antidark  states  in  the  case  where  the  nonlinearity  is  of  the  Kerr  type.  Our  analysis  indicates 
that  these  solitons  always  involve  a  set  of  discrete  bound  states  as  well  as  a  continuum  of  odd  and  even 
radiation  modes.  Experimental  results  and  computer  simulations  indicate  that  the  instability  affecting 
an  antidark  soliton  can  be  totally  eliminated  by  properly  increasing  the  incoherence  of  its  background 
beam  above  the  threshold  of  incoherent  MI.  These  antidark  solitons  are  the  only  known  solitons  that 
have  nonzero  intensity  everywhere  in  space  yet  they  propagate  in  a  stable  fashion  in  a  self-focusing 
medium.  We  published  these  findings  in  Physical  Review  Letters,  vol.  84,  pp.  2374  (2000). 


Stimulated  Raman  Scattering 

Stimulated  Raman  interactions  in  massive  wavelength-division-multiplexed  systems 

In  addition  to  our  work  in  photorefractive  crystals,  we  have  explored  stimulated  Raman  interactions 
(SRS)  in  massive  wavelength-division-multiplexed  (WDM)  systems.  This  technology  shows  consider¬ 
able  promise  in  terms  of  exploiting  the  large  spectral  window  offered  by  today’s  low  loss  optical  fibers. 
In  general,  the  SRS  power  exchange  among  WDM  channels  is  described  by  a  rather  involved  system 
of  differential  equations  of  the  Lotka-Volterra  type.  This  problem  is  of  considerable  complexity  since 
its  interconnectivity  increases  with  the  square  of  the  number  of  channels  involved.  For  example,  in  a 
100-channel  system,  the  number  of  SRS  cross-talk  flow  streams  is  approximately  5000!  Yet,  this  rather 
involved  problem  can  exhibit  a  closed  form  solution  as  demonstrated  by  our  group  (Optics  Letters, 
vol.  24,  pp.  735,  1999).  In  this  paper  we  found  that  solitary-like  Raman  states  are  possible  in  such 
massive  wavelength-division-multiplexed  systems.  These  self-similar  states  can  propagate  undistorted 
among  channels  at  constant  velocity.  The  evolution  of  these  states  under  the  action  of  noise  and  during 
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collisions  was  also  investigated.  In  another  study  (to  be  published  in  IEEE  Photonics  Technology  Let¬ 
ters,  October  issue,  1999),  we  have  proved  that  spectral  inversion  techniques  can  totally  eliminate  SRS 
cross-talk  effects  which  are  known  to  severely  limit  the  performance  of  high  capacity  WDM  systems. 
The  result  was  obtained  by  exploiting  the  reflection  symmetries  of  the  underlying  evolution  equations 
and  it  is  therefore  general.  This  work  was  done  in  collaboration  with  Lucent  Technologies  and  Tyco 
Submarine  Systems. 
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2.  Computational  methodology /algorithms 


The  equation  describing  light  propagation  in  nonlinear  media  is  obtained  by  the  solution  of  Maxwell’s 
equations.  By  the  solution  of  the  Maxwell  set  the  evolution  of  incoherent  light  inside  a  nonlinear 
medium  is  governed  by 


+  ^x,j 


dx 


^  dy  )  2  \  dx^ 


+  F(/)[/,=0, 


(1) 


where  Uj  are  the  so-called  coherent  density  components,  and  0x,j,  dyj  are  the  direction  angles  of 
each  component  with  respect  to  the  propagation  2:  axis.  The  total  intensity  of  the  beam  is  obtained 
from  I  =  where  N  is  the  number  of  coherent  components.  F{I)  determines  the  change 

in  the  refractive  index  caused  by  the  material  nonlinearity.  For  a  Kerr  medium  F{I)  =  /,  whereas 
for  saturable  nonlinearities  (as  is  the  case  for  photorefractive  crystals)  F{I)  =  -1/(1  +  7).  The 
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initial  contribution  of  each  component  to  the  total  intensity  is  given  by  the  normalized  angular  power 
spectrum  GN{0^,dy)  of  the  incoherent  source.  Most  often  Gjv(0x,^y)  is  assumed  to  have  a  Gaussian 
or  Lorentzian  distribution. 

In  order  to  solve  Eqs.  (1),  we  first  discretize  the  {x,y)  domain.  We  then  use  the  split-step 
Fourier  method  to  solve  the  corresponding  discretized  equations.  To  demonstrate  the  split-step  Fourier 
technique,  we  write  the  equation  describing  each  coherent  component  as: 

^  =  {L  +  N)U,  (2) 

oz 

where  L  is  a  linear  differential  operator  given  by 

^  ^  2  \dx^  5y2  J 

and  N  is  the  nonlinear  operator 

N  -  iF{I). 

The  exact  solution  of  eq.  (2)  can  be  found,  by  integrating  with  respect  to  z 

U{z  +  h,x,y)  =  exp[h{L  + N)]U{z,x,y).  (3) 

where  h  is  the  step  size.  According  to  the  Cambel-Baker-Hausdorff  formula,  the  operator  on  the  right 
hand  side  of  this  equation  satisfies  the  following  relation 

/  h?  \ 

exp{hL)  exp{hN)  =  exp  lh{L  +  N)  +  —[L,  N]  +  —[L  -  N,[L,N]]-\ - J  , 

where  the  commutation  operator  is  defined  as  [a,  b]  =  ab  -  ba.  Although  in  principle  the  operators  L 
and  N  do  not  commute,  when  the  step  h  is  small  enough  the  higher  order  terms,  i.e.,  Oih?),  are  small 
perturbations  in  this  expansion.  Then,  eq.  (3)  can  be  approximated  by 

U{z  +  h,x,y)  =  exp{hL)exp{hN)U{z,x,y).  (4) 

Now  these  two  operators  can  be  applied  separately  into  U{z,x,y)  in  order  to  propagate  a  single  step 
in  z.  The  action  of  the  nonlinear  operator,  N,  in  Eq.  (4)  is  straightforward,  i.e., 

Ui{z,x,y)  =  exp{hN)U{z,x,y). 

On  the  other  hand,  the  linear,  differential  part  of  Eq.  (4)  can  be  solved  using  a  two-dimensional  discrete 
Fourier  transform 

U{z  +  h,x,y)  =  exp{hL)U{z)  -  F~^ {exp[hL{iqa;,iqy)]F[Ui{z,x,y)]},  (5) 

where  L{iqx,iqy)  is  the  Fourier  transform  of  operator  L,  and  F,  F~^  denote  the  operators  of  the 
forward  and  the  inverse  Fourier  transform,  respectively.  This  method  is  implemented  on  the  CRAY 
T3E  at  the  Pittsburgh  Supercomputer  center.  A  sample  code  copy  can  be  found  in  Appendix  A. 
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A 


Sample  code 


!  2D  Incoherent  bright  solitons 

!  Fourier  definitions  and  main  arrays 

parameter (Np=101) 
parameter (NRA=1024) 
parameter (Pi=3 . 141593) 
complex  SIG(NRA,Np),  B(NRA) 
real  ak(NRA) 

real  Pow(NRA) ,B1(NRA) ,AU(NRA) 
real  WFFTC(4*NRA+15) ,  CPY(2*NRA) 
c 

c  /Program  Definitions/ 

integer  StepNum,  SaveNum 
integer  Nx,  Npp 
c 

real  IntX,IntY,beta,h 
real  ThetaO,  Strech,  rt 
real  alpha2 
real  xw,  xf 
c 

c  /Loop  counters  / 

integer  st,  jj,  ECL 
c 

c  /Save  definitions  / 

integer*4  sh5,  sh6,  first_r 
character*!  sheets,  sheets, list3 
integer  grid_ste,  dat_ste 
c 

external  fftcf,  fftcb 
c 

c  /Common  blocks/ 

common  /discr/  nx 

common/components/Npp 

common/discr2/IntX 

common  /nonl/beta 

common /ArbUnits/  ThetaO,  Strech 

common/TotInt/rt 

common/steps/h , dz 

common/LogPar/ alpha2 

common/widths/ xw , xf 

common/outr/sh5 , shS , f irst_r , sheets , sheets , lists 
c 

c  /Constants/ 
c 
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-)  IntX  and  IntY  are  the  WHOLE  interval  lengths.  :-) 


Npp=Np 

Nx=NRA 

IntX=200.e0 

call  InitParam 
print*,  ’alpha2  =’,  alpha2 

print*, ’t=  ’,  h 

stepNuin=95000 

grid_ste=1900 

dat_ste=200 

SaveNuni=stepNuin/grid_ste 

list3=’e’ 

pow( : )=0. 

B1=0. 

call  fftci(NRA,WFFTC) 

call  store_init  !  Initializing  the  storage  parameters 

call  getBigSignal(SIG)  !  ...  and  getting  the  signal  -  new  version 

call  getPower(SIG,Pow) 

B=Pow 

call  fftcf (NRA,B,B) 

B=B/sqrt (real (Nx) ) 

Bl=abs(B) 

AU=0. 

AU ( 1 : nx/2) =B1 (Nx/2+1 : nx) 

AU(nx/2+l :nx)=Bl(l :nx/2) 

call  f illreaKPow,  AU,IntX,0,dz) 
call  getak(ak) 

write(*,*)  ’Initializing  FFT. .. ’ 

Print*,  ’Good  Luck  ! ! ! ’ 

The  external  loop  using  the  (E)xternal  (C)ounter  (L)oop  SaveNum 
do  105  ECL  =  1 , SaveNum 

:-)  The  first  half  step 

do  45  jj=l,Np 

call  fftcf (NRA,SIG(: ,jj) ,SIG(: ,jj) .WFFTC,CPY) 
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call  dispersion(SIG(: , jj) ,h/2. ,ak) 

! 

call  fftcb(NRA,SIG(: , jj) ,SIG( : , j j) ,WFFTC,CPY) 

45  continue 

call  Power(SIG,Pow) 

call  nonlinearity(SIG,h,Pow,Np) 

! 

do  65  st=2,grid_ste 

!  /Action  /  The  main  portion  of  the  Difr.-  Nonl.  loop,  grid_ste 
! 

do  85  jj=  1,  Np  !  The  component  loop 

! 

call  fftcf (NRA,SIG(: ,jj).SIG(: ,jj),WFFTC,CPY) 
call  dispersion(SIG( : , j j) ,h,ak) 

! 

call  fftcb(NRA,SIG(:,jj),SIG(:,jj),WFFTC,CPY) 

85  continue 

! 

call  Power (SIGjPow) 

call  nonlinearity(SIG,h,Pow,Np) 

! 

65  continue 

j  The  last  half  step 

! 

do  35  jj=l,Np  !  the  component  loop 

j 

call  fftcf (NRA,SIG( : ,jj) ,SIG(: ,jj) ,WFFTC,CPY) 

call  dispersion(SIG(: , jj) ,h/2. ,ak) 

call  f ftcb(NRA,SIG( : , j j) ,SIG( : , j j) ,WFFTC,CPY) 

35  continue 

call  Power (SIG,  Pow) 

B=Pow 

call  fftcf (NRA,B,B) 

B=B/sqrt (real (Nx) ) 

Bl=abs(B) 

AU=0. 

AU(l:nx/2)=Bl(Nx/2+l:nx) 

AU (nx/2+1 : nx) =B1 (1 : nx/2) 
call  f illreal (Pow , AU , IntX , ECL*gr id_ste , dz) 


105 

continue 

stop 

end 

c 

c  :-) 

c 

;-)  :-) 

:-)  :-)  :-)  :-) 
PROCEDURES 

:-)  :-) 

c  :-) 

:-)  :-) 

:-)  :-) 
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c 


subroutine  dispersion(b,d,ak) 
c 

c 

c  -  Variable  declarations  - 

c 

common  /discr/  nx 
integer  nx 
complex  b(nx) 
real  alc(nx) 
complex  c 

real  d  !  d  is  the  step  size 

! 

c=cmplx(0. ,1.) 
c 

c  -  calculation  of  the  dispersion 

b= (b/sqrt (real (nx) ) ) *cexp (c+d^ak) 
return 
end 

subroutine  getak(ak) 
c 

c 

c  -  Variable  declarations  - 

c 

parameter (Pi=3 . 141593) 
common  /discr/  nx 
common/di s cr2/ IntX 
integer  nx 
real  IntX 
real  ak(nx) 


real  dkx 
real  kx 
integer  nx2,i 
integer  kx2c 
c 

nx2=nx/2+l 
dkx=2 . eO*Pi/IntX 
c 

c  -  Begin  nested  loops  to  calculate  the  frequency  distribution 

c 

do  67  i=l,nx 

if  (i  <=  nx2)  kx2  =  (i-1) 
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if  (i  >  nx2)  kx2  =  (-nx  -  1  +  i) 
kx=real (kx2) *dkx 
ak(i)=((-5.e-l)*(kx**2  )) 

67  continue 

return 
end 
c 

subrout ine  nonlinearity (b , d , pow , Npp) 

common  /discr/  nx 

integer  nx.Npp 

common  /nonl/beta 

real  beta 

common  /Totint/rt 

real  rt 

complex  b(nx,Npp) 
real  pow(nx) 
real  d 

real  nonlin(nx) 
integer  i 
complex  c 

c 

c=cmplx(0. ,1.) 

! 

nonlin=(beta)*(l .+rt)/ (pow+1 . ) 
nonl in=nonl in*d 
do  60  i=l,  Npp 

b( : , i)=b( : , i)*cexp(c*nonlin( : )) 

60  continue 

return 
end 

subroutine  fillreal  (a,Bl,IntX,rec,step) 

common/discr/nx 

integer  nx.rec 

real  a(nx),  Bl(nx),IntX 

common/outr/ sh5 , sh6 , f irst_r , sheets , sheet6 , listS 

common/TotInt/rt 

integer*4  first_r,  sh5,  sh6 

character*!  sheets, sheet6, lists 

character* 10  alfabet 

character*8  filename 

integer  ii 

real  rt,x,step 

real  field,  fmax 

c - 
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c 

open(uiiit=29,file=’m’ //lists//' .dat’ ,status= 'UNKNOWN’ , 
$position= ’ APPEND ’ ) 
f max=0 . eO 

do  225  mx=nx/2-20,  nx/2+20 
f ield=a(mx) 

if  (f ield.ge.fmax)  fmax=field 
225  continue 

print* ,  fmax 

write (29,16)  real (rec) *step/l . e-4 , fmax 
close  (;mit=29) 

16  format(f 10.5, f 10.5) 

c  fmax=rt 

alfabet=' 1234567890’ 

Q - - - begin - C 

c 

if (f irst_r .eq.O)  then 
first_r=l 
goto  101 
end  if 

if  (sh6.lt. 9)  then 
sh6=sh6+l 

sheet6=alf abet (sh6 : sh6) 
else 

sheet6='0' 

sh6=0 

if  (sh5.lt. 9) then 
sh5=sh5+l 

sheet5=alf abet (sh5 : sh5) 
else 

print*,  'I  can  not  save  that  many  files' 


c 

C=======END  IF  SHEET3=========C 

end  if 
c 

C======eND  if  SHEET4==========C 

end  if 


c  /  saving  the  INTENSITY  distrubution  as  a  grid  file  / 
c 

101  filename  =  ’i’//list3//sheet5//sheet6//’ .dat’ 
open(unit=31,  f ile=f ilename,  status=’ unknown’) 
print*, 'writing  ’,  filename 
do  122  ii=nx/2,  nx 
x=-IntX/2 . +real (ii-1) *IntX/real (nx) 
write(31,*)  x,  a(ii) 

122  continue 

close (unit=31) 
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filename  =  ’mV/list3//sheet5//sheet6//' .dat’ 
open(unit=32,  file=f ilename,  status=’ unknown') 
print* , 'writing  ',  filename 
do  ii=nx/2+l,  nx 

x=-IntX/2 . +real (ii-1) *IntX/real (nx) 
write(32,*)  x,  Bl(ii) 
end  do 

close (unit=32) 
c 

!!  open(unit=71,  f ile=’w'//list3//' .dat' ,  status= ' UNKNOWN ' , 

!  !  $position=' APPEND') 

!!  write(71,10)  a(nx/2-39:nx/2+40) 

10  format (80f 10.5) 

! !  close (unit=71) 

! 

!!  open(unit=71 ,  file='s'//list3//' .dat' ,  status= ' UNKNOWN ' , 

!!  $position=' APPEND' ) 

!!  write(71,10)  Bl(nx/2+l:nx/2+80) 

! !  close (unit=71) 

c 

return 

end 


subroutine  store_init 

common/outr/sh5 ,  sh6,  first_r,  sheets, sheets, list3 

integer*4  sh5,  shS,  first_r 
character*!  sheets,  sheets, list3 

sh5=0 

shS=0 

f irst_r=0 

sheet5='0' 

sheetS='0' 

return 

end 


c  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  =-)  =-)  =-)  =-) 

subroutine  getBigSignal (data) 

common/components/Npp 

common/discr/nx 

common/discr2/IntX 

integer  Nx,Npp 


15 


NC  is  the  niomber  of  components  in  x  and  y 
< —  counters  for  ThetaX  and  ThetaY 


common/ArbUnits/  ThetaO,  Strech 
common/TotInt/rt 
common/widths/xw , xf 
real  ThetaO,  Strech,  rt 
complex  data(nx,Npp) 
real  xw,  xf 
complex  c 
real  xste,  IntX 
real  x,  ampl,fi 
integer  i,jcomp  ! 

integer  itx  ! 

real  p 
real  thetaX 
real  InThX,  dThX 
real  rt,  rHat,  RGN 
real  sum  ,suml 
real  pow,fiO,  aO,al,pi 

!  aO  is  the  radius  of  the  supergauss ian  background  beam 

real  fx,  ftx  !  dummy  arguments  for  the  functions  below 

f (f x)=exp(-0 . 5* ( (f x/ (70 . *xw) ) **6) ) 

f (f x) = (2 . / (exp (f x/6 . ) +exp (-f x/6 . ) ) ) **32 
f  (f  x)  =2 .  /  (exp  (f  x/xw)  +exp  (-f  x/xw) ) 

GN (f tx) =exp (- (f tx**2  ) /Theta0**2) 

Print*,  ’Getting  the  signal 


!  Intensity  normalization 


!  ! 
I  I 


!c 
!  c 
! 


!  ! 


! 


c=cmplx(0. ,1.) 
pi=3. 141593 
suml=0. 
fi0=0. 

aO=(3.25e-l)*IntX 
al=(4.5e-l)*IntX 
xste=IntX/real (nx-1) 
print*,  ’xw  =’,  xw 
print*,  ’rt  =  ’ ,  rt 

InThX=Strech*ThetaO 


If(Npp.ne.l)  then  !  incoherent  case 

! 

dThX= (2 . *InThX) /real (Npp-1) 

! 

sum=0. 

do  41  i=l,  Npp 

p=-InThX+real(i-l)*dThX 
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siiin=sum+GN  (p) 

!!  print*,  ’GN  =  GN(p),  p 

41  continue 

!!  print*,  ’  sum  =  sum 

rHat  =  (rt/sum) 
print*,  ’rhat  =  rhat 

! 

else  !  coherent  case 

InThX=0 . ;  dThX=0 . ; rhat=rt ; 
end  if 

jcomp=l  !  Initialization  of  the  component  counter 

do  32  itx=l,Npp 

! 

ThetaX  =  -InThX+dThX*real(itx-l) 

!!  Print*,  ’ThetaX,  ThetaY  ThetaX,  ThetaY 

RGN=sqrt (rHat*GN (ThetaX) ) 

! 

suml=suml+RGN**2  !  suml  is  being  used  for  energy  check  only 

! 

do  20  i=l,Nx 

x=-IntX/2.eO  +(real(i-l))*xste 

! 

ampl=RGN*f (x) 

! ! !  Getting  the  phase  of  the  dark  beam  (s)  (-: 


f i=ThetaX*x 

! 

data(i, jcomp)=ampl*cexp(c*f i) 

! 

19  continue 

20  continue 

! 

!!  call  fillgrid(data(: , jcomp)) 

jcomp=jcomp+l 

32  continue 

print* , ’suml=’ ,  suml 

return 

end 

j  - 

subroutine  spectrum(a,b) 
real  a(nx) 
common  /discr/nx 
integer  nx,i 
complex  b(nx) 

! 
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do  i  =1,  nx 
b(nx)=a(nx) 
end  do 

! 

return 

end 

c  :-)  :-)  :-)  :-)  :-)  :-)  :-)  =-)  :-)  =-)  •-'>  ■-'>  •") 

subroutine  InitParam 
parameter (Pi=3. 141593) 
common/ArbUnits/  ThetaO,  Strech 
common/TotInt/rt 
common/ steps/h , dz 
common/LogPar/  alpha2 
common/nonl/beta 
common/widths/ xw , xf 

real  ThetaO,  Strech,  rt,dz  !  ThetaO  is  the  width  of  Int.  Distr. 

real  alpha2 
real  xw,xf 

real  xO,  h,  lng,Ld,LDi, ratio, iratio 
!  Ing  is  the  crystal  lengthh  in  mikrometers 

! 

!  Strech  determines  the  half  interval 
!  width  for  thetaX  and  thetaY.  InThX  =  Stech*ThetaO. 

!  ThetaO  is  in  degrees.  All  distance  are  in  mikrometers. 

!  rt  is  the  ratio  between  the  beam  intensity  and  the  background  intensity 
! 

real  lambda,  nO 

lambda=0 .5145 

n0=2 . 35 

x0=10. 

xw=l . 

xf=0. 

! 

Theta0=0.55  !  degrees 

Strech=2 .7 
rt=10. 
beta=-37 . 

!  Transformation  of  ThetaO  in  radians  and  in  a.  u. 

!  In  fact  the  normalized  value  of  ThetaO  is  the  V  parameter 

!  that  occurs  in  the  logarithmic  paper 

! 

ThetaO  =(4 . *Pi*Pi*xO*ThetaO*nO) / (lambda*360 . ) 
print*,  ’ThetaO  =  ’,  thetaO 

alpha2=ThetaO**2+l .  !  The  incoherence  parameter  in  normalised  units 

lng=l.e+4  !  1  cm  in  mikrons 

dz=lng/95000. 
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Ld=(nO*2.*Pi/lambda)*(xO**2)  !  the  coherent  difraction  length 

!  since  the  equation  is  normalized  to  it 
Ldi=Ld/sqrt(alpha2)  !  The  incoherent  diffraction  length 

ratio=dz/Ld 

iratio=dz/Ldi  !  The  fraction  of  the  ’incoherent  difraction  length 

! 

if  (iratio.gt. 0.025)  then 

print*,  ’your  step  is  too  big’,  h,  iratio 
end  if 
h=ratio 

print*,  ’h,  Ld,  ratio,Ldi, iratio  =  ’ 
print*,  h,  Ld,  ratio,  Ldi,  iratio 

! 

print* ,  ’ alpha2  =  ’ ,  alpha2 
return 
end 

! 

!  :-)  ;-)  :-)  ;-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  =-) 


subroutine  getPower (data,  PW) 

common/ components/Npp 

common/discr/nx 

integer  nx,Npp 

complex  data(nx,Npp) 

real  PW(nx) 

integer  comp,i 

! 

PW(:)=0. 
do  comp=l,Npp 

PW(:)  =  PW(:)+cabs(data(: ,comp))**2 
end  do 
return 
end 

!  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  =-)  =-) 

!  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  =-)  =-) 

- - 

subroutine  Power (data,  PW) 

common/components/Npp 

common/discr/nx 

integer  nx,Npp 

complex  data(nx,Npp) 

real  PW(nx) 

integer  comp,i 

! 

PW(:)=0. 

data=data/ sqrt (real (nx) ) 
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do  comp=l,Npp 

PW(:)  =  PW(;)+cabs(data(: ,comp))**2 
end  do 
return 
end 

!  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  :-)  =-)  =-) 
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Abstract 


This  thesis  investigates  the  dynamical  and  soliton-like  behavior  of  spatially-partially 
incoherent  light  beams  in  non-instantaneous  nonlinear  materials.  Self-trapping  of 
incoherent  beams  and  their  coherence  properties  are  studied  in  great  detail.  Two 
complimentary  approaches  are  developed  in  order  to  explain  the  observed  incoherent 
self-focusing  behavior  in  biased  photorefractive  media;  namely  the  coherent  density 
approach  and  the  self-consistent  incoherent  multimode  method. 

It  is  shown  that  under  appropriate  initial  conditions,  bright  as  well  as  dark¬ 
like  incoherent  solitons  are  possible  in  saturable  nonlinear  media.  Our  numerical 
simulations  have  demonstrated  that  the  coherence  properties  of  these  beams  are 
signi  cantly  affected  by  the  self-trapping  process.  More  speci  cally,  in  the  case 
of  bright  beams,  we  have  found  that  the  coherence  length  remains  approximately 
constant  around  the  center  of  the  beam,  whereas  it  increases  at  the  margins.  An 
exact  (one/two  dimensional)  Gaussian  solution  was  also  obtained  in  the  case  where 
the  nonlinearity  is  of  the  logarithmic  type.  These  incoherent  Gaussian  solitons  can 
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exist  as  long  as  their  spatial  width  is  appropriately  related  to  the  strength  of  the 
nonlinearity  and  the  width  of  the  angular  power  spectrum  of  the  incoherent  source. 
Moreover,  two  dimensional  incoherent  solitons  can  be  elliptical  or  circular  depending 
on  whether  the  angular  power  spectrum  of  the  incoherent  source  is  symmetric  or 
not. 

Our  analysis  has  demonstrated  that  dark  incoherent  solitons  are  possible  only 
if  a  TT-phase  shift  is  initially  imposed  on  the  wave  front,  and  that  they  are  in  fact 
gray.  In  this  case,  the  coherence  length  has  been  found  to  be  higher  within  the  dark 
notch,  with  a  depression  at  the  center.  These  dark  soliton  entities  involve  radiation 
modes  as  well  as  bound  states.  Depending  on  the  initial  conditions,  an  even  or  an 
odd  number  of  incoherent  dark-like  structures  can  be  obtained  in  a  self-defocusing 
environment.  It  was  found  that  over  a  wide  range  of  parameters,  the  Y-splitting  is 
approximately  the  same  irrespective  of  spatial  coherence.  It  is  important  to  note 
that  the  behavior  of  bright/dark  incoherent  solitons  is  fundamentally  different  from 
that  of  their  coherent  counterparts.  For  example,  unlike  a  coherent  gray  soliton, 
an  incoherent  dark  soliton  does  not  exhibit  a  transverse  velocity  in  spite  of  its 
grayness.  Furthermore,  it  is  shown  that  the  evolution  of  incoherent  dark  solitons  in 
non-instantaneous  nonhnear  media  is  associated  with  strong  phase-memory  effects 
which  are  otherwise  absent  in  the  linear  regime. 
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Chapter  1 


Introduction 


Solitons  represent  ubiquitous  entities  in  the  general  eld  of  nonlinear  science.  Thus 
far,  they  have  been  successfully  identi  ed  in  several  branches  of  physics,  such  as  for 
example  hydrodynamics,  nonlinear  optics,  plasma  and  solid  state  physics,  magma 
ow  and  high  energy  physics  [1],[2].  By  de  nition,  solitons  are  stationary  solutions 
to  nonUnear  dispersive  partial  differential  equations  which  also  retain  their  identity 
after  a  collision  event.  Losely  speaking,  the  term  can  also  be  used  in  conjuction  with 
self-trapped  wave-packets  which  are  described  by  non-integrable  systems  where  the 
outcome  of  a  collision  is  typically  more  involved. 

Within  the  discipline  of  nonlinear  optics,  solitons  are  typically  subdivided  into 
two  basic  categories:  spatial  and  spatio-temporal  solitons  [3].  In  the  case  of  spatio- 
temporal  solitons,  perhaps  the  best  known  example  is  the  nonlinear  Schroedinger 
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soliton  which  occurs  in  optical  bers  [4] .  Fiber  solitons,  which  may  play  a  crucial  role 
in  tomorrow  s  high-speed  ber  optic  networks,  have  been  intensively  investigated 
in  the  last  two  decades  at  both  the  experimental  and  theoretical  level  [5].  Other 
types  of  spatio-temporal  solitons  can  also  be  mentioned.  These  include  for  example, 
self- induced-transparency  solitons,  three-wave  mixing  and  Raman  solitons,  ampli  er 
solitons  [3]  as  well  as  Bragg  or  gap  solitons  [6]. 

An  optical  beam  naturally  tends  to  broaden  during  propagation  in  a  homogenous 
medium.  This  broadening  can  occur  in  time,  in  space,  or  in  both.  In  the  temporal 
domain,  the  broadening  is  due  to  chromatic  dispersion  where  each  frequency  com¬ 
ponent  of  the  pulse  (wave-packet)  bears  different  phase  velocities.  On  the  other 
hand,  in  the  spatial  domain  the  broadening  is  due  to  diffraction  which  is  similar  to 
dispersion.  A  spatial  beam  can  be  viewed  as  a  linear  superposition  of  plane  waves 
(akin  to  frequency  components  of  a  pulse)  all  having  the  same  wavenumber,  yet  each 
propagating  at  a  different  angle,  let  say  a,  with  respect  to  the  propagation  axis.  As 
for  temporal  pulses,  each  plane-wave  component  accumulates  a  different  phase,  since 
its  propagation  constant  is  proportional  to  cos(a).  This  results  in  phase  building 
differences  between  each  spatial  frequency  and  thus  the  beam  diffracts.  In  general, 
narrower  beams  have  broader  plane-wave  specta  and  diverge  faster.  Waveguiding  is 
a  commonly  used  method  to  eliminate  diffraction.  When  the  interference  between 
total  internal  re  ections  from  the  boundaries  of  a  waveguide  is  constructive,  the 
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beam  gets  trapped  inside  the  waveguide  thus  forming  a  guided  mode  . 

An  astonishing  possibility  to  counteract  diffraction  effect  manifests  itself  when 
the  material  is  nonlinear,  as  rst  suggested  by  Chiao,  Garmire  and  Townes  [7]. 
While  diffraction  causes  the  spatial  beam  to  spread  and  decrease  in  amplitude,  the 
nonlinear  effects  tend  to  make  the  beam  steepen  and  become  narrower.  In  the 
right  situation,  these  opposing  effects  will  complement  each  other,  leaving  the  beam 
propagating  without  any  distortion.  As  a  matter  of  fact,  the  optical  beam  creates  its 
own  waveguide  as  it  travels  through  the  material.  Thus  a  balance  between  diffraction 
and  nonlinearity  is  reached  and  as  a  result  the  beam  remains  localized  and  preserves 
its  shape  forming  a  spatial  soliton  . 

It  is  interesting  to  note,  that  the  rst  solitons  to  be  investigated  in  the  eld 
of  nonlinear  optics  were  in  fact  spatial  solitons.  Shortly  after  the  discovery  of  the 
laser,  Chiao,  Garmire  and  Townes  suggested  the  possibility  of  self- trapped  optical 
beams  or  solitons  in  nonlinear  Kerr  media  [7].  In  this  case,  the  self-trapping  of  such 
a  needle  of  light  in  the  bulk  of  a  nonlinear  medium  is  made  possible  by  coun¬ 
teracting  the  hnear  process  of  diffraction  through  the  lensing  effect  of  self-focusing 
[8], [9].  Like  their  spatio-temporal  counterparts,  spatial  solitons  also  show  consider¬ 
able  promise  in  terms  of  practical  applications.  Several  applications,  ranging  from 
all-optical  switches  to  dynamic  optical  interconnects  and  routers  may  be  envisioned. 
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operating  at  strategic  locations  in  a  high-speed  optical  network  [10], [11].  Until  re¬ 
cently,  however,  all  efforts  aimed  to  harness  this  potential  were  met  with  failure. 
The  reason  behind  this  difficulty  is  primarily  twofold.  First  of  all,  the  observation 
of  Kerr  spatial  solitons  always  requires  exceedingly  high  power  levels  [12].  This  is 
because  the  nonlinear  index  change  in  this  case  is  proportional  to  the  hght  inten¬ 
sity.  For  example,  a  soliton  of  approximately  10  microns  in  diameter,  will  require 
more  than  200  kWatts,  if  it  is  to  be  observed  say  in  a  highly  transparent  glass 
medium.  This  situation  often  becomes  further  complicated  when  other  nonhnear 
processes  are  invited  in  (because  of  high  power)  such  as  stimulated  Raman  and 
Brillouin  scattering  as  well  as  two-photon  absorption  which  tend  to  spoil  the  soh- 
ton  itself  [12].  Moreover,  because  of  these  high  power  requirements,  lasers  can  only 
be  used  under  pulsed  conditions.  The  second  reason  behind  this  difficulty  is  even 
more  fundamental.  More  sped  cally,  Kerr  solitons  in  two  transverse  dimensions 
are  highly  unstable  since  they  will  either  catastrophically  collapse  or  break  up  in 
multiple  laments  [13], [14].  Thus  far,  bright  Kerr  solitons  were  only  observed  in 
one-dimensional  con  gurations  where  slab  waveguides  have  been  used  to  compen¬ 
sate  for  the  other  dimension  [10], [11].  This  explains  why  two-dimensional  spatial 
solitons  have  remained  elusive  for  the  last  thirty  years  or  so. 

This  picture  totally  changed  when  it  was  rst  predicted  that  spatial  solitons  may 
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in  fact  be  possible  in  biased  photorefractive  crystals  [15], [16].  At  the  time,  photore- 
fractives  have  been  traditionally  regarded  as  the  material  of  choice  for  real  time 
holography,  phase  conjugation  and  optical  data  storage  [17],  [18].  Nevertheless,  any 
notion  of  solitons  in  this  material  system  was  entirely  unheard  of.  Shortly  after, 
the  rst  successful  observation  of  what  are  by  now  better  known  as  quasi-steady- 
state  photorefractive  spatial  solitons  has  been  reported  in  strontium  barium  niobate 
(SBN)  crystals  [19], [20].  As  their  name  implies,  these  solitons  occur  in  transient, 
i.e.  before  the  space  charge  eld  settles  down.  Even  more  importantly,  these  soliton 
states  were  found  to  be  stable  or  robust  in  both  transverse  dimensions  and  they  were 
observed  at  extremely  low  power  levels  ,  i.e.  microwatts.  Stability  of  two  dimen¬ 
sional  self-trapping  is  due  to  the  fact  that  the  nonlinearity  happens  to  be  saturating 
which  means  that  the  induced  index  change  can  not  exceed  a  maximum  value.  In 
turn,  the  run-away  interaction  that  gives  rise  to  catastrophic  collapse  in  Kerr 
media  can  be  fully  arrested  since  the  induced  waveguide  eventually  becomes  wider. 
As  expected,  the  event  revived  considerable  interest  within  the  optics  community 
in  the  hope  that  spatial  solitons  may  now  be  within  reach  [21], [22].  Subsequently, 
the  so-called  steady  state  screening  solitons  have  been  independently  predicted  by 
Christodoulides  et  al  and  Segev  et  al  [23], [24]  and  were  experimentally  observed 
in  both  one  and  two  transverse  dimensions  in  SBN  [25], [26].  These  latter  entities 
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appear  after  the  space  charge  has  reached  its  steady-state.  At  the  same  time,  photo¬ 
voltaic  spatial  solitons  were  predicted  [27]  and  observed  in  lithium  niobate  [28], [29]. 
In  all  cases,  photorefractive  solitons  were  found  to  be  stable  against  spatial  pertur¬ 
bations.  Apart  from  fundamental  differences  associated  with  the  mechanism  behind 
their  existence,  what  greatly  distinguishes  the  class  of  photorefractive  solitons  from 
their  Kerr  counterparts,  is  the  required  power  level.  Unlike  Kerr  solitons  which 
typically  require  hundreds  of  kWatts  [10], [11], [12],  photorefractive  solitons  can  be 
observed  using  microwatt  or  even  lower  optical  power  levels  [19], [20].  Evidently, 
there  is  a  difference  of  approximately  11  orders  of  magnitude.  Even  for  the  newly 
discovered  family  of  chi-2  solitons  [30],  which  exhibit  stable  self-trapping  in  both 
transverse  dimensions,  one  still  has  to  use  hundreds  of  kWatts.  Another  important 
distinction,  arises  from  the  fact  that  photorefractive  solitons  can  write  permanent 
waveguides  in  the  bulk  of  a  photorefraetive  crystal,  which  in  turn  can  be  used  to 
guide  other  more  intense  beams  at  less  photosensitive  wavelengths.  In  this  manner, 
a  microwatt  soliton  in  the  blue  can  guide  and  control  a  1  Watt  beam  at  1-1.55 
micrometers,  intended  for  telecommunication  purposes  [28]- [34]. 

Very  recently,  for  the  rst  time,  incoherent  light  was  found  to  self- trap  itself  in 
a  biased  photorefractive  crystal  [35], [36].  This  very  interesting  twist  of  events,  now 
suggests  that  nonlinear  optics  no  longer  needs  to  belong  to  the  exclusive  domain  of 
powerful  lasers.  This  is  the  rst  time  we  know  of,  that  incoherent  light  determined 
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the  outcome  of  a  nonlinear  optics  experiment.  This  development  is  expected  to 
have  important  scienti  c  and  technological  implications.  Furthermore,  this  is  the 
rst  time  in  Physics  that  an  incoherent  wave-packet,  that  is  a  wave-packet  whose 
phase  varies  randomly  in  space  and  time,  has  been  self-trapped. 

This  dissertation  is  dedicated  to  theoretical  investigation  of  a  number  of  funda¬ 
mental  issues  associated  with  this  newly  discovered  class  of  incoherent  solitons. 
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1.1  Spatial  solitons  in  saturable  media 

Photorefractive  solitons  are  generally  classi  ed  into  three  generic  categories  each 
resulting  from  a  different  nonlinear  mechanism  which  is  essentially  saturable.  Sat¬ 
uration  of  nonlinearity  indicates  that  the  optically  induced  refractive  index  change 
can  not  exceed  a  maximum  value.  Saturating  nonlinearities  lead  to  more  compli¬ 
cated  systems  as  opposed  to  Kerr  nonlinearities.  However,  as  mentioned  before, 
the  presence  of  saturation  solves  the  problem  of  instabihty  associated  with  the  Kerr 
nonlinearity  in  (2-1- 1)D.  In  what  follow  we  provide  an  overview  of  the  photorefractive 
soliton  types  by  giving  special  emphasis  to  the  screening  solitons. 

1.1.1  Screening  solitons 

Screening  solitons  occur  in  biased  photorefractive  media  after  the  space  charge  eld 
has  reached  steady  state  [23], [24].  Basically  the  principle  behind  screening  solitons  is 
rather  simple.  At  relatively  high  bias  eld  strengths,  the  current  in  the  photorefrac¬ 
tive  sample  is  dominated  by  drift  whereas  the  rather  weak  diffusion  component  can 
be  neglected.  In  this  case,  the  photo  excited  carriers  migrate  in  such  a  way  that  the 
conductivity  of  the  photorefractive  crystal  becomes  higher  in  the  illuminated  regions 
whereas  is  lower  in  the  dark.  As  a  result,  the  absolute  value  of  the  space  charge 
eld  is  lower  (screened  out)  in  the  illuminated  regions  [23], [24].  In  turn,  through 
the  electrooptic  or  Pockels  effect,  the  refractive  index  is  modi  ed  in  proportion  to 
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the  local  space-charge  eld.  When  the  proportionality  coefficient  is  fairly  large  and 
negative,  this  leads  to  a  large  index  change  in  the  dark  regions,  whereas  the  refrac¬ 
tive  index  in  the  illuminated  regions  remains  almost  unaffected.  Consequently,  this 
leads  to  a  graded-index  waveguide  which  can  then  support  a  nonlinear  mode  or 
what  is  better  known  as  a  bright  spatial  screening  soliton.  Similarly,  dark  spatial 
solitons  are  possible  in  biased  photorefractives  provided  the  polarity  of  the  bias  eld 
is  reversed.  In  this  latter  case,  the  refractive  index  is  lower  in  the  illuminated  regions 
whereas  is  higher  within  the  dark  notch. 

Screening  solitons  were  independently  predicted  by  Segev  et  al.  and  Christodoulides 
et  al.  [23], [24].  In  these  studies,  bright,  dark  and  gray  one-dimensional  spatial 
screening  solitons  were  found  to  be  possible  in  biased  photorefractive  crystals.  More¬ 
over  it  was  shown  that  for  a  given  physical  system,  these  states  depend  only  on  two 
parameters,  namely  the  strength  of  the  applied  bias  eld  Eq  as  well  as  the  so-called 
intensity  ratio  r  or  p,  i.e.  the  ratio  of  the  optical  peak  intensity  with  respect  to 
the  dark  irradiance  of  the  crystal.  The  symbol  r  represents  the  maximum  in¬ 
tensity  ratio  of  a  bright  soliton  whereas  p  stands  for  the  maximum  dark  intensity 
ratio.  Shortly  after,  the  existence  of  bright  screening  solitons  was  experimentally 
con  rmed  in  SBN  crystals  by  Segev  s  group  [25], [26].  These  solitons  were  found  to 
be  stable  in  both  one  and  two  transverse  dimensions.  Concurrently,  dark  screening 
solitons  and  Y-sphtting  gray  solitons  were  observed  in  a  reverse  biased  SBN  crystal 
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[32],  [33], [34].  The  evolution  equations  of  bright  sohtons  were  also  investigated  [37] 
under  arbitrary  input  conditions.  In  this  study  [37],  various  issues  were  considered 
including  the  stability  of  these  states,  the  effects  arising  from  loss,  bias  and  the 
intensity  ratio  as  well  as  the  interaction  forces  among  in-phase  and  out  of  phase 
coherent  solitons.  In  another  theoretical  study,  the  effects  of  diffusion  on  these  soli- 
tons  were  considered  and  it  was  predicted  that  self-bending  should  occur  [38], [39]. 
This  self-bending  effect  was  later  successfully  veri  ed  [26].  The  guiding  properties 
of  the  circular  waveguides  permanently  induced  by  such  spatial  solitons  were  also 
studied  experimentally  [40],  along  with  the  rst  experiments  on  collisions  between 
pairs  of  photorefractive  solitons. 

At  the  same  time,  screening  solitons  were  also  investigated  by  Iturbe-Castillo 
et  al.  in  Bismuth  Titanium  Oxide  (BTO)  crystals  [41].  In  their  studies,  induced 
modulational  instability  [42]  was  observed  as  well  as  bright  and  dark  solitons  were 
studied  in  biased  BTO  crystals  [43].  A  theoretical  model  capable  of  describing  the 
process  of  modulational  instability  in  biased  photorefractives  was  later  developed 
by  Carvalho  et  al.  [44].  The  effects  arising  from  the  presence  of  optical  activity 
on  solitons  in  BTO  crystals  were  also  theoretically  investigated  and  it  was  found 
that  in  fact  soliton  behavior  in  this  system  can  only  be  claimed  at  relatively  short 
distances  and  at  small  intensity  ratios  [45], [46]. 
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In  addition  to  single  screening  solitons,  one-dimensional  vector-screening  soli- 
tons  have  been  predicted  in  crystals  of  appropriate  symmetry  [47].  Unlike  their 
Kerr  counterparts,  these  solitons  can  now  be  both  self-coupled  or  cross-coupled  de¬ 
pending  on  the  properties  of  the  electrooptic  tensor  [47].  Other  studies  have  also 
addressed  vector  interactions  among  solitons  [48]  as  well  as  the  possibility  of  ob¬ 
serving  bright-dark  vector  soliton  pairs  [49].  In  a  recent  study,  incoherently  coupled 
soliton  pairs  have  also  been  proposed  [50].  These  pairs  are  possible  provided  that  the 
two  carrier  beams  share  the  same  polarization  and  are  mutually  incoherent.  They 
can  appear  in  bright-bright,  dark-dark  as  well  as  in  dark-bright  con  gurations.  The 
existence  of  such  incoherent  soliton  pairs  has  been  recently  con  rmed  experimentally 
[51], [52], [53].  In  these  experiments,  bright-bright  and  dark-bright  pairs  have  been 
observed  in  accord  with  theory.  In  the  same  spirit,  incoherent  collisions  among  soh- 
tons  have  been  investigated  in  both  one  and  two  transverse  dimensions  [54], [55], [56]. 
These  latter  experiments  have  opened  up  a  new  avenue  in  terms  of  controlling  the 
interaction  forces  between  solitons. 

Another  important  development  in  this  area  includes  the  rst  experimental 
demonstration  of  self-trapping  in  InP:Fe  [57].  This  material  is  compatible  with 
existing  semiconductor  technologies  and  moreover  it  has  the  potential  to  respond 
much  faster  than  any  of  the  previously  studied  oxides  (it  responds  approximately 
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1000  times  faster  at  the  same  intensity).  Evidently,  this  may  have  important  ap¬ 
plications  in  high-speed  interconnects  and  in  all-optical  switching.  Photorefractive 
solitons  in  semiconductors  are  possible  in  spite  of  the  very  small  electrooptic  coeffi¬ 
cient  in  these  materials.  This  is  due  to  a  large  enhancement  of  the  space-charge  eld 
in  the  region  of  the  optical  beam  which  may  be  up  to  10  times  the  applied  eld.  The 
enhancement  results  from  a  resonance  which  occurs  when  the  dark  generation  rate 
of  electrons  is  equal  to  the  photo-generation  rate  of  holes.  This  effect  was  observed 
by  Chauvet  et  al.[57],[58].  An  alternative  way  to  achieve  such  desirable  fast  response 
times  is  also  possible  by  employing  photorefractive  solitons  of  high  intensity  [59]. 

1.1.2  Quasi-steady-state  solitons 

Unlike  steady-state  screening  solitons,  quasi-steady-state  solitons  can  only  exist  in 
transient,  i.e.  before  the  space  charge  eld  has  reached  equilibrium.  Interestingly 
enough,  quasi-steady-state  solitons  were  the  rst  to  be  predicted  [15], [16]  and  ob¬ 
served  in  biased  photorefractives  [19], [20].  These  solitons  can  only  be  observed 
during  a  time  window  which  may  be  seconds  or  minutes  depending  on  the  material 
and  the  absolute  beam  intensity.  The  theory  of  quasi-steady-state  solitons  in  biased 
photorefractives  was  rst  developed  by  Segev  et  al.[15],[16].  In  this  formalism,  the 
self-trapping  effect  arises  from  the  phase  coupling  among  the  spatial  spectral  com¬ 
ponents  of  an  optical  beam.  This  phase  coupling  mechanism  is  modelled  using  the 
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two- wave  mixing  response  function  [17] ,  [18]  under  external  bias  conditions.  The  evo¬ 
lution  equation  describing  these  solitons  is  then  obtained  by  truncating  to  rst  order 
the  Taylor  series  expansion  of  the  two-wave  mixing  response  function.  By  doing  so, 
the  end  result  in  turn  suggests  that  in  this  case  the  light  induced  nonlinear  index 
change  now  follows  the  spatial  derivatives  of  the  beam  intensity  distribution.  Along 
these  lines  several  results  have  been  predicted.  First,  the  solitons  in  this  case  (which 
can  be  either  Gaussian  or  hyperbolic-secant  like)  are  independent  of  the  absolute 
beam  intensity  as  long  as  the  light  intensity  is  well  above  the  dark  irradiance  of  the 
crystal  [15], [16].  Second,  these  solitons  are  only  possible  when  the  bias  eld  strength 
lies  within  a  certain  range  of  values.  These  properties  were  veri  ed  experimentally 
in  strontium  barium  niobate  [19], [20].  In  this  experiment,  stable  propagation  of  two 
dimensional  solitons  was  observed  which  was  also  independent  of  the  beam  inten¬ 
sity.  Subsequent  experimental  efforts  conducted  in  the  quasi-steady-state  regime 
have  also  demonstrated  one-dimensional  solitons  as  well  as  planar  dark  and  vortex 
solitons  [19], [20], [60].  On  the  theoretical  side,  the  stability  of  these  class  of  solitons 
has  been  investigated  [61]  and  the  dynamical  evolution  of  Gaussian  beams  in  biased 
photorefractives  has  been  solved  exactly  [62].  In  [62],  spatial  compression  and  self¬ 
bending  was  found  to  be  possible  when  the  bias  eld  exceeds  the  allowed  range  of 
values  necessary  to  establish  a  Gaussian  soliton.  This  behavior  is  in  agreement  with 
the  experimental  results  reported  thus  far  [19], [20]. 
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In  another  recent  study  it  has  been  also  shown  that  quasi-steady-state  solitons 
can  be  employed  to  write  permanent  waveguides  intended  for  other  beams  at 
less  photosensitive  wavelengths.  Even  though  these  waveguides  were  established 
at  microwatt  levels,  they  were  in  fact  able  to  guide  other  intense  beams  (in  the 
bulk  of  the  crystal)  at  power  levels  of  approximately  1  Watt  [31].  Note  that  no 
other  method  exists  today  that  is  capable  of  establishing  waveguides  in  the  bulk  of 
a  material.  Alternative  methods  currently  employed  in  Integrated  Optics  such  as 
sputtering,  liquid  phase  and  molecular  beam  epitaxy,  ion  implantation,  diffusion  and 
ion  exchange  are  extremely  involved  and  expensive  and  can  only  provide  waveguides 
at  the  surface  of  a  substrate. 

1.1.3  Photovoltaic  solitons 

Photovoltaic  solitons  are  possible  in  photorefractive  media  with  appreciable  photo¬ 
voltaic  coefficients  under  steady-state  conditions.  These  solitons  were  rst  predicted 
by  Valley  et  al.  [27]  and  observed  by  Taya  et  al.  [28].  Unlike  screening  solitons  or 
quasi-steady-state  solitons,  photovoltaic  solitons  do  not  require  an  external  bias. 
Instead,  in  this  case,  the  required  nonlinear  index  change  is  set  up  by  a  strong  pho¬ 
tovoltaic  current.  So  far,  the  theory  of  photovoltaic  solitons  has  been  developed  in 
one  transverse  dimension  [27]  and  by  appropriately  incorporating  the  photovoltaic 
nonlinearity  into  the  wave  equation.  By  doing  so,  explicit  soliton  solutions  (dark 
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and  bright)  have  been  obtained  that  resemble  those  in  Kerr  saturable  media  [63]. 
Prom  the  mathematical  point  of  view,  the  evolution  of  1-D  photovoltaic  spatial  soli- 
tons  is  governed  by  an  equation  very  similar  to  that  of  photorefractive  screening 
solitons  where  this  time  nonlinearity  constant  is  proportional  to  the  photovoltaic 
eld  constant  [27]. 

The  rst  experimental  observation  of  photovoltaic  solitons  was  carried  out  in 
lithium  niobate  by  Taya  et  al  [28].  This  crystal  exhibits  a  negative  nonlinear  index 
change  and  it  is  thus  suited  for  dark  solitons.  In  these  experiments  dark  and  gray 
solitons  have  been  observed  and  the  tensorial  nature  of  the  underlying  photovoltaic 
effect  also  became  evident.  In  a  subsequent  experiment,  Y-junction  waveguides  have 
been  formed  by  utilizing  higher-order  dark  soliton  splitting.  These  Y-junctions  were 
then  used  to  guide  and  split  other  more  intense  beams  at  higher  wavelengths  [29]. 

Recently,  2-D  photovoltaic  vortex  solitons  were  successfully  generated  [64].  This 
is  very  surprising  because  the  photovoltaic  current  in  lithium  niobate  is  inherently 
a  polar  current,  i.e.  the  electrons  are  accelerated  along  the  ferroelectric  axis  only. 
Nevertheless,  the  2-D  nature  of  the  optical  beam  enforces  a  two-dimensional  current 
and  in  turn  a  2-D  space-charge  eld  capable  of  supporting  2-D  spatial  solitons. 
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1.2  Incoherent  solitons 


In  a  recent  experimental  study,  Mitchell  et  al.  have  reported  for  the  rst  time  self¬ 
trapping  of  partially  spatially  incoherent  light  in  a  biased  SBN  crystal  [35], [36].  As 
previously  noted,  this  represents  an  important  discovery  not  only  in  the  area  of  pho- 
torefractives  but  in  the  entire  eld  of  Optics  in  general.  The  experiment  was  carried 
out  using  laser  light  after  a  rotating  diffuser.  After  the  diflfuser,  the  incoherent  light 
beam  (which  can  be  considered  to  be  quasi-thermal),  was  then  launched  in  a  biased 
SBN  specimen.  In  the  absence  of  any  external  bias,  the  input  beam  (30  micrometers 
in  diameter)  was  found  to  expand  to  102  microns  after  6  millimeters  of  propagation. 
Instead,  at  a  bias  voltage  of  550  Volts  the  incoherent  beam  remained  invariant  by 
coimteracting  diffraction  effects.  At  even  higher  voltages,  spatial  compression  was 
found  to  occur.  These  results  clearly  demonstrate  that  incoherent  self-trapping  ef¬ 
fects  were  present  and  to  some  extent  they  suggest  that  spatial  incoherent  solitons 
may  be  possible.  Of  all  the  processes  discussed  in  this  section,  this  latter  one  used 
to  be  the  least  understood. 
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Chapter  2 


Background 


In  a  recent  experiment  [35],  it  has  been  demonstrated  that  incoherent  light  can 
self-trap  itself  in  photorefractives.  More  speci  cally,  incoherent  hght  was  found 
to  self-trap  in  a  biased  Strontium  Barium  Niobate  (SBN:75)  photorefractive  crys¬ 
tal.  This  observation  was  made  possible  through  the  nonlinearity  associated  with 
photorefractive  screening  solitons.  This  nonlinearity  has  two  properties  which  are 
critical  for  the  self-trapping  of  incoherent  light  beams:  Finite  response  time  which 
can  be  much  longer  than  the  phase  uctuation  time  across  the  incoherent  beam;  sat¬ 
uration  which  is  essential  to  create  a  multimoded  waveguide.  Therefore,  to  further 
investigate  the  incoherent  self-trapping  process,  it  is  important  to  understand  the 
photorefractive  nonlinearity  and  coherent  screening  solitons.  In  this  chapter  we  will 
rst  give  a  brief  explanation  of  the  photorefractive  effect.  We  will  later  formulate 
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the  theory  of  coherent  photorefractive  solitons  and  discuss  the  properties  of  partially 
spatially  incoherent  beams. 

2.1  Photorefractive  nonlinearity 

The  optical  properties  of  photorefractive  (PR)  crystals  have  been  a  topic  of  consid¬ 
erable  interest  in  the  last  three  decades  or  so.  The  photorefractive  effect  was  rst 
discovered  accidentally  in  Lithium  Niobate  {LiNbOz)  in  1966  by  Ashkin  et  al.  [65]. 
Since  then,  it  has  been  observed  in  many  electro-optic  crystals  such  as  Barium  Ti- 
tanate  (BaTiOz),  Strontium  Barium  Niobate  (SBN,  SrxBai-xNb^Oz),  Potassium 
Niobate  (KNbOz),  Bismuth  Silicon  Oxide  (BSO,  Bzi2Si02o),  Bismuth  Titanium 
Oxide  (BTO,  BiuTiOzo),  and  several  other  materials  such  as  polymers.  Tradition¬ 
ally,  photorefractive  crystals  have  been  employed  in  various  applications  including 
real-time  holography,  phase  conjugation,  and  optical  data  storage  [17], [18].  How¬ 
ever,  the  prediction  of  optical  solitons  in  photorefractive  media  [15],  [16]  has  opened 
up  a  new  exciting  avenue  for  different  applications. 

In  general,  self-trapped  optical  beams  or  solitons  are  possible  when  the  process 
of  diffraction  is  exactly  balanced  by  nonhnear  self-focusing.  In  the  case  of  biased 
photorefractives,  the  nonlinearity  arises  from  the  dominant  drift  conduction  term. 
Currently  the  most  widely  accepted  model  describing  the  photorefractive  effect  is 
the  one  initially  suggested  by  Vinetskii  and  Kukhtarev  in  1975  [66], [67].  In  this 
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model,  a  PR  material  involves  both  donor  and  acceptor  centers.  All  donor  impurities 
are  assumed  to  be  identical  and  have  exactly  the  same  energy  level  somewhere  in 
the  middle  of  the  band  gap.  These  neutral  donor  impurities  No  can  be  ionized 
(resulting  to  electron  donation)  by  absorbing  photons  of  sufficient  energy.  In  turn 
the  ionized  ones  are  capable  of  recapturing  the  donated  photoelectrons.  On  the 
other  hand,  the  acceptor  impurities  Na  (whose  density  is  typically  much  smaller 
than  that  of  the  donors)  are  assumed  to  be  fully  ionized  at  all  times.  Therefore, 
they  only  participate  indirectly  in  the  PR  effect  by  providing  charge  neutrality  when 
no  optical  illumination  is  present.  In  the  absence  of  light,  the  ionized  donor  density 
iVp  is  taken  to  be  equal  to  the  acceptor  density  Na-  Nonuniform  illumination  of 
the  material  with  light  of  suitable  wavelength  ionizes  the  donors  and  generates  free 
carriers  (either  electrons  or  holes,  or  both).  The  free  electrons,  generated  in  the 
bright  regions  of  the  incident  optical  beam,  diffuse  and  drift  away.  As  a  result,  a 
charge  redistribution  between  illuminated  aud  dark  areas  is  induced.  This  space- 
charge  density  is  associated,  through  Gauss  s  law,  to  a  low  frequency  space-charge 
electric  eld  which,  in  turn,  causes  a  refractive  index  pattern  via  the  linear  electro¬ 
optic  Pockels  effect.  In  other  words,  non-uniform  illumination  of  the  photorefractive 
medium  results  in  a  non-uniform  change  in  the  refractive  index.  The  local  change 
in  the  refractive  index  is  given  by  A(l/n^)ij  =  VijkEk  where  Ek  are  space  charge 
electric  eld  components  and  ry  are  the  linear  electro-optic  coefficients  [17], [18]. 
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Even  though  the  electro-optic  tensor  involves  a  number  of  independent  coefficients, 
this  number  decreases  as  the  symmetry  increases. 


2.2  Steady-state  photorefractive  screening  solitons 

In  this  section,  we  develop  the  theory  of  steady-state  photorefractive  screening  soli¬ 
tons.  Under  steady-state  conditions,  the  evolution  equations  of  one  dimensional 
spatial  optical  beams  in  photorefractive  media  are  obtained.  As  an  example,  we 
will  consider  screening  solitons  in  a  SEN  photorefractive  crystal  where  most  of  the 
experiments  have  been  carried  out. 

2.2.1  Space-cheirge  eld 

To  study  the  propagation  of  1-D  planar  spatial  beams  in  this  material  system,  we 
consider  a  standard  con  guration  where  the  optical  beam  propagates  along  the  z  axis 
and  it  is  allowed  to  diffract  only  along  the  x  coordinate.  The  optical  c-axis  of  the  SEN 
crystal  is  assumed  to  be  oriented  in  the  x  direction.  Furthermore,  an  external  electric 
eld  is  applied  along  the  same  coordinate.  In  this  case,  the  perturbed  extraordinary 
refractive  index  along  the  optical  c  axis  is  given  by  —  n^VasEac  where 

is  the  electro-optic  coefficient  involved,  is  the  unperturbed  extraordinary  index 
of  refraction,  and  Esc  is  the  total  electric  eld.  The  induced  space-charge  eld  E^c 
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can  be  obtained  from  the  Kuktarev-Vinetskii  equations  [66], [67], [68].  For  a  one 
dimensional  con  guration  and  under  steady-state  conditions  these  equations  take 
the  following  form: 


lB.Nin  =  Si{I  +  h)  [Nd-N^) 

(2.1) 

J  =  constant 

(2.2) 

(  ^  KBTdn\ 

J-eii\nEac+  g 

(2.3) 

(2.4) 

The  rst  equation  is  the  charge  recombination  equation  and  the  second  one  is  the 
continuity  equation.  The  third  one  shows  that  the  current  density  J  is  the  sum  of 
the  drift  and  diffusion  terms.  Moreover,  I  =  I  {x,  z)  is  the  power  density  of  the 
optical  beam  and  is  proportional  to  the  square  of  the  absolute  magnitude  optical 
eld  envelope,  is  the  carrier  recombination  rate,  Sj  is  the  photo-excitation  cross 
section,  e  is  the  electron  charge,  fi  is  the  electron  mobility,  n  is  the  free  electron 
density,  and  Kb  is  Boltzmann  constant;  T  is  the  absolute  temperature,  GoGr  is  the 
total  permittivity,  Np  is  the  donor  concentration,  Na  is  the  trap  density,  is  the 
ionized  donor  density.  Id  is  the  dark  irradiance  of  the  crystal.  Furthermore,  let  us 
assume  that  the  optical  intensity  can  remain  constant  at  the  tails  far  away  from  the 
beam,  i.e.  I{x  ±(Xi,z)  =  loo-  It  can  then  be  shown  that  the  space-charge  eld 
is  also  constant  at  the  tails,  i.e.,  >  ±00,2:)  =  Eq.  When  the  beam  is  well 
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con  ned  within  the  crystal  s  width  W,  the  constant  value  Eo  is  approximately  equal 
to  V/W  where  V  is  the  .external  bias  voltage.  As  usual,  we  assume  that  the  donor 
concentration  and  the  trap  density  is  much  bigger  than  the  free  electron  density 
and  that  the  ionized  donor  density  is  approximately  the  same  as  the  trap  density. 
Under  these  assmnptions,  the  steady-state  space  charge  eld  can  be  obtained  from 
the  above  Kuktarev-Vinetskii  transport  equations  and  given  by  [24], 


E.,. 


p  (/oo  +  h) 

KbT  GpGr 
e  cNa 


1  + 


1  + 


Go^r  dEsc 
cNa  dx 

Go^r  9Ea( 


KbT 


di 


e  {I  +  Id)  dx 


(2.5) 


-1 


dx"^ 


eNa  dx 

When  the  bias  voltage  V  is  high,  the  space-charge  eld  Esc  is  dominated  by  drift 
and  typically  in  photorefractive  crystals  |(Go£r  l&KjC)  dEsc/dx\  1.  In  that  case, 
the  space-charge  eld  is  approximately  given  by  [23], [24] 


Esc «  E, 


(/qO  +  Id) 
(/  -I-  Id) 


(2.6) 


2.2.2  Evolution  equation 

Starting  from  Maxwell  s  equations  and  by  employing  standard  procedures,  we  nd 
that  the  envelope  of  the  optical  eld  obeys  the  following  differential  equation  [23] ,  [24] : 

(2.7) 


+  — —  -  Y  [nlrasEsc)  0  =  0. 


dz  2k  dx^  2 

In  this  equation,  0  is  a  slowly  varying  envelope.  In  turn,  the  optical  eld 

£  =  x0  (x,  z)  exp  (ikz) 
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is  expressed  in  terms  of  this  envelope,  k  =  (27r/Ao)ne  and  Aq  is  the  free-space 
wavelength.  Considering  only  the  dominant  drift  process,  the  evolution  equation  for 
the  envelope  </>  can  be  obtained  by  substituting  Eq.  2.6  into  Eq.  2.7.  For  simplicity, 
this  equation  can  be  further  simpli  ed  by  adopting  normalized  coordinates  and 
quantities,  i.e.  ^  =  z/  (kx^),  s  =  x/xq,  ^  =  {2r}Jd/nef'^  U  where  xq  is  an  arbitrary 
spatial  scale  and  the  beam  intensity  has  been  scaled  with  respect  to  dark  irradiance. 
In  that  case,  the  normalized  evolution  equation  is  given  by 


.dU  l&^U 
+ 


;0(l  +  p) 


rC/  =  0 


(2.8) 


di  '  2  38-^  I  ^  |t7|2 

where  p  =  loo/h  and  the  nonlinearity  constant  P  is  proportional  to  effective  elec¬ 
trooptic  coefficient  of  the  photorefractive  crystal  and  the  external  bias  eld,  i.e. 
P  =  {koxof  n‘^r33Eo/2.  Note  that  P  can  be  positive  or  negative  depending  on  the 
polarity  of  the  bias  eld  Eq.  As  it  has  been  shown  [24],  bright,  dark  and  gray 
one-dimensional  spatial  screening  solitons  can  be  possible  in  biased  photorefractive 
crystals.  Moreover  it  was  shown  that  for  a  given  physical  system,  these  states  de¬ 
pend  only  on  two  parameters,  namely  the  strength  of  the  applied  bias  eld  Eq  as 
well  as  the  so-called  intensity  ratio  r  or  p,  i.e.  the  ratio  of  the  optical  peak  intensity 
with  respect  to  the  dark  irradiance  of  the  crystal.  The  symbol  r  represents  the 
maximum  intensity  ratio  of  a  bright  soliton  whereas  p  stands  for  the  maximum  dark 
intensity  ratio.  Bright  and  dark  soliton  beam  pro  les  are  shown  in  Fig.2.1(a)  and 
(b).  The  dependence  of  the  normalized  spatial  intensity  FWHM  of  these  beams 
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respect  to  r  or  p  in  units  of  x/xq  is  given  Fig.  Fig.  2. 1(c)  and  (d).  Here  Xo  is  a  nor¬ 
malization  parameter  and  given  by  xq  =  l/(A:on.gyT33F/o/2).  It  is  important  to  note 
that,  bright  soliton  states  require  positive  bias,  Eq  >  0,  whereas  the  dark  and  gray 
soliton  states  are  possible  only  when  the  bias  eld  is  negative  or  the  photorefractive 
material  is  reverse  biased,  Eq  <  0- 

Photorefractive  solitons  require  very  low  powers,  as  low  as  1  fiWatt,  because  the 
nonlinear  refractive  index  change  depends  on  the  I /Id  ratio  rather  than  the  optical 
intensity  I.  In  this  way,  even  if  the  total  power  is  very  low,  a  very  strong  nonlin¬ 
earity  can  be  induced.  Moreover,  the  nonlinearity  can  be  adjusted  by  changing  the 
external  bias  voltage  thus  introducing  an  extra  free  parameter  to  vary  the  size  of 
the  soliton  beam.  In  addition,  the  nonlinear  response  of  photorefractive  materials  is 
wavelength  dependent.  This  means  that,  one  can  use  a  soliton  induced  waveguide 
to  control  powerful  beams  at  wavelengths  at  which  the  material  is  less  photosensi¬ 
tive.  Moreover,  the  response  time  of  the  photorefractive  material  can  change  from 
nanoseconds  to  hours  depending  on  the  optical  intensity.  For  example,  the  response 
time  is  in  the  order  of  seconds  for  a  typical  soliton  formation  in  a  SEN  crystal.  In 
fact,  this  property  is  one  of  the  requirements  for  incoherent  soliton  formation. 
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Figure  2. 1 :  Field  profiles  of  (a)  bright  and  (b)  dark  screening  spatial  solitons. 
Normalized  intensity  full-width  half-maximum  versus  intensity  ratio  for  (c)  bright  and 
(b)  dark  screening  solitons. 


2.3  Spatially  incoherent  light  £ind  coherence  length 

For  a  coherent  beam,  the  phases  at  all  points  of  the  beam  varies  in  unison  with 
time.  On  the  other  hand,  the  phases  at  different  points  across  an  incoherent  beam 
vary  in  a  statistical  manner.  This  randomly  changing  phase  distribution  creates  an 
intensity  pattern  known  as  speckle.  Speckle  size  is  an  indication  of  the  degree  of 
coherence  for  a  spatially  partially  incoherent  beam.  The  diffraction  characteristics 
of  a  beam  are  related  to  the  speckle  size  or  the  degree  of  coherence:  the  smaller 
the  speckles,  the  larger  the  diffraction.  In  other  words,  beams  of  higher  coherence 
diffract  less.  When  the  beam  is  fully  coherent,  every  point  of  the  beam  is  well 
correlated  with  all  other  points.  Therefore,  the  coherence  length  across  a  coherent 
beam  is  in  nite.  Each  point  of  a  fully  incoherent  beam  ,  however,  correlates  only 
with  itself  and  it  is  totally  uncorrelated  with  other  and  thus  the  correlation  length 
is  zero.  Between  these  two  extreme  cases,  is  the  class  of  partially  incoherent  beams. 
Most  of  the  sources  in  nature  fall  in  this  category.  In  this  case,  each  point  across 
the  beam  correlates  with  other  points.  This  nite  correlation  distance  is  called 
coherence  length. 

The  coherence  of  a  source  can  be  described  by  using  a  time  averaged  temporal 
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correlation  between  the  electric  eld  at  different  points  across  the  beam.  This  prop¬ 
erty  is  given  by  a  normalized  spatial  complex  coherence  function  [69], [70]  |u(xi,X2) 


A*12(®1)  ^2) 


{Eixr)E*{x2)) 


(2.9) 


^I{Xi)I{x2) 

where  (.)  denotes  the  time  average  and  I  is  the  optical  beam  intensity.  In  general, 
the  complex  coherence  function  can  be  an  arbitrary  function  of  both  xi  and  X2. 
Thus,  it  is  often  more  convenient  to  describe  the  coherence  properties  of  an  optical 
beam  through  a  coherence  length  4  which  is  de  ned  as  follows: 


C» 

,(a;)=  J  \fj,i^{x,S)f  dS 


(2.10) 


To  understand  the  complex  coherence  function  and  its  relation  to  optical  intensity, 
consider  two  point  sources,  let  say  1  and  2.  The  optical  eld,  E,  at  an  observation 
point  will  be  superposition  of  the  elds,  Ei  and  E2,  from  these  two  points.  In 
this  case  the  total  eld  at  the  observation  point  will  be  given  by  E  =  Ei  +  E2. 
Furthermore,  if  we  assume  that  the  eld  is  stationary,  i.e.  the  coherence  function 
only  depends  on  the  separation  between  these  two  points  r  =  r2  —  ri  but  not 
their  placement,  Eq.2.9  can  be  written  as  fJ-i2{r)  =  |a^i2  (''')[  ^^  ('’■)]■  Here, 

0  <  |Ati2('r)|  <  1  and  </)i2('r)  =  arg  [yu^a  ('^)]  is  a  time  averaged  phase  difference 
between  two  elds.  Therefore  the  time  averaged  total  intensity  can  be  given  as 


It  =  ^{EE*)=h  +  l2  +  2{hl2f^'R-e[fi,2{r)]  (2.11) 

=  -I-  72  -f  2  |Mi2  (■’’)!  [*^12  Wl 
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where  Jiand  I 2  are  the  intensities  that  would  be  produced  at  the  observation  point  by 
two  elds  independently.  Eq.  2.11  shows  that,  depending  on  the  degree  of  coherence 
between  two  elds,  different  intensity  patterns  are  formed.  For  example  if  the  beams 
are  completely  coherent  among  each  other,  i.e.|//i2  (t)|  =  1,  a  fringe  pattern  with  a 
maximum  contrast  will  be  observed.  When  the  beames  are  completely  imcorrelated, 
i.e.  I/X12  =  0)  no  fringe  pattern  will  be  observed.  This  latter  case  corresponds  to 

complete  incoherence,  i.e.  there  is  no  de  nite  phase  relationship  between  the  waves 
arriving  at  the  observation  point  from  two  different  point  sources.  When  these 
beams  are  partially  correlated,  0  <  \^ll2  (r)!  <  1,  the  contrast  of  the  fringe  pattern 
will  change  accordingly  .  This  latter  case  corresponds  to  partially  incoherent  beams 
[69],[70]. 

The  time  averaged  intensity  of  a  partially  spatially  incoherent  beam  is  different 
than  its  instantaneous  intensity.  If  we  observe  the  intensity  of  an  incoherent  light 
at  any  time,  we  will  see  a  speckled  intensity  pattern  as  shown  in  Fig.2.2.  What 
our  eyes  or  a  slow  responding  material  see,  in  fact,  will  be  an  average  of  these 
randomly  changing  speckled  patterns.  If  the  averaging  time  is  long  enough,  the 
resulting  intensity  will  be  smooth  as  shown  in  Fig.2.2.  In  this  latter  case,  the  third 
term  (interference  term)  in  Eq.  2.11  averages  out  to  zero  since  the  nonlinearity  re¬ 
sponds  slower  than  the  rate  of  random  phase  changes.  Therefore,  the  total  intensity 
that  nonlinearity  sees  will  be  the  summation  of  contributing  intensities  from  each 
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individual  point,  i.e.  It  ^h  +  h- 
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Figure  2.2:  Illustration  of  the  temporal  avaraging  of  many  instantaneous  speckled 
patterns  due  to  the  noninstantaneous  response  of  nonlinearity. 
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Chapter  3 

Theory  of  incoherent  becun 
propagation  in  nonlinear  media 

In  this  chapter,  the  coherent  density  approach,  describing  propagation  and  self- 
focusing  of  partially  spatially-incoherent  light  beams  in  nonlinear  media  is  devel¬ 
oped.  It  is  shown  that  this  process  is  effectively  governed  by  an  in  nite  set  of  coupled 
nonlinear  Schrddinger-like  equations  provided  that  they  are  initially  appropriately 
weighted  with  respect  to  the  incoherent  angular  power  spectrum  of  the  source.  The 
particular  case  of  spatially  partially  incoherent  beam  propagation  in  biased  photore- 
fractive  media  is  considered  in  detail.  Numerical  simulations  indicate  that  spatial 
compression  as  well  as  self-trapped  states  are  possible  under  appropriate  conditions. 
Our  results  are  in  good  agreement  with  recent  experimental  observations. 
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3.1  Introduction 


Optical  self-focusing  has  been  a  subject  of  considerable  interest  in  the  last  three 
decades  or  so.  Over  the  years,  this  process  has  been  systematically  investigated  in 
all  states  of  matter  with  the  aid  of  laser  sources  [12].  Thus  far,  several  physical 
systems  have  been  identi  ed  that  can  lead  to  optical  self-trapping.  These  include 
for  example,  or  Kerr-like  media  [10],  [11],  the  x^^^  family  of  materials  [71], [72] 
as  well  as  the  class  of  photorefractives  [19]- [25], [73].  Highly  relevant  to  this  topic 
is  of  course  the  very  existence  of  optical  spatial  solitons  [8], [9].  These  latter  en¬ 
tities  can  occur  provided  that  diffraction  effects  are  exactly  balanced  by  optical 
self-trapping.  At  this  point,  it  may  be  fair  to  say  that  both  coherent  self-focusing  as 
well  as  the  coherent  excitation  of  solitons  are  by  now  in  principle  well  understood. 
In  a  very  recent  study  however,  a  successful  observation  of  incoherent  self-trapping 
has  been  reported  for  the  rst  time  [35].  In  this  experiment,  partially  spatially  in¬ 
coherent  light  was  found  to  self-trap  in  a  biased  strontium  barium  niobate  (SEN) 
photorefractive  crystal.  Apart  from  its  possible  scienti  c  and  technological  impli¬ 
cations,  this  observation  [35]  in  turn  poses  a  new  fundamental  challenge.  More 
sped  cally,  a  theory  of  incoherent  self- focusing  and  possibly  of  incoherent  solitons 
now  needs  to  be  developed.  Unlike  their  coherent  beam  counterparts,  for  which  the 
phase  at  all  points  varies  in  unison  with  time,  the  phases  at  different  points  across 
an  incoherent  beam  vary  in  an  uncorrelated  manner  [69], [70].  This  introduces  an 
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important  new  element  in  the  nonlinear  theory  of  self-focusing.  Even  though  the 
properties  of  speckle-inhomogeneous  elds  have  been  considered  in  the  past  in  con¬ 
nection  with  optical  phase  conjugation  [74],  the  propagation  behavior  of  partially 
incoherent  (multimode)  beams  in  nonlinear  media  has  not  yet  been  explored. 

In  this  chapter,  a  theory  of  incoherent  self-focusing  in  non-instantaneous  non¬ 
linear  media,  and  in  particular  in  biased  photorefractive  crystals  is  presented.  We 
consider  nonlinear  media  with  response  times  much  bigger  than  the  characteristic 
phase  actuation  time  across  the  optical  beam  which  in  turn  will  only  experience 
the  time-averaged  intensity.  In  the  case  of  stationary  spatial  source  actuations,  we 
show  that  this  process  can  be  effectively  described  by  means  of  an  in  nite  set  of 
coupled  nonlinear  Schrodinger-hke  equations  provided  that  they  are  appropriately 
weighted  with  respect  to  the  source  incoherent  angular  power  spectrum.  Our  nu¬ 
merical  computations  demonstrate  that  spatial  compression  as  well  as  self-trapped 
states  are  possible  under  appropriate  conditions.  In  these  cases,  self-trapping  can  be 
intuitively  understood  as  fusion  of  multi-particles.  Pertinent  examples  are  provided 
to  further  elucidate  this  behavior.  A  possible  soliton  solution  to  this  system  is  also 
presented. 
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3.2  Coherent  density  approach 


To  start,  let  us  consider  for  example  a  SBN  crystal  with  its  optical  c-axis  oriented 
in  the  x  direction.  Let  us  also  assume  that  the  optical  beam  propagates  along  the 
z  axis  and  it  is  allowed  to  diffract  only  along  the  x  direction.  For  simplicity,  we 
limit  our  analysis  to  one  transverse  dimension  (x)  and  we  assume  uniformity  in 
y.  Furthermore,  an  external  bias  electric  eld  is  applied  along  x  (i.e.,  the  c-axis), 
in  which  case  the  perturbed  extraordinary  refractive  index  is  given  by  [23], [24], [59] 
(n'  )2  =  ng  -  njrssEsc,  where  is  the  unperturbed  index  of  refraction,  is  the 
electrooptic  coefficient  involved,  and  Esc  is  the  static  space  charge  eld  in  this  pho- 
torefractive  crystal.  Under  strong  bias  conditions  (0.2-4  kV/cm)  and  for  relatively 
broad  bright-like  beam  con  gurations  (  25Ao/ne),  the  steady-state  space  charge 
eld  is  approximately  given  by  [23],  [24], [59]  Esc  =  EoId[Id  +  ^  ^  shown 

in  Chapter  2.  Here  I  =  I{x,z)  is  the  power  density  of  the  optical  beam.  Id  is  the 
so-called  dark  irradiance  of  the  crystal  and  Eq  is  the  value  of  the  space  charge  eld 
at  X  — >  ±oo.  If  the  spatial  extent  of  the  optical  wave  involved  is  much  less  than  the 
a: —width  W  of  the  crystal,  then  for  a  constant  bias  voltage  V,  Eo  is  approximately 
given  by  ±V/W  [24]. 

Next,  it  is  important  to  consider  the  diffraction  behavior  of  this  incoherent  beam. 
Let  us  assume,  as  in  the  experiment  reported  [35],  that  the  incoherent  wavefront 
results  from  a  quasi-thermal  quasi-monochromatic  source  such  as  for  example  laser 
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light  after  passing  through  a  rotating  diffuser  [35], [70], [75].  Let  the  x  polarized 
electric  eld  component  of  the  optical  wave  be  expressed  in  terms  of  a  slowly  varying 
envelope  4>{x,z)  i.e.  E{x,z)  =  x(l>{x,z)ex-p{ikz)  where  k  =  (27r/Ao)ne.  In  the 
paraxial  limit  {k^^/k  <  1)  and  in  the  linear  regime,  the  envelope  (f)  is  known  to 
evolve  according  to 

<l>{x,  =  ^  I  dkMkx)  exp  {i[k,x  -  {kl/2k)z]}  (3.1) 

—  CO 

where  is  the  Fourier  transform  of  the  eld  right  at  the  input,  i.e.  at  z  =  0. 

Let  the  optical  eld  at  the  origin  also  be  written  as  <j){x,z  =  0)  =  f{x)<po{x) 
where  f{x)  is  a  spatial  modulation  function  and  (f>o{x)  is  the  eld  before  modu¬ 
lation  which  implicitly  contains  all  the  spatial  statistical  properties  of  the  source. 
If  the  source  uctuations  obey  a  stationary  random  process  [76],  then  the  spatial 
statistical  autocorrelation  function  of  (I)q{x)  is  given  by  (^(l>o{x)4>l{x)'^  =  R{x  -  x). 
In  turn,  the  autocorrelation  function  of  the  source  spectrum  can  be  obtained,  i.e. 

=  ‘2'K6{k^-k'^)G{k'^)  where  ^o{ka:)  and  G{kJ  are  the  Fourier  trans¬ 
forms  of  ^o(x)  and  R{x)  respectively  and  6{x)  is  a  delta  function.  Physically  the  real 
function  [76]  G{k^)  represents  the  incoherent  angular  power  spectrum  of  the  source. 
Using  the  frequency  convolution  theorem  and  the  fact  that  ^{k^)  is  the  Fourier 
spectrum  of  the  product  f{x)<f)Q{x),  one  can  readily  derive  the  following  result: 

OO 

(Hk^mk'j)  =  (l/27r)  I  dvG{ii)F{k,  -  v)F*ik',  -  v)  (3-2) 

— cx> 
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where  F{kx)  is  the  Fourier  transform  of  the  spatial  modulation  function  f{x).  Keep¬ 
ing  in  mind  that  the  intensity  of  this  wave  is  given  by  I{x,z)  =  and  by 

employing  Eq.3.1  and  3.2  we  nally  obtain 

oo  oo  ^ 

I{x,z)=  j  dBGN{6)  ^  J  dkxF{kx)exp[ikx{x -9z)]exp[-ikl{z/2k)]  (3.3) 

—  OO  -OO 

where  Gn{6)  in  Eq.3.3  has  been  normalized  for  convenience,  and  9  =  k^/k  represents 
an  angle  (in  radians)  with  respect  to  the  2  axis.  At  this  point  it  may  be  useful 
to  make  few  remarks.  First,  the  quantity  in  the  brackets  of  Eq.3.3  represents  in 
fact  the  intensity  pro  le  resulting  from  an  otherwise  coherent  beam  when  its  initial 
eld  pro  le  is  f{x).  Furthermore,  this  coherent  component  propagates  at  an 
angle  9  with  respect  to  the  2-axis  by  obeying  the  paraxial  equation  of  diffraction: 
i{Ux  +  9Ux)  +  (l/2A:)C4x  =  0  where  C4  =  dU/dz  etc.  In  essence,  Eq.3.3  leads 
to  the  following  important  conclusion:  the  diffraction  behavior  of  an  incoherent 
beam  can  be  effectively  described  by  the  sum  of  the  intensity  contributions  from  all 
its  coherent  components  provided  that  their  eld  pro  les  at  the  origin  have  been 
appropriately  scaled  with  respect  to  the  incoherent  angular  power  spectrum  Gn{9), 
i.e.  U{x,z  =  0)  =  G^^^{9)f{x).  As  one  may  anticipate,  in  the  limit  Gm{9)  — > 
6{9),  the  result  of  Eq.3.3  reduces  to  that  of  the  coherent  case.  Note,  that  similar 
arguments  have  been  previously  employed  in  connection  with  incoherent  imaging 
[77],  theory  of  speckle-inhomogeneous  elds  [74],  and  incoherent  wave  propagation 
in  dispersive  media  [78]. 
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Up  to  this  point  our  treatment  is  quite  general.  On  the  other  hand,  when  an 
incoherent  optical  beam  propagates  in  a  slowly-responding  nonlinear  medium,  one 
should  also  expect  that  each  of  these  coherent  components  or  quasi-particles  will  be 
in  uenced  by  the  nonlinearity  involved.  In  turn,  an  intensity-dependent  nonlinearity 
will  follow  the  incoherent  (intensity)  superposition  of  all  these  components.  In 
the  particular  case  of  a  biased  photorefaractive  crystal,  the  nonlinear  index  change 
An  =  {nl/2)  r^^Esc  can  then  be  readily  incorporated  into  the  underlying  equations 
of  motion  by  following  the  procedure  of  [23] ,  [24] ,  [59] .  By  doing  so  and  by  discretizing 
the  diffraction  integral  of  Eq.3.3  (i.e.  6  j A6>),  then  under  steady-state  conditions 
we  are  nally  led  to  the  following  in  nite  set  of  coupled  nonlinear  Schrbdinger-like 
equations: 


(dUj 


1  ko 


2k  dx^ 


nlrssEo 


Uj 


1  +  I{x,z) 


=  0 


(3.4) 


where  I(x,z)  is  the  intensity  pro  le  of  the  incoherent  beam  which  is  given  by 


OO 

7(x,z)=  (3.5) 

j  =  — OO 

and  the  discrete  index  j  =  0,  ±1,  ±2, ...  .  Effectively,  Eqs.3.4  and  3.5  describe  the 
process  of  incoherent  self-focusing  in  biased  photorefractives  in  the  limit  — >  0. 

Moreover,  in  these  equations  each  coherent  fragment  has  been  scaled  with  respect 
to  Id,  and  at  the  origin  we  have  assumed  that  Uj{x,  0)  oc  {j A9)  f  {x) . 
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3.3  Incoherent  bright  photorefractive  solitons 


As  an  example  let  us  consider  a  biased  SBN:75  crystal.  Here,  the  parameters 
used  will  be  very  similar  to  those  reported  in  [35].  In  particular,  let  =  1022 
pm/V,  Ao  =  488  nm,  =  2.3,  W  =  6  mm  and  let  the  length  of  propagation 

be  6  mm.  Let  us  also  assume  that  the  input  intensity  pro  le  as  well  as  the  in¬ 
coherent  angular  power  spectrum  are  Gaussian,  i.e.  f{x)  =  exp  (— x^/2a:Q)  and 
^  exp  In  our  numerical  computations  we  use  201  com¬ 

ponents  (—100  <  j  <  100)  equidistantly  spanning  the  range  ±2.5  9o.  Prom  our 
previous  discussion,  the  eld  of  each  coherent  component  at  the  origin  is  set  to  be 
Uj{x, 0)  =  exp  (iA0)^  /20o]  where  r  is  an  appropriate  constant  which  is 

related  to  the  maximmn  intensity  ratio  rr  (with  respect  to  Id)  of  the  input  incoher¬ 
ent  beam.  More  speci  cally  ry  =  f  ^exp  [—  (jA0/0o)^]-  Equations  3.4  and  3.5  are 
solved  by  means  of  standard  split-step  Fourier  methods.  The  accuracy  of  our  results 
was  then  checked  against  the  conservation  laws  of  Eqs.  3.4  and  3.5  and  by  increasing 
the  number  of  coherent  components.  As  a  rst  example  let  us  consider  an  incoherent 
Gaussian  beam  with  Xo  =  18  jum,  in  which  case  its  input  intensity  FWHM  is  30 
nm.  Moreover  let  the  width  of  the  source  angular  power  spectrum  be  9q  =  9.56 
milliradians  or  0.548°.  In  this  case,  the  beam  linearly  diffracts  to  a  FWHM  of 
102  /um  after  6  mm  of  propagation  as  shown  in  Figure  3.1.  Note  that  if  this  beam 
was  spatially  coherent,  it  would  have  diffracted  to  35.4  //m  after  6mm,  as  has  been 
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observed  in  [35].  Once  the  crystal  is  appropriately  biased,  self-trapping  effects  start 
to  emerge.  Fig.3.2(a)  depicts  the  intensity  evolution  of  this  beam  (xq  =  18  /xm, 
9o  =  0.548°)  when  the  applied  voltage  is  400  V  and  rr  =  3.  After  6  mm,  the  beam 
has  developed  a  rectangular-like  pro  le  with  a  FWHM  of  ~  34.7  nm  as  shown  in 
Fig.3.2(b).  Evidently  self-focusing  played  an  important  role  in  this  example  even 
though  it  was  not  enough  to  balance  diffraction  effects.  Figure  3.3  on  the  other 
hand  shows  what  would  happen  if  the  bias  is  increased  to  550  V .  In  this  case,  the 
incoherent  beam  propagates  almost  undistorted  and  it  behaves  like  a  quasi-soliton. 
In  other  words,  all  the  coherent  components  or  quasi-particles  have  appropriately 
fused  together  thus  producing  a  stationary  beam.  Figs.3.4  and  3.5  show  intensity 
evolution  of  two  coherent  fragments  of  this  self-trapped  beam  at  ^  »  0.27°  and 
9  «  0.54°  respectively.  At  an  even  higher  voltage,  V  =  1000  V,  the  beam  starts  to 
exhibit  considerable  spatial  compression.  Figure  3.6  shows  that  cycles  of  compres¬ 
sion  and  expansion  are  now  possible  during  propagation.  At  the  output,  i.e.  z  =  6 
mm,  its  intensity  FWHM  is  ~  18  //m.  Note  that  behavior  of  this  sort  (including  the 
rectangular-like  beam  features  of  Figs. 3. 2  and  3.6)  is  consistent  with  previous  ex¬ 
perimental  observations  (In  Reference  [35]  only  the  self-trapped  case  was  reported. 
Nevertheless,  rectangular-like  beam  behavior  at  low  bias  and  beam  compression  at 
high  bias  voltages  has  been  observed  as  well).  Moreover,  had  this  beam  (30  /xm 
FWHM)  been  fully  coherent,  it  would  have  disintegrated  as  a  result  of  speckle  noise 
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instabilities.  Next,  let  us  consider  what  will  happen  at  higher  intensity  ratios  tt- 
Figure  3.7  shows  the  propagation  of  an  incoherent  Gaussian  beam  when  V  =  550 
V,  9o  =  0.548°  and  =  40.  At  z  =  6  mm,  the  beam  now  expands  from  30  fim 
FWHM  to  87  fim.  As  in  the  coherent  limit  [23], [24], [59],  this  expansion  is  attributed 
to  the  saturation  of  the  photorefractive  nonlinearity  when  ^  1.  In  our  simula¬ 
tions  the  angular  width  Oo  was  also  found  to  play  an  important  role.  Fig.3.8  depicts 
the  propagation  of  such  an  incoherent  Gaussian  beam  (30  fim  FWHM  at  the  input) 
when  V  =  550  V,  =  3  and  =  0.8°.  As  the  gure  shows,  the  beam  expands 
in  a  rectangular-like  fashion  to  a  FWHM  of  51.3  fim  after  6  mm  of  propagation. 
Clearly,  for  =  0.8°  a  higher  bias  voltage  is  required  to  overcome  diffraction  effects 
and  indeed  at  ~  950  V  self-trapping  is  reestabhshed. 


3.4  Solitary  wave  solution 

We  would  also  like  to  point  out  that  the  in  nite  system  of  coupled  nonlinear  partial 
differential  equations  of  Eq.3.4  does  in  fact  admit  solitary  wave  solutions.  To  obtain 
such  a  solution  we  write  Uj{x,z)  =  Uj{x,  z)  exp  {kz/2)  —  jAOkx^^  and 

Uj  =  ry^Q{x)  exp{iiJ,z),  where  Q{x)  is  a  normalized  function,  i.e.  0  <  Q{x)  <  1. 


Direct  substitution  of  these  latter  forms  in  Eq.3.4  leads  to  the  following  ordinary 


0.15H 


Figure  3.5:  Intensity  evolution  of  20th  component  \U2o(x,  6o~0.54°)l^  of  the  quasi- 
soliton  shown  in  Fig.3.3. 
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differential  equation 


^-2^kQ-{kV,r,A)j:^=0  (3.6) 

which  is  known  to  allow  bright  solitons  [23], [24], [5 9]  when  >  0  and 

lj,=  -  (^knlr33Eo/2rT)  ln(l  +  rr)  ■ 

In  Eq.3.6,  ry  is  the  total  intensity  ratio,  i.e.  rr  =  We  would  like  to  emphasize 

j 

however  that  physically  this  solitary  wave  solution  has  a  limited  range  of  applica¬ 
bility.  More  speci  cally  the  employed  Uj  -  Uj  transformation  implies  in  reality  that 
all  the  coherent  components  propagate  in  parallel  along  z.  Thus  a  solitary  solution 
of  this  sort  may  only  be  applicable  when  6o  is  quite  small  in  which  case  all  quasi¬ 
particles  may  tend  to  propagate  almost  in  parallel  after  fusion.  In  fact,  the  above 
mentioned  solution  represents  a  generalization  of  the  so-called  incoherent  coupled 
photorefractive  soliton  pairs  previously  discussed  in  the  literature  [50] ,  [34] . 

3.5  Conclusions 

In  conclusion,  a  theory  of  incoherent  self-trapping  in  biased  photorefractives  has 
been  developed.  It  has  been  shown  that  this  process  can  be  effectively  described  by 
an  in  nite  set  of  coupled  nonlinear  Schrbdinger-like  equations  provided  they  have 
been  initially  weighted  with  respect  to  incoherent  angular  power  spectrum  of  the 
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source.  Relevant  examples  have  been  provided.  Our  numerical  computations  were 
in  close  agreement  with  recent  experimental  data.  We  wish  to  emphasize  that  our 
theoretical  approach  applies  not  only  to  photorefractive  media  but  also  to  any  other 
nonlinear  material  whose  temporal  response  time  is  much  bigger  than  the  phase 
uctuation  time  of  an  optical  beam. 
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Chapter  4 


Incoherent  spatial  solitons  in 
saturable  nonlinear  media  of  the 
logarithmic  type 


By  using  the  coherent  density  approach,  it  is  shown  that  partially  incoherent  spatial 
solitons  are  possible  in  saturable  nonlinear  media.  This  is  demonstrated  by  means 
of  an  exact  Gaussian  solution  in  the  case  where  the  saturable  nonlinearity  is  of  the 
logarithmic  type.  The  conditions  necessary  to  establish  these  soliton  states  as  well 
as  their  associated  characteristics  are  discussed  in  detail.  Pertinent  examples  are 
provided  to  further  elucidate  their  behavior. 


53 


4.1  Introduction 


In  the  last  chapter,  a  theory  was  developed  in  order  to  explain  the  observed  inco¬ 
herent  self- focusing  behavior  in  biased  photorefractive  media  [79] .  More  speci  cally, 
it  was  found  that  this  process  can  be  effectively  described  by  an  in  nite  set  of 
nonlinear  Schr6dinger-like  equations  provided  that  they  are  initially  appropriately 
weighted  with  respect  to  the  incoherent  angular  power  spectrum  of  the  source.  The 
numerical  simulations  resulting  from  this  theoretical  approach  [79]  were  found  to 
be  in  good  agreement  with  the  experimental  observations  of  Ref.  [35].  It  is  also 
important  to  note  that  these  experimental  [35]  and  theoretical  [79]  investigations 
have  demonstrated  that  a  stationary  incoherent  beam  or  a  quasi-soliton  could  exist 
under  appropriate  conditions.  Nevertheless,  there  is  no  solid  evidence  (in  terms  of 
an  analytical  solution)  as  to  whether  such  incoherent  sohtons  exist  or  not. 

In  this  chapter,  we  show  that  partially  incoherent  spatial  sohtons  are  indeed 
possible  in  saturable  nonlinear  media.  This  is  demonstrated  by  means  of  an  exact 
Gaussian  sohton  solution  in  the  case  where  the  saturable  nonlinearity  is  the  loga¬ 
rithmic  type.  The  conditions  necessary  to  establish  these  incoherent  soliton  states 
as  well  as  their  associated  characteristics  are  discussed  in  detail.  Our  analysis  in¬ 
dicates,  that,  for  a  given  strength  of  the  logarithmic  nonlinearity,  these  solitons 
exist  only  within  a  certain  range  of  incoherent  angular  spectral  widths.  Relevant 
computer  simulations  are  provided  to  further  elucidate  their  behavior. 
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4.2  Theoretical  formulation 


Let  us  consider  a  saturable  nonlinear  medium,  the  refractive  index  of  which  varies 
logarithmically  with  intensity  /,  i.e.  , 

n^I)  =  nl  +  n2ln{I/It)  (4.1) 

where  no  is  the  linear  refractive  index  of  this  material,  n2  is  a  positive  dimensionless 
coefficient  associated  with  the  strength  of  the  nonlinearity  and  It  is  a  threshold  in¬ 
tensity.  Note  that  a  logarithmic  nonlinearity  of  this  sort  has  been  recently  employed 
by  Snyder  and  Mitchell  in  their  study  of  mighty  morphing  spatial  solitons  and 
bullets  [80].  To  avoid  any  unphysical  singularities,  we  may  assume  that  the  In  (// It) 
nonlinearity  results  from  the  more  realistic  ln(l  + 1 /It)  model  in  the  limit  I  ^  It- 
Let  us  also  assume  that  the  light  beam  propagates  along  the  2;  axis  and  diffracts 
only  in  the  x  direction,  i.e.,  for  simplicity  we  are  dealing  with  a  planar  incoherent 
beam.  Moreover,  we  will  make  the  important  assumption  that  the  nonlinear  ma¬ 
terial  involved  responds  much  slower  than  the  characteristic  phase  actuation  time 
across  the  optical  beam  as  in  the  case  of  photorefractives.  As  a  result,  the  mate¬ 
rial  will  experience  only  the  time  averaged  intensity.  Under  these  conditions  and 
by  following  the  procedure  given  in  the  last  chapter,  one  may  then  show  that  the 
normalized  intensity  I^  =  I/It  of  this  incoherent  light  beam  obeys  the  following  set 
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of  equations: 


+  0 


df\  1  d'^f  kon2 .  (T  t  \\  f  n 


where 

c» 

In{x,z)=  J  \f  {x,z,e)f  dO 

— oo 


and  at  z  =  0, 


f{z  =  o,x,e)  =  r'i^a]i^(e)<l,„(x) 


(4.2) 


(4.3) 


(4.4) 


In  the  above  equations,  /  represents  the  so-called  coherent  density,  9  is  an  angle 
in  radians  with  respect  to  the  z  axis,  k  =  kouo  and  ko  =  27r/Ao  .  Gn  (9)  is  the 
normalized  angular  power  spectrum  of  the  incoherent  source,  (I>q  (x)  is  the  input 
spatial  modulation  function  and  r  is  an  intensity  ratio,  that  is  r  =  max  (/at)  • 
Note  that,  in  the  last  chapter,  the  process  was  described  in  terms  of  the  so-called 
coherent  components  which  are  in  fact  a  discrete  version  of  the  coherent  density 
employed  here.  Let  us  now  assume  that  the  input  intensity  pro  le  as  well  as  the 
incoherent  angular  power  spectrum  are  Gaussian,  i.e.,  (^o  (^)  ~  exp  (— a:^/2xo)  and 
Gat  (9)  =  exp  (^-9^/91^  In  this  case,  the  input  intensity  pro  le  can  be 

obtained  from  Eqs.4.3  and  4.4  and  it  is  given  by 


/at  (x,  z  =  0)  =  r  exp  {-x^/xfj  .  (4.5) 

For  convenience,  we  will  also  employ  the  following  dimensionless  coordinates  and 
variables;  ^  =  z/kx'^,  s  =  x/xq,  a  =  kxQ9  and  /?  =  (^2/2)  (koXo)^ . 
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At  this  point  let  us  also  assume  for  a  moment  that  the  intensity  of  the  inco¬ 
herent  beam  remains  invariant  or  stationary  during  propagation.  If  we  let  /  = 


F  exp{i/3^\nr),  we  then  obtain  from  Eqs.4.2  and  4.5 


dF  OF 

■^+“■57 


d^F 

2  as2 


^s'^F  =  0 


(4.6) 


where  at  ^  =  0,  F  =  gxp  (-6>V2^o)  (“•572)  and  f  =  1  Equation 

4.6  can  now  be  solved  using  a  Gaussian  wavepacket  solution[81],  i.e. 


F  =  A{^)  exp 


V 


2^2  (0 


exp  [5  (0  +  VP  iO  +  V^Q  iO]  ]  (4-7) 


where  77  =  s—v  (^).  w  (^)  represents  the  variable  spot-size  of  the  coherent  density  and 
V  (^)  is  the  displacement  of  this  wavepacket  along  the  x  axis.  Equation  4.7  satis  es 
Eq.  4.6  provided  that  the  following  relations  hold  true:  A^w  =  AqWq,  v=  oi+p,  q  —w 
/2w,  P=  -2l3v,  9=  (1/2)  (p2  -  w~^  -  20v^)  and  q=  -  (1/2)  {Aq‘^  -h  2/?  -  w"^)  where 
v=  dv/d^  etc.  The  initial  conditions  (at  C  =  0)  associated  with  these  equations  are 
given  by:  Aq  =  exp  ^^/20o)>  wo  =  !>  Po  =  =  0  and  w=  0.  In  that  case,  it 

can  be  readily  shown  that 


V  (0  =  a  (20)  sin  [(2^)^^^ 


(4.8) 


=  +  .  (4.9) 

The  rest  of  the  variables  A,  p,  g,  and  q  can  then  be  deduced  by  directly  employing 
the  results  of  Eqs.  4.8  and  4.9. 
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In  turn,  from  Eqs.  4.7,  4.8  and  4.9,  the  normalized  intensity  In  can  be  obtained 


and  it  is  given  by 


Q 


(4.10) 


where 

■  (4'“) 

Seeking  a  soliton  solution,  we  must  require  that  the  intensity  In  remains  station¬ 
ary  during  propagation.  Prom  Eqs.  4.5  and  4.10  it  is  evident  that  this  is  possible 
only  if  (5  =  1  for  all  values  of  z  or  Combining  the  results  of  Eqs.  4.9  and  4.11 
we  then  readily  obtain  the  condition  necessary  to  establish  a  partially  incoherent 
spatial  soliton,  which  is  2,0  —  1  =  or 


Xo 


1 

kosjna-  nl6l 


(4.12) 


Therefore,  given  the  strength  of  the  logarithmic  saturable  nonlinearity  n2  and  the 
width  of  the  incoherent  angular  power  spectrum  of  the  source  0O)  Eq.  4.12  uniquely 
determines  the  spatial  width  of  the  partially  incoherent  spatial  sohton  state  that 
can  propagate  in  this  material.  In  this  case,  the  partially  incoherent  soliton  beam 
maintains  an  invariant  Gaussian  intensity  pro  le,  i.e.  In  =  rexp  (— x^/xq),  even 
though  the  associated  coherent  density  function  /  evolves  considerably  as  a  function 
of  ^  or  z.  In  the  special  case  where  0o  0  or  Gn  {0)  =  ^{0),  our  result  correctly 
reduces  to  the  coherent  spatial  soliton  limit  (xq  =  previously  obtained  by 
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Snyder  and  Mitchell  [80].  On  the  other  hand,  Eq.  4.12  shows  that  as  9q  increases, 
the  soliton  width  xq  increases.  Nevertheless,  it  is  clear  from  Eq.  4.12,  that  there  is 
an  upper  limit  for  the  incoherent  angular  spectrum  width  6q  above  which  a  soliton 
beam  is  not  allowed.  This  upper  value  ^omax  is  equal  to  n2  /tiq. 

4.3  Results  and  discussion 

We  will  now  illustrate  our  results  by  means  of  relevant  examples.  Let  us  consider  a 
logarithmically  saturable  nonlinear  medium  with  no  =  2  and  n2  =  10“‘*.  Moreover 
let  the  optical  wavelength  Ao  be  0.5  jum.  For  this  set  of  parameters,  the  maximum 
source  angular  spectrum  width  do  that  can  be  self-trapped  is  ^omax  =  5  mrad.  Let 
us  also  assume  that  the  angular  power  spectrum  width  of  the  incoherent  source 
is  6o  =  4.5  mrad.  In  this  case  xq  ^  18  /xm  and  thus  the  intensity  FWHM  of  this 
partially-incoherent  spatial  soliton  is  ~  30  /xm.  Figure  4.1  shows  the  stationary 
propagation  of  this  soliton  as  obtained  by  numerically  solving  Eqs.  4.2  and  4.3.  In 
all  our  computations  we  have  used  the  better  behaved  In  (1  -|-  //7f)  nonlinearity  with 
an  intensity  ratio  of  one  thousand  (r  =  10^).  Note  that  had  this  same  incoherent 
beam  been  launched  in  a  linear  material  (no  =  2,  n2  =  0)  it  would  have  diffracted 
to  a  FWHM  of  ~  90  /xm  after  a  distance  of  1  cm.  Fig.4.2  on  the  other  hand,  depicts 
the  coherent  density  associated  with  this  soliton  beam  when  0  =  3  mrad. 

Evidently,  the  density  periodically  expands  and  contracts  during  propagation 
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Figure  4.1:  Stationary  propagation  of  a  partially-incoherent  spatial  soliton  when 
xo=18  Hm,  6 0=4. 5  mrad  and  n2=10-‘'. 
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Figure  4.3:  Propagation  of  a  partially  incoherent  Gaussian  beam  in  logarithmically 
saturable  nonlinear  medium  when  xo=18  tlm,  do=6.0  mrad. 
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Figure  4.4:  Propagation  of  a  partially  incoherent  Gaussian  beam  in  logarithmically 
saturable  nonlinear  medium  when  X6=18  jJm,  do=4.0mmd. 
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and  oscillates  around  the  center  of  the  beam.  In  spite  of  this,  the  overall  intensity 
pro  le  In  of  the  incoherent  light  beam  remains  stationary  as  long  as  Eq.  4.12  is 
satis  ed.  Figure  4.3  shows  what  would  happen  to  this  Gaussian  incoherent  beam  if 
6o  =  Q  mrad  {9o  is  above  the  upper  bound  0omax  =  5  mrad)  provided  that  again 
Xo  =  18  nm.  In  this  case,  the  beam  rst  expands  and  then  compresses  and  at 
z  =  l  cm  its  intensity  FWHM  is  ~  52  /xm.  If  is,  on  the  other  hand,  lower  than 
the  that  prescribed  by  Eq.  4.12,  say  6q  =  A  mrad  and  for  xq  =  18  fim,  then  the 
beam  exhibits  cyclic  compressions  and  expansions  as  shown  in  Fig.4.4.  Even  though 
our  analysis  does  not  directly  apply  to  biased  photorefractive  crystals  (where  the 
saturable  nonlinearity  is  different),  it  provides  valuable  insight  into  the  considerable 
more  complex  incoherent  dynamics  in  this  material  system.  In  fact,  the  scenario 
depicted  in  Figs.4.1,  4.3  and  4.4  is  quite  reminiscent  of  the  incoherent  wave  dynamics 
previously  encountered  in  biased  photorefractives  in  the  previous  chapter  [35], [79]. 

4.4  Conclusions 

In  conclusion  we  have  shown  that  partially  incoherent  spatial  solitons  are  possible 
in  saturable  nonlinear  media.  An  exact  Gaussian  solution  was  obtained  in  the  case 
where  the  nonlinearity  is  of  the  logarithmic  type.  These  incoherent  Gaussian  solitons 
can  exist  as  long  as  their  spatial  width  is  appropriately  interrelated  with  the  strength 
of  the  nonlinearity  and  the  width  of  the  incoherent  angular  power  spectrum  of  the 
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source.  The  behavior  of  these  beams  above  and  below  this  limit  was  also  investigated 
using  computer  simulations. 
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Chapter  5 

Multimode  Incoherent  Spatial 
Solitons  in  Non-Instantaneous 
Kerr  Media 

In  this  chapter,  by  using  the  self-consistent  multimode  approach,  we  show  that  in¬ 
coherent  sohtons  are  possible  in  non-instantaneous  Kerr-like  nonlinear  media.  The 
coherence  characteristics  of  multimode  incoherent  spatial  solitons  in  this  material 
system  are  investigated.  Other  properties  of  these  incoherent  solitons  are  also  dis¬ 
cussed  as  a  function  of  their  modal  composition. 
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5.1  Introduction 


In  Chapter  3,  the  coherent  density  approach,  describing  propagation  and  self-focusing 
of  partially  spatially-incoherent  light  beams  in  nonlinear  media  is  developed.  In  the 
last  chapter,  this  approach  was  found  to  lead  to  closed  form  solutions  for  logarith¬ 
mic  type  of  nonlinearities  [82].  It  has  been  also  shown  [83],  that  a  self-consistent 
multimode  method  can  also  be  used  to  investigate  incoherent  soliton  formation  in 
nonlinear  materials.  This  latter  model  seeks,  in  a  self-consistent  fashion,  multimode 
soliton  solutions  whose  total  intensity  can  be  obtained  via  intensity  superposition 
of  all  the  modes  guided  within  the  nonlinearly  induced  waveguide.  This  method 
is  capable  of  identifying  incoherent  spatial  soliton  states,  their  range  of  existence 
and  coherence  properties  [83].  Another  approach  for  describing  broad  incoherent 
bright  solitons  was  also  suggested  by  Snyder  and  Mitchell  [84].  This  ray  model  is 
to  some  extend  related  to  the  Vlasov  transport  description  previously  suggested  by 
Hasegawa  in  the  theory  of  random-phase  solitons  in  plasmas  [85],  [86],  [87]. 

It  is  worth  noting  that  incoherent  solitons,  as  viewed  from  the  perspective  of  the 
self-consistent  modal  theory,  are  in  fact  related  to  the  so-called  incoherently  cou¬ 
pled  solitons  in  photorefractives  [32]-[34],[50]  or  to  the  vector  solitons  in  Kerr  media 

[88] .  Within  the  context  of  incoherently  coupled  solitons  [32] -[34],  [50],  Vysloukh  et  al 

[89] , [90]  has  shown  that  multi-component  coupled  soliton  modes  can  be  incoherently 
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superimposed  in  weakly  saturating  photorefractive  crystals.  We  would  like  to  em¬ 
phasize  however  that  there  is  a  subtle  difference  between  self-trapping  an  incoherent 
beam  and  creating  a  multi-component  soUton.  This  difference  is  due  to  the  statisti¬ 
cally  varying  modal  weights  of  the  incoherent  beam.  Consider  rst  a  situation  where 
self-trapping  occurs  from  an  incoherent  light  source.  During  self-trapping,  the  beam 
continuously  excites  several  modes  of  the  jointly  induced  waveguide,  established  by 
the  time  averaged  intensity  of  the  speckled  beam  itself.  If  the  time  averaged  inten¬ 
sity  is  decomposed  into  the  modes  of  the  induced  waveguide,  each  mode  will  have  a 
certain  coupling  amount  or  weight.  On  the  other  hand,  when  multiple  laser  beams 
are  incoherently  superimposed  (engineering  the  beam  pro  les),  the  relative  modal 
weights  are  time  independent  and  only  the  relative  phase  between  each  pair  of  modes 
changes  with  time.  As  far  as  the  crystal  could  tell,  the  two  situations  are  exactly 
the  same.  The  difference  lies  in  the  time  dependence  of  the  mode  occupancies. 

In  this  chapter  we  demonstrate  that  incoherent  spatial  solitons  are  possible  in 
non-instantaneous  Kerr-like  media.  Closed  form  solutions  are  obtained  using  the 
self-consistent  multimode  approach  [83] .  It  is  shown  that  the  intensity  pro  le  of  these 
incoherent  soliton  states  is  of  the  sech^  {x/xq)  type.  Moreover,  our  analysis  indicates 
that,  in  this  case,  the  peak  intensity  and  spatial  width  of  these  incoherent  spatial 
solitons  are  related  to  the  number  of  allowed  modes.  The  coherence  properties  of 
these  soliton  states  are  investigated  in  detail.  Relevant  examples  are  provided. 
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5.2  Self-consistent  multimode  approach 


We  consider  a  spatially  partially-incoherent  optical  beam  which  propagates  in  a 
nonlinear  self-focusing  Kerr-like  medium  along  the  z-axis.  For  simplicity  we  assume 
that  this  beam  is  planar,  that  is,  it  diffracts  only  along  x  the  direction.  The  re¬ 
fractive  index  of  this  Kerr-like  material  varies  linearly  with  optical  intensity  /  i.e. 
let  ri^  =  nl  +  712/,  where  no  is  the  linear  part  of  the  refractive  index  and  n2  is  the 
nonlinear  Kerr  coefficient.  We  also  make  the  important  assumption  that  the  nonlin¬ 
earity  responds  much  slower  that  the  characteristic  phase  uctuation  time  across  the 
incoherent  beam,  so  as  to  avoid  speckle-induced  lamentation  instabilities  [91].  In 
this  regime,  the  material  will  experience  only  the  time  averaged  beam  intensity.  For 
example,  such  noninstantaneous  Kerr  nonlinearities  can  be  encountered  in  biased 
photorefractives  in  the  low-intensity  regime,  (i.e.,  when  the  so-called  dark  irradiance 
is  much  larger  than  the  intensity  of  the  optical  beam  [23], [24])  or  in  materials  with 
thermal  nonlinearities  [92], [93].  Furthermore,  let  the  total  electric  eld  E  of  this 
spatially  incoherent  beam  be  written  in  terms  of  a  slowly  varying  envelope  U,  that 
is,  let  E  =  U  exp{ikz),  where  the  wavevector  k  is  given  hy  k  —  koUQ  —  (27r/Ao)no 
and  Ao  is  the  free-space  wavelength.  In  this  case,  it  can  be  readily  shown  that  the 
envelope  U  evolves  according  to: 


^  dz  ^  2k  dx'^  ^ 


kln2l 

2k 


U  =  0 


(5.1) 
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Let  us  now  assume  that  the  incoherent  spatial  soliton  in  this  Kerr-like  medium  has 
a  sech^  {x/xq)  intensity  pro  le,  that  is 

I  =  Jo  sec  h^  {x/xq)  (5.2) 


where  represents  the  peak  intensity  of  the  optical  incoherent  multimode  beam  and 
xq  is  associated  with  its  spatial  extent.  In  this  case,  Eq.  5.1  takes  the  form 


.dU  Id^U  u2/xrr  n 


(5.3) 


In  this  equation  we  have  used  normalized  coordinates  and  quantities,  i.e.,  rj  =  x/xq, 
i  =  z/kxl,  and  =  klxln2h.  The  incoherent  spatial  sohton  solutions  of  this 
equation  can  then  be  obtained  by  expressing  the  optical  envelope  through  a  super¬ 
position  of  all  the  modes  involved,  that  is,  17  a  (??)  exp  (^),  where  are 

m 

the  mode-occupancy  coefficients  that  vary  randomly  with  time,  w™  {rf)  is  the  pro  le 
of  the  mth-order  mode,  and  (3^  its  phase  constant.  The  discrete  index  n  stands  for 
the  total  number  of  modes  allowed  in  the  system.  By  substituting  this  form  of  U 
into  Eq.  5.3,  we  obtain,  for  each  mode,  the  following  ordinary  differential  equation: 


+  [a^  sec  h^  (r/)  -  2^  J  <  =  0  .  (5.4) 

The  mode  pro  les  u^(r])  and  their  phase  constants  /3„  can  then  be  found  by  solving 
Eq.  5.4.  In  order  to  nd  the  mode  occupancy  coefficients,  the  self-consistency  of  this 
solution  must  be  satis  ed,  i.e.,  the  total  intensity  of  this  incoherent  beam  is  given 
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by  Eq.  5.2.  The  self-consistency  method  relies  on  the  basic  concept  that  a  soliton 
itself  is  a  mode  of  the  self-induced  (via  the  nonlinearity)  waveguide.  This  was  rst 
suggested  by  Askar  yan  [94]  and  later  developed  into  a  self-consistency  methodology 
by  Snyder  at  al  [95] -[98]  as  means  to  analyze  coherent  solitons. 


5.3  Results  and  Discussion 


To  nd  the  modal  structure  of  these  multimode  incoherent  spatial  solitons,  we  rst 
adopt  the  transformation  t  =  tanh  (r?),  in  which  case  Eq.5.4  takes  the  form: 


<1  _  fA 
<  )  dt^ 


2t 


d< 

dt 


+ 


2/?„ 


(l-t2)J 


<  =  0 


(5.5) 


It  can  now  be  easily  shown  that,  when  =  n{n  +  1),  where  n  =  1, 2, (an  integer 
determined  by  the  amount  of  nonlinearity)  and  m  =  1,  2,  the  allowed  modes 
can  be  expressed  in  terms  of  associated  Legendre  functions  [99],  i.  e., 


<  (77)  =  Pr  (taiih(77)) 


(5.6) 


and  Pjn  =  where  P™  (x)  are  associated  Legendre  functions  of  the  rst  kind. 

In  general,  {x)  =  (1  -  d^P„  (x)  /dx^  where  P„  (x)  denote  Legendre  poly¬ 

nomials.  The  value  of  the  integer  n  (associated  with  or  the  degree  of  non¬ 
linearity)  represents  the  number  of  allowed  modes  in  the  self-induced  waveguide 
and  m  the  mode  index  for  this  case.  Note  that  the  lowest  order  mode  occurs 
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when  m  =  n  whereas  the  highest  at  m  =  1.  For  a  given  value  of  n,  the  statis¬ 
tically  varying  optical  eld  of  the  incoherent  spatial  soliton  is  given  by  t/  (r?,  0 
E  (tanh(77))exp(im2^/2).  Note  that  the  self-consistency  condition  also  re- 

m=l 

quires  that  the  mode-occupancy  coefficients  also  depend  on  the  number  of  modes  n. 
The  time  averaged  intensity  of  this  beam  can  then  be  evaluated  from 

(|c/p>  oc  /  =  i:  (ICI')  IJT  (tanh(,))l^  ,  (6.7) 

m=l 

where  we  made  use  of  the  fact  that  under  incoherent  excitation  the  time-average 

of  cross-interference  terms  among  modes  is  zero,  oc  Sij.  Finally,  Eq.  5.7 

can  be  used  to  obtain  the  mode-occupancy  coefficients  that  self-consistently  lead 

to  the  incoherent  spatial  soliton  I  =  lo  sech^  [r])  assumed  in  the  very  beginning  of 

this  analysis.  It  proves  convenient  to  normalize  the  eigenmodes  in  the  following 

fashion:  {v)  =  or  equivalently  Iq  E  =  hsech^  {rf). 

L  “*0  J  m—l 

The  functional  form  of  the  normalized  eigenfunctions  {rj)  is  given  in  Appendix  A. 
Figure  5. 1  depicts  the  normalized  mode  pro  les,  u™  (77)  for  a  three-component  (n  = 
3)  multimode  incoherent  soliton  whereas  Fig.  5.2  for  a  four  component  structure 
(n  =  4). 

Having  found  closed  form  solutions  for  these  spatial  incoherent  solitons,  the 
question  naturally  arises  as  to  what  factors  contribute  to  their  spatial  FWHM  (full- 
width  half-maximum).  To  answer  this  question  one  has  to  consider  the  expression 
=  n{n  +  l).  More  speci  cally,  since  =  A;gx§n2/o,  then  given  the  number  of 


guided  modes  n,  and  the  peak  intensity  Jo,  the  spatial  width  xq  can  be  determined. 
In  turn  one  quickly  nds  the  following  result:  Xq  =  [n  (n  +  1)  /A:on2/o]  '  •  Using  Eq. 
5.2  and  provided  that  n  is  known,  the  intensity  FWHM  of  these  soliton  states  u;„, 


can  be  calculated  from  Wn  =  1.76xo,  and  thus 

1.76  [n  (n  +  1) 


Note  that  when  n  =  1  (  single-mode  soliton),  this  latter  expression  is  in  full 
agreement  with  the  results  previously  obtained  for  coherent  Kerr  spatial  solitons 
[100].  Equation  5.8  can  also  be  written  in  terms  of  a  normalized  FWHM  (wn/^o) 
as  function  of  the  maximum  induced  nonlinear  index  change  (An  =  n2/o))  i-e., 
Wn/Ao  =  (1.76/27r)  [n  (n  -f  1)  /An]^^^.  Figure  5.3  shows  the  dependence  of  the  nor¬ 
malized  intensity  FWHM  of  these  incoherent  optical  solitons  on  An  for  different 
values  of  n.  This  gure  clearly  shows  that  for  a  given  nonlinear  index  change,  the 
intensity  FWHM  tends  to  increase  with  the  number  of  modes  involved.  This  should 
have  been  anticipated  since  the  incoherency  increases  with  n.  It  is  also  important 
to  note  that  in  Fig.5.3,  the  existence  curves  are  discrete  as  opposed  to  continuum 
range  found  in  Ref.  [83].  This  is  due  to  the  speci  c  sech^  {x/xq)  intensity  shape 


assumed  at  the  beginning  of  our  analysis. 

We  now  investigate  the  coherence  properties  of  these  spatial  solitons.  This  is  done 
by  evaluating  the  complex  coherence  factor  fj,i2  of  these  soliton  states  [69] ,  which  is 
given  by /Xi2  {'n,’n  +  <5)  =  {U  (77,0  U*  (t?  ■+  5,0)  /  [(l*^  (OP)  {\U  (v  +  5)1^)]  ,  where 
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Figure  5. 1 :  Normalized  mode-profiles  as  a  function  of  Tj  for  n=3. 
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Figure  5.2:  Normalized  mode-profiles  as  a  function  of  t]  for  n=4. 
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6  represents  the  normalized  distance  between  two  points.  Prom  Eq.  5.2,  and  by  using 
the  modal  expansion  of  the  optical  envelope  U,  then  for  a  given  n,  we  obtain 

E  iv)  Un  iV  +  S) 

fj'12  iVi  ^7  +  ^)  =  '""V  /  s - .  /  c<-  >  (5-9) 

^  sec  h  (ry)  sec  h  (77  +  o) 

where  again  we  made  use  of  the  fact  that  the  time-average  of  cross-interference  terms 
is  zero  under  incoherent  excitation.  The  complex  coherence  factors  /ii2  (r?,  rj  +  6)  up 
to  n  =  5  have  been  tabulated  in  Appendix  B.  When  n  =  1,  which  corresponds  to 
a  coherent  optical  soliton,  ^^2  (?7)  ^  =  1  ^  should  have  been  anticipated.  This 

in  turn  implies  an  in  nite  correlation  length.  As  the  number  of  modes  n  increases, 
the  coherence  properties  of  these  incoherent  soliton  states  become  more  complicated. 
Figure  5.4  provides  a  two-dimensional  plot  of  fj.12  as  a  function  of  normalized  position 
T]  and  point  separation  6  when  n  —  2.  As  expected,  /Xj2  1  when  6  =  0  and  this  true 
for  every  value  of  n.  Note  that  any  cut  {ij  =constant)  of  these  two  dimensional  plots 
provides  the  correlation  (magnitude  and  phase)  between  the  two  points  located 
at  77  and  77  +  6.  Moreover,  as  it  is  well  known  [69],  the  magnitude  of  the  complex 
coherence  factor  is  always  I/X12I  ^  !•  In  general,  two  points  are  mutually  coherent  if 
/ij2  —  1  nnd  mutually  incoherent  (poor  correlation)  if  fii2  ~  0.  Figure  5.5  depicts  a 
cross-section  of  the  complex  coherence  ^12  as  obtained  right  at  the  center  (77  =  0)  of 
an  incoherent  soliton  when  n  —  2.  In  this  case,  n^2  (0;  ~  sech(6),  which  implies 

that  correlation  between  the  center  and  any  other  point  of  this  soliton  decreases  with 
6.  This  latter  expression  also  demonstrates  that  the  tails  of  this  beam  (6  — >  ±00) 
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and  the  center  are  totally  uncorrelated.  Fig.5.6  on  the  other  hand  shows  how  a  point 
at  the  tails  (at  77  =  -3)  correlates  with  the  rest  of  the  beam  when  again  n  =  2.  For 
example,  when  5  ft;  3,  in  which  case  the  second  point  is  close  to  the  beam  center, 
A*i2  ~  0  (poor  correlation),  in  agreement  with  Fig.5.5.  As  6  — >  —00,  (i.e.  both  points 
are  on  the  left  tail  of  the  beam)  there  is  maximum  correlation  [11^2  —  !)•  This  high 
degree  of  coherence  is  due  to  the  following  fact:  at  the  tails  there  is  essentially 
only  one  mode,  (the  highest-order  mode,  m  =  1)  which  is  coherent  in  itself.  For 
^  +00,  /ii2  ~  -1  and  this  tt  phase-shift  is  due  to  the  anti-symmetric  character  of 

this  highest-order  mode.  Fig.5.7,  provides  the  same  information  at  the  intermediate 
point  rj  =  -0.4,  when  n  =  2.  Note  that  all  these  cross-sections  are  very  different  in 
character  from  those  found  in  the  case  of  logarithmically  saturable  nonlinear  media 
(see  Chapter  6).  In  Log  systems,  the  statistical  process  is  everywhere  stationary, 
i.e.  /Lii2  depends  only  on  6  and  it  is  Gaussian  in  nature.  As  shown  here  for  the 
Kerr  case,  the  coherence  curves  depend  also  on  the  position  77  and  do  not  always 
go  to  zero  as  6  ^  ±00.  This  feature  was  also  encountered  in  saturable  Kerr  media 
of  the  type  //(!  +  /)  (as  will  be  seen  in  Chapter  7).  This  important  difference  is 
due  to  the  nite  number  of  modes  associated  with  the  Kerr  incoherent  solitons.  In 
particular,  the  induced  waveguides  in  Kerr  and  Kerr-saturable  media  exhibit  cutoffs 
(  nite  number  of  modes)  whereas  the  logarithmic  one  does  not.  Figure  5.8  depicts 
the  coherence  function  for  tt  =  3  .  Similar  cuts  are  also  provided  in  Figs.  5.9, 
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5.10,  and  5.11.  Fig.5.9  shows  how  the  center  correlates  with  rest  of  the  beam  for 
n  =  3.  Interestingly  enough,  in  this  case,  there  is  a  weak  correlation  between  the 
center  (rj  =  0)  and  the  tails  (5  -»  ±oo).  This  is  because  the  highest-order  mode 
(m  =  1),  which  survives  at  the  tails,  happens  to  be  even  and  it  is  thus  nite  at 
the  center.  Fig.5.10  gives  the  correlation  between  a  point  at  the  tails  (77  =  —3)  and 
the  rest  of  the  beam.  Note  the  high  degree  of  coherence  (//12  1)  between  the  two 

tails.  Unlike  the  case  shown  in  Fig.5.6,  the  correlation  is  in  phase  since  again  in 
this  case  the  highest-order  mode  is  even.  Finally  Fig.5.11  gives  the  coherence  factor 
at  r;  =  sech“^  (2/a/5)  ~  ±0.4812  when  n  =  3.  Using  the  coherence  factor  given 
in  Appendix  B,  it  can  be  readily  shown  that  these  particular  points  are  completely 
uncorrelated  with  the  tails  of  this  incoherent  soliton.  Figures  5.12  and  5.13  depict 
the  coherence  factors  for  n  =  4  and  5,  respectively.  It  is  interesting  to  observe  the 
similarity  between  the  7/12  surfaces  when  n  is  odd  or  even.  Clearly,  as  the  order  of 
the  soliton  increases,  the  coherence  surfaces  become  more  involved. 

5.4  Conclusions 

In  conclusion,  we  have  shown  that  multimode  incoherent  spatial  solitons  are  pos¬ 
sible  in  non-instantaneous  Kerr-like  media.  Closed  form  solutions  were  obtained 
using  the  self-consistency  approach  provided  that  their  intensity  pro  le  is  of  the 
sech^  (x/xo)  type.  The  coherence  properties  of  these  incoherent  soliton  states  were 
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Figure  5.4:  Spatial  coherence  function  as  a  function  of  T]  and  8  when  n=2. 
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Figure  5.6:  Cross-section  of  the  coherence  function  when  n=2  and  r]=-3. 
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Figure  5.8:  Spatial  coherence  function  as  a  function  of  rj  and  5  when  n=3. 
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Figure  5.12:  Spatial  coherence  function  as  a  function  of  t]  and  S  when  n-4. 
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Figure  5.13:  Spatial  coherence  function  as  a  function  of  rj  and  5  when  «=5. 
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also  investigated  in  detail  and  explained  by  means  of  their  modal  composition. 
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Chapter  6 


Multimode  incoherent  solitons  in 
saturable  nonlinear  systems 


In  this  chapter,  we  show  that  multimode  incoherent  spatial  sohtons  are  possible  in 
log-type  saturable  nonhnear  systems.  The  mode-occupancy  function  associated  with 
these  soliton  states  is  found  to  obey  a  Poisson  distribution.  Our  analysis  indicates 
that  two  approaches,  i.e.  the  dynamic  coherent  density  description  as  well  as  the 
static  self-consistent  multimode  method  lead  to  exactly  the  same  results.  Closed 
form  solutions  are  obtained  for  (l-t-l)  D  as  well  as  for  (2-t-l)  D  circular  and  elhptical 


incoherent  solitons. 


6.1  Introduction 


Using  coherent  density  approach  it  has  been  shown  that  (1+1)  D  incoherent  spatial 
solitons  are  in  fact  possible  in  saturable  nonlinear  media  of  the  logarithmic  type. 
It  is  important  to  note  however,  that  in  general,  the  coherent  density  model  is 
by  nature  better  suited  to  describe  dynamic  evolution  and  thus  it  does  not  readily 
lend  itself  toward  identifying  stationary  solutions  such  as  incoherent  spatial  solitons. 
With  the  exception  of  the  saturable  logarithmic  model  (perhaps  the  only  exactly 
soluble  model  in  the  coherent  density  description),  the  question  naturally  arises  as  to 
whether  incoherent  spatial  solitons  also  exist  (or  not)  in  other  saturable  nonhnear 
systems  (such  as  in  biased  photorefractives).  In  order  to  address  this  problem, 
we  developed  an  alternative  theory  which  is  essentially  a  self-consistent  incoherent 
multimode  approach  [83].  In  brief,  in  this  procedure,  an  incoherent  spatial  soliton 
is  sought,  which  intensity- wise  is  a  superposition  of  all  the  modes  self-consistently 
guided  in  its  nonlinearly  induced  waveguide.  By  doing  so,  (1+1)  D  photorefractive 
incoherent  spatial  solitons,  the  conditions  necessary  for  their  existence  as  well  as 
their  coherence  characteristics  were  obtained  [83].  Furthermore,  in  the  last  chapter, 
by  using  this  approach,  we  showed  that  multimode  incoherent  spatial  solitons  are 
possible  in  non-instantaneous  Kerr-like  media.  At  a  rst  sight,  the  two  theories 
[79], [83]  may  seem  to  have  very  little  in  common.  Thus  it  is  important  to  ask 
whether  the  two  approaches  are  mutually-consistent  or  lead  to  the  same  results. 
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Moreover,  it  is  interesting  to  know  whether  (2+1)  D  incoherent  spatial  solitons  are 
also  possible  in  saturable  nonlinear  media. 

In  this  chapter  we  show  that  multimode  incoherent  spatial  solitons  are  possi¬ 
ble  in  logarithmically  saturable  nonlinear  media.  The  mode-occupancy  function  of 
these  soliton  states  is  found  to  obey  a  Poisson  distribution.  Our  analysis  demon¬ 
strates,  that  in  this  case,  the  two  theories,  i.e.  the  dynamic  coherent  density  de¬ 
scription  and  the  static  self-consistent  multimode  method  lead  to  identical  results. 
Two-dimensional  circular  and  elliptical  incoherent  spatial  solitons  along  with  their 
characteristics  are  also  obtained  in  closed  form.  Even  though  the  saturable  logarith¬ 
mic  nonlinearity  differs  from  the  photorefractive,  it  provides  nevertheless  a  platform 
(perhaps  the  only  platform)  upon  which  the  equivalence  of  the  two  previously  men¬ 
tioned  approaches  can  be  established  in  closed  form.  Finally,  as  recently  argued  by 
Snyder  and  Mitchell  [101]  and  by  Shen  [102],  simpli  cations  of  this  sort  in  terms 
of  accessible  nonlinearities,  can  provide  valuable  insight  and  still  maintain  the 
characteristic  features  of  the  underlying  physical  process. 


6.2  Incoherent  multimode  bright  solitons 

We  begin  our  analysis  by  considering  a  saturable  nonlinear  medium  of  the  loga¬ 
rithmic  type,  similar  to  that  previously  considered  by  Snyder  and  Mitchell  in  their 
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(6.1) 


study  of  mighty  morphing  spatial  solitons  and  bullets  [80],  i.e. 

n2  (/)  =n2  +  n2ln(7//t) 


where  no  is  the  Unear  refractive  index  of  the  material,  n2  is  a  positive  dimension¬ 
less  coefficient  associated  with  the  strength  of  the  nonlinearity  and  It  is  a  threshold 
intensity.  To  avoid  any  logarithmic  singularities  we  assume  that  the  ln(///t)  non¬ 
linearity  results  from  the  more  realistic  ln(l  -|- 1 /It)  model  in  the  limit  I  ^  It  ■ 
Let  the  optical  beam  propagate  in  the  z  direction.  We  also  make  the  important  as¬ 
sumption  that  the  nonlinearity  responds  much  slower  than  the  characteristic  phase 
actuation  time  across  the  incoherent  beam  so  as  to  avoid  speckle-induced  lamen¬ 
tation  instabilities  [91].  Thus,  in  this  regime  the  material  will  experience  only  the 
time  averaged  beam  intensity  as  in  the  case  of  photorefractives.  Furthermore  let 
the  electric  component  of  the  optical  eld  be  written  in  terms  of  a  slowly  varying 
envelope  U,  i.e.,  E  =  U e^{ikz)  where  k  =  fco«o  and  ko  =  27r/Ao.  In  that  case, 
it  can  be  readily  shown  that  the  envelope  U  in  this  nonlinearly  induced  waveguide 
evolves  according  to 


.dU  1  fd'^U  d^U\  a 


+ 


-I- 


+  —\n{I:,)U  =  0 


(6.2) 


2  y  ds^  dr/^ 

where  In  =  I /h  is  a  normalized  intensity  and  in  Eq.  6.2  we  have  used  normalized 
coordinates  and  quantities,  that  is,  s  =  x/wq,  r]  =  y/wo,  ^  =  z/kwl,  o?  =  n2{koWo)‘^ 
where  Wo  is  an  arbitrary  spatial  scale  or  spot  size. 
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6.3  One  dimensional  incoherent  solitons 


We  rst  employ  the  self-consistent  incoherent  multimode  description  in  the  case 
when  the  incoherent  beam  is  one-dimensional  or  planar.  For  the  time  being,  let  us 
assume  that  an  incoherent  spatial  soliton  of  a  Gaussian  intensity  pro  le  exists,  i.e. 
In  =  rexp(— s^)  where  r  is  an  intensity  ratio  with  respect  to  It.  From  these  last 
assumptions,  Eq.  6.2  takes  the  form 


^  2  ds^ 


+  ^[ln(r)-s^]Cf  =  0 


(6.3) 


The  allowed  modes  in  this  parabolic  waveguide  can  then  be  easily  obtained  in  terms 
of  Gauss-Hermite  functions  [103]  and  thus  the  optical  eld  can  be  expressed  through 

OO 

superposition,  i.e.,  U  =  CmUm  where 

m=0 

=  Hm  exp  (-q;sV2)  exp  {iPrnO  (6-4) 


and  Cm  are  the  mode-occupancy  coefficients  that  vary  randomly  in  time.  In  Eq. 
6.4,  Hm{x)  are  Hermite  polynomials,  /?„,  =  (1/2)  [Q;^ln(r)  —  (2m -I-  l)a]  and  m  = 
0, 1,2, ...  .  The  time-average  intensity  of  this  beam  can  then  be  obtained  from 

lN  =  {\Uf)='E{\Cmf)\Umf  (6-5) 

m=0 

where  we  made  use  of  the  fact  that  the  time  average  of  the  cross-interference  terms 
is  zero  ((cjcj^  oc  6ij)  under  incoherent  excitation  [83], [104].  Yet,  at  this  point,  we 
still  do  not  know  the  mode-occupancy  function  that  self-consistently  leads 
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to  the  incoherent  Gaussian  soliton  In  -  rexp{-s^)  assumed  in  the  very  beginning 
of  this  analysis.  Let  the  mode-occupancy,  as  in  the  case  of  quantum  mechanical 
coherent  or  Glauber  states  [69], [103],  be  described  by  a  Poisson  distribution  [105], 


Prom  Eq.6.7  it  is  evident  that  in  this  case,  the  positive  Poisson  parameter  p  must 

be  below  unity  (p  <  1)  for  this  incoherent  soliton  to  exist.  Or  alternatively,  from 

1  /2 

Eq.6.9  ,  a>  1  which  implies  that  at  a  given  wavelength  Ao,  the  quantity  na  wq  has 
a  cut  off,  below  which  no  incoherent  spatial  solitons  are  allowed. 

To  establish  the  equivalence  between  this  result,  Eqs.6.8  and  6.9,  with  that 
previously  obtained  using  the  coherent  density  approach  in  Chapter  4  one  must 
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resort  to  diffraction  data.  More  speci  cally,  during  linear  diffraction  (n2  =  0),  the 
intensity  of  each  Gauss-Hermite  mode  evolves  according  to  [103]  (see  Appendix  C) 


:  exp  - 


as 


I 

l  +  a^e) 


.VT+^. 


(6.10) 


By  keeping  in  mind  that  the  mode-occupancy  function  satis  es  the  Poisson 

distribution  of  Eq.6.6  and  by  employing  again  Eqs.6.7  and  6.9,  we  nd  that  upon 
diffraction  the  normalized  intensity  pro  le  Id  (see  Appendix  D)  is  given  by 


In  = 


(6.11) 


Equation  6.11  clearly  shows  that  during  diffraction  this  incoherent  beam  remains 
Gaussian.  Given  the  fact  that  the  initial  modulation  function  [79]  is  Gaussian,  i.e. 
,^p(s)  =  exp(-s^/2)  and  that  the  diffraction  behavior  of  this  statistically  stationary 
beam  obeys  a  convolution  integral,  Eq.3.3  of  Chapter  3,  it  then  becomes  immediately 
apparent  that  the  angular  power  spectrum  Gn{9)  of  the  incoherent  source  must 
also  be  Gaussian.  In  this  regard,  let  Gff{6)  =  exp(— where  9o  is 

the  width  of  the  incoherent  angular  power  spectrum.  Thus,  had  we  treated  the 
diffraction  regime  within  the  coherent  density  approach  [71]  (see  Appendix  E),  we 
would  have  arrived  at  the  following  result 

where  V  =  kwo9o.  Since  diffraction-wise,  the  results  of  Eqs.  6.11  and  6.12  are  the 
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same,  we  conclude  that  =  1  +  or  that 


Wq 


1 

ko\Jn2-niel 


(6.13) 


The  result  of  Eq.6.13  is  exactly  what  has  been  previously  obtained  in  Chapter  4, 
Eq.4.12,  using  the  coherent  density  approach.  At  a  more  fundamental  level,  our 
analysis  demonstrates  that  the  two  theories,  i.e.  the  dynamic  coherent  density 
description  as  well  as  the  self-consistent  multimode  method  lead  in  this  case  to 
exactly  the  same  results.  In  the  coherent  Umit  0  or  a  ^  1,  the  soliton  spot- 
size  reduces  to  wo  =  as  previously  found  by  Snyder  and  Mitchell  [80].  In 

this  same  limit,  Eq.  6.9  suggests  that  the  Poisson  parameter  p  ^  0  and  thus  the 
beam  itself  is  single-moded  or  fully  coherent.  Moreover,  the  width  6q  of  the  angular 
power  spectrum  of  the  incoherent  source  and  the  Poisson  parameter  p  are  now 
related  through  the  following  expression  p  =  [(1-1-  +  l]. 

Of  course  these  incoherent  spatial  soliton  states  exist  as  long  as  p  <  1  or  a  >  1  or 
equivalently  Oq  <  jn^. 


6.4  Two  dimensional  incoherent  solitons 

In  the  same  vein,  one  may  show  that  two-dimensional  incoherent  spatial  solitons 
also  exist  in  this  nonlinear  system.  To  do  so,  let  us  assume  in  general  an  elliptic 
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Gaussian  incoherent  beam, 


Iff  =  r  exp  — (s^  +  <T^77^)j 


(6.14) 


where  the  parameter  a  is  associated  with  its  degree  of  ellipticity.  Substituting  Eq. 
6.14  in  Eq.  6.2  we  get 


dU  1  fd'^U  d‘^U\  aV,  ,  ,  2  2  2l  rr  n  fa 

5  +  a?  j  +  T  C)  -  *  -  "  ”  1  ^  •  <«■ 


15) 


The  optical  eld  U  in  this  elliptic  paraboloid  potential  can  be  obtained  as  a  super- 

OO  OO 

position  of  the  allowed  Gauss-Hermite  modes,  i.e.  U  =  Cm,nUm,n  where 

m=0  n=0 


Um,n  =  Hm  Hn  exp  [-  {a/2)  (s^  -t  W)]  exp  {i/3rn,n^)  (6.16) 

where  „  =  {a/2)  [Q!ln(r)  -  (2m  -I- 1)  -  a{2n  +  1)].  The  overall  intensity  of  this 
beam  can  then  be  found  from 


CX>  OO 

(6.17) 

m=0  n=0 

where  we  have  assumed  again  that  under  incoherent  excitation  {cijch)  a  SikSji.  If 
we  now  allow  the  mode-occupancy  ^|cm,n|^)  to  obey  a  double  Poisson  distribution 


(p/2)-(g/2)" 

m\n\ 


(6.18) 


then  the  two  dimensional  incoherent  spatial  sohton  Iff  =  rexp  [— (s^  -t-  a^ff))  can 
be  self-consistently  recovered  using  Eqs.  6.7,  6.17,  and  6.18  provided  that  the  two 
Poisson  parameters  p,  q  satisfy 


P  = 


a  —  1 
a  +  1 


(6.19) 
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a  —  a 


(6.20) 


^  a  +  a 

and  moreover  f  =  r  [(1  —  p^)  (1  —  exp  [(p  +  g)  /2] .  These  incoherent  solitons 
can  be  elliptical  (a  7^  1)  or  circular  (a  =  1)  depending  on  whether  the  angular 
power  spectrum  of  the  incoherent  source  (exciting  this  structure)  is  symmetric  or 
not.  Moreover,  these  states  exist  as  long  as  0  <  (p,  9)  <  1  or  a  >  max(l,cr). 
The  possibility  of  generating  (2+1)  D  elliptical  solitons  in  isotropic  nonlinear  media 
seems  to  be  unique  to  incoherent  multi-mode  solitons,  since  it  has  been  shown  that 
their  coherent  counterparts  change  their  widths  periodically  during  propagation  [80]. 

6.5  Coherence  properties 

Having  found  the  modal  composition  of  these  solitons,  their  coherence  properties 
[69]  can  then  be  described  by  following  procedures  similar  to  those  employed  in 
Refs.  [104]  and  [107].  Using  Eqs.6.7,  and  6.16-6.20,  one  can  show  that  the  complex 
coherence  factor  Pi2  of  these  solitons  is  given  by 

Pi2(s,?7;s  +  <5,p  +  e)  =  exp|-  ^ 

This  latter  result  demonstrates  that  the  coherence  function  pi2  is  independent  of 
position  (s,  p)  within  the  soliton  beam  and  instead  it  depends  only  on  the  deviation 
distances  S  and  e.  For  a  circular  incoherent  soliton  (p  =  9,  cr  =  1),  the  actual 
correlation  length  Ic  (the  distance  where  Pi2  to  its  e"^  value)  can  be  obtained 
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from  Eq.6.21  and  it  is  given  by  Ic  =  iyo(l  “  Evidently,  Ic—^oo  when  p  =  0 

(single  mode  case)  whereas  h  0  when  p  ->  1.  These  results  are  in  agreement 
with  our  previous  discussion.  In  fact,  for  one  dimensional  solitons,  the  complex 
coherence  factor  Pi2  is  identical  to  the  statistical  autocorrelation  function  of  the 
source  R{xi  -  X2),  R{xi  -  X2)  can  be  obtained  from  Gn  (0)  via  an  inverse  Fourier 
transform  [79] ,  which  is  given  by  R{xi  —  X2)  =  exp  where  X1  —  X2  =  woS. 

This  is  true  since,  from  Eqs.  6.13  and  6.19,  p(l  -p)"^  =  V^/4.  In  other  words, 
the  correlation  length  of  a  logarithmic  incoherent  soliton  remains  invariant  during 
propagation  and  it  is  equal  to  that  of  the  exciting  incoherent  source.  Moreover,  as  in 
the  fully  coherent  case  [80],  the  dynamics  of  these  two-dimensional  solitons  (during 
compression  or  expansion,  i.e.  when  p  and  q  deviate  from  the  prescribed  values  of 
Elqs.6.19  and  6.20)  can  also  be  described  in  closed  form.  We  would  like  to  emphasize 
that  by  its  very  nature  the  logarithmic  nonlinearity  leads  to  an  in  nite  set  of  modes 
which  are  in  turn  related  through  a  Poisson  distribution.  Of  course  this  is  also  true 
for  the  ln(l  -t-  In)  nonlinearity  provided  that  the  maximum  intensity  at  the  center 
of  the  beam  is  much  greater  than  If  In  particular,  computer  simulations  show  that 
intensity-wise  the  ln(/7v)  model  nicely  approximates  the  ln(l  -f  In)  nonlinearity  as 
long  as  the  intensity  ratio,  max(7jv)  >  100.  However  for  relatively  low  intensity 
ratios  the  mode  occupancy  (in  the  ln(l  -b  In)  model)  is  not  exactly  Poissonian  and 
the  correlation  length  (even  though  at  in  the  middle)  tends  to  increase  at  the  low 
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intensity  boundaries  of  the  beam.  This  is  in  full  agreement  with  the  results  we 
have  previously  obtained  [83]  for  incoherent  photorefractive  spatial  solitons  where 
An  {1  +  I m)~^-  Finally  at  very  high  intensity  ratios,  the  correlation  length  becomes 
again  constant  across  the  Gaussian  beam  as  suggested  by  Eq.  6.21. 

6.6  Results  and  discussion 

As  an  example,  let  us  consider  a  logarithmically  saturable  nonlinear  medium  with  a 
ln(l  +  In)  nonhnearity,  no  =  2  and  n2  =  10“^.  The  free  space  optical  wavelength  is 
taken  to  be  Aq  =  0.5  fim  and  the  width  of  the  angular  power  spectrum  of  the  source 
is  ^0  =  0.258°.  Fig.6.1(a)  shows  stationary  propagation  of  a  1-D  incoherent  soliton 
when  r  =  10^  and  xq  ~  18  /im  (intensity  FWHM  is  ~  30  /xm)  as  obtained  numer¬ 
ically  from  the  coherent  density  method.  Fig.6.1(b)  on  the  other  hand,  illustrates 
the  intensity  pro  le  of  the  rst  four  modes  involved  in  this  Gaussian  soliton  as  ob¬ 
tained  from  beam  propagation  methods.  Evidently,  these  individual  modes  remain 
invariant  as  a  function  of  distance.  As  expected,  the  smn  of  the  individual  modes 
provides  the  intensity  pro  le  of  Fig.6.1(a).  These  results  are  in  excellent  agreement 
with  our  theoretical  analysis. 
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6.7  Conclusions 


In  conclusion  we  have  shown  that  incoherent  spatial  solitons  are  possible  in  log¬ 
arithmically  saturable  nonlinear  media.  These  solitons  can  exist  as  long  as  their 
mode-occupancy  function  obeys  a  Poisson  distribution.  We  have  found  that  two 
approaches,  that  is  the  dynamic  coherent  density  description  as  well  as  the  self- 
consistent  multimode  method  lead  in  this  case  to  exactly  the  same  results.  Two 
dimensional  (circular  and  elliptical)  multimode  spatial  solitons  have  also  been  ob¬ 
tained  in  closed  form. 
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Chapter  7 


Dynamics  of  incoherent 
self-trapped  beams 


In  this  chapter,  using  the  coherent  density  approach,  we  study  the  propagation  dy¬ 
namics  of  incoherent  bright  and  dark  beams  in  biased  photorefractive  crystals.  We 
show  that  under  appropriate  initial  conditions,  bright  as  well  as  dark-like  incoherent 
quasi-solitons  can  be  established  in  this  material  system.  Our  numerical  simulations 
demonstrate  that  the  coherence  properties  of  these  beams  can  be  signi  cantly  af¬ 
fected  by  the  self-trapping  process. 


105 


7.1  Introduction 


In  Chapter  3,  a  theory  based  on  the  so-called  coherent  density  approach  has  been 
developed  and  good  agreement  was  found  between  the  experimentally  observed  be¬ 
havior  and  the  theoretical  results  of  this  model  [79] .  In  Chapter  4,  using  the  coherent 
density  approach,  an  exact  analytical  solution  has  been  obtained  for  (1+1)  D  in¬ 
coherent  spatial  solitons  in  saturable  nonlinear  media  of  the  logarithmic  type  [82]. 
We  have  also  developed  an  alternative  method  which  is  capable  of  identifying  static 
incoherent  spatial  soliton  states  [83].  This  latter  model,  seeks  in  a  self-consistent 
fashion  multimode  sohton  solutions  whose  total  intensity  can  be  obtained  via  in¬ 
tensity  superposition  of  all  the  modes  guided  in  the  nonUnearly  induced  waveguide. 
In  this  manner,  we  have  found  (1+1)  D  photorefractive  incoherent  spatial  solitons, 
the  conditions  necessary  for  their  existence  as  well  as  their  spatial  coherence  char¬ 
acteristics  [83].  In  Chapter  5,  by  using  this  latter  approach,  we  have  shown  that 
incoherent  solitons  are  possible  in  non-instantaneous  Ker-like  nonlinear  media.  In 
Chapter  6,  on  the  other  hand,  it  is  demonstrated  that  the  log  type  incoherent  soli¬ 
tons  can  also  be  described  by  using  self-consistent  multimode  method  [114].  These 
studies  have  shown  that  incoherent  spatial  bright  solitons  are  indeed  feasible  in  any 
non-instantaneous  nonhnear  media.  Yet  at  this  point,  there  are  still  several  other 
issues  that  merit  further  investigation.  For  example,  it  is  not  at  all  clear  whether 
incoherent  dark  sohtons  are  also  possible.  Moreover,  it  is  of  interest  to  know  how  the 
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coherence  properties  or  uctuation  statistics  of  these  beams  are  in  uenced  during 
propagation  by  the  photorefractive  nonlinearity. 

In  this  chapter,  we  study  the  dynamics  of  incoherent  bright  and  dark  beams  in  bi¬ 
ased  photorefr actives,  which  represent  an  example  of  a  non-instantaneous  saturable 
nonlinear  material  system.  We  employ  the  dynamic  coherent  density  approach  and 
numerically  solve  the  underlying  nonlinear  Schrodinger-like  integro-differential  equa¬ 
tion.  We  predict  that  under  odd  initial  beam  conditions,  a  dark-like  incoherent 
spatial  soliton  will  emerge  and  after  a  certain  distance  it  takes  the  form  of  a  gray 
structure.  As  in  the  case  of  coherent  sohtons,  incoherent  gray  soliton  pairs  (dou¬ 
blets)  and  triplets  are  also  possible  under  even  and  odd  initial  beam  conditions 
respectively.  In  all  cases  (dark  and  bright),  we  show  that  the  photorefractive  nonli- 
narity  signi  cantly  alters  the  coherence  properties  of  these  beams.  For  bright  beams, 
we  nd  that  the  correlation  length  across  the  beam  nally  evolves  according  to  the 
predictions  of  Ref.  [83].  Moreover,  we  emphasize  that,  our  results  are  of  a  gen¬ 
eral  nature,  i.e.  they  can  be  used  to  describe  incoherent  beam  dynamics  in  any 
non-instantaneous  nonlinear  medium. 
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7.2  Evolution  dynamics  of  incoherent  beams  in 
biased  photorefractive  media 


7.2.1  Evolution  equation 


In  a  standard  biased  Strontium  Barium  Niobate  (SBN)  photorefraotive  crystal  arrange¬ 
ment  [23], [24],  the  normalized  intensity  In  =  I /Id  of  an  incoherent  beam  obeys  the 
following  nonlinear  integro-differential  equation  (see  Chapters  3  and  4): 


dz  dx 


+ 


2k  dx"^ 


-^nlvs^Eo  (1  +  p) 


/ 


1  +  In  {x,  z) 


=  0 


(7.1) 


where 


and  at  z  =  0 


-f-OO 

In{x,z)=  J  \f  {x,z,9)f  d9 

—  OO 


f{z  =  0,x,9)  =  {r,pf^G]i\9)(t>,{x) 


(7.2) 


(7.3) 


In  the  above  equations,  /  represents  the  so-called  coherent  density,  0  is  an  angle 
(in  radians)  with  respect  to  the  z  axis,  k  —  kong,  and  ko  =  27r/Ao.  Id  is  the  dark 
irradiance,  Ug  the  extraordinary  refractive  index  and  rss  the  electrooptic  coefficient 
involved.  Eq  =  ^.V/W  where  V  is  the  applied  bias  and  W  the  x-width  of  the  crystal. 
Gn{9)  is  the  normalized  angular  power  spectrum  of  the  incoherent  source  and  <f)o{x) 
is  the  input  complex  spatial  modulation  function.  Moreover,  r  is  the  maximum 
intensity  of  a  bright  incoherent  beam  whereas  p  is  the  normalized  intensity  of  a 
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dark  beam  at  a:  ±oo.  Equation  7.3  restates  the  fact  that  the  coherent  density 

has  been  initially  weighted  with  respect  to  the  source  angular  power  spectrum. 


7.2.2  Coherence  length 


The  coherence  properties  of  these  evolving  partially  incoherent  beams  can  be  ob¬ 
tained  through  a  modi  ed  version  of  the  Van  Cittert-Zernike  theorem  [74],  particu¬ 
larly  suited  for  the  coherent  density  approach.  The  complex  coherence  factor  /ijj  is 
obtained  from 


Mi2  {xi,X2) 


_ 1 _ 

y/lNlxi,z)lN{x2,z) 

-foo 

X  j  fixi,z,  6)  f*  (X2, 2, 9)  exp  [ikO  (xi  -  a;2)]  d9 


(7.4) 


In  general,  the  complex  coherence  factor  can  be  an  arbitrary  function  of  both  Xi 
and  X2.  Thus,  it  is  often  more  convenient  to  describe  the  coherence  properties  of  an 
optical  beam  through  the  so-called  correlation  length  Ic  which  is  de  ned  as  follows 

lc{x)=  j  \fii2(x,x)f  dx  .  (7.5) 
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7.3  Dynamics  of  incoherent  bright  and  dark  beams 


and  their  coherence  properties 

As  an  example  let  us  consider  a  biased  SBN:75  crystal  with  rie  =  2.3,  =  1022 

pm/V,  Ao  =  488  nm,  W  =  6  mm,  similar  to  that  used  in  Ref.  [35].  Let  us  also 
assume  that  at  the  input,  the  source  uctuations  obey  a  stationary  random  process 
(i.e.  fj,i2  depends  only  on  X1—X2)  and  that  the  normalized  angular  power  spectrum  of 
the  incoherent  source  is  Gaussian,  i.e.  Gff{d)  =  exp(— where  Oq  is 

associated  with  its  angular  width  [79].  Equations  7.1  and  7.2  are  numerically  solved 
by  appropriately  discretezing  them  in  6  and  by  using  beam  propagation  methods. 
In  turn,  the  correlation  length  is  obtained  from  Eqs.  7.4  and  7.5. 

7.3.1  Incoherent  bright  beeuns 

Let  us  rst  consider  an  incoherent  Gaussian  bright  beam  whose  spatial  modulation 
function  is  given  by  <^o  (®)  =  6xp[— a:^/(2xo)]  in  which  case  In  =  rexp(— s^/xq) 
where  Xq  is  related  to  the  intensity  FWHM.  At  the  input,  the  intensity  FWHM  of 
this  beam  is  30  pm,  and  its  maximum  normalized  intensity  is  r  =  3.  The  width  of 
the  angular  power  spectrum  is  taken  here  to  be  60  =  0.55°.  Since  we  have  assumed 
stationary  source  uctuations,  it  can  then  be  readily  shown  that  at  2  =  0,  Pi2{S)  = 
exp  [-5^/5o]  where  6  =  X1-X2  and  60  -  2/{k0o).  Thus,  the  incoherent  beam  at  2  = 
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0  exhibits  a  constant  coherence  length  Ic  =  \/2^/  (kOo)  which  is  independent  of  the 
spatial  modulation  function  ^o(®)  or  In{x,z  =  0).  For  the  above  given  parameters, 
the  coherence  length  of  the  source  is  approximately  8.8  fim.  It  can  be  shown,  that 
during  linear  diffraction  the  FWHM  of  this  bright  beam  increases  to  165  nm  after  1 
cm.  The  coherence  length  lc{x,  z)  of  this  Gaussian  beam  also  increases  with  z  while 
remaining  constant  across  x.  Note  that,  this  latter  feature  is  particular  to  Gaussian- 
Schell  sources.  If  the  input  complex  envelope  (^)  ^  Gaussian,  the  correlation 
length  Ic  is  no  longer  constant  across  the  beam  during  diffraction.  On  the  other 
hand,  by  activating  the  photorefractive  nonlinearity,  by  applying  a  bias  voltage,  the 
coherence  properties  change  dramatically.  For  example,  by  setting  the  applied  bias 
to  550  V,  this  incoherent  bright  beam  propagates  almost  undistorted  and  behaves 
like  a  quasi-soliton  as  shown  in  Fig.7.1.  This  incoherent  quasi-soliton  tends  to 
slightly  oscillate  in  an  effort  to  adjust  to  its  new  statistical-phase  environment.  In 
this  case,  the  coherence  length  remains  approximately  constant  around  the  center  of 
the  beam  whereas  it  increases  at  the  margins  as  depicted  in  Fig.7.2  within  a  spatial 
window  of  ±40  fim.  This  behavior  is  in  good  agreement  with  the  static  predictions  of 
Ref.  [83] .  The  existence  of  this  bright  quasi  soliton  can  be  understood  on  the  basis  of 
the  multimode  approach  [83].  The  increase  of  the  coherence  length  at  the  margins  of 
the  beam  can  be  explained  from  the  fact,  that  away  from  the  center,  only  the  higher 
order  modes  contribute  to  the  coherence.  Our  numerical  simulations  also  show  that 
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further  increase  of  the  bias  voltage  sets  the  beam  into  cycles  of  compression  and 
expansion.  In  this  latter  case,  the  coherence  length  decreases  at  the  center  of  the 
beam  during  compression  and  it  increases  when  the  beam  expands. 

7.3.2  Incoherent  dzirk  bezims 

As  a  second  example,  we  consider  the  dynamics  and  coherence  properties  of  an 
incoherent  dark-like  beam  under  odd  and  even  initial  conditions.  In  particular,  let 
the  spatial  modulation  function  at  the  input  be  (j>o{x)  =  tanh(x/a:o)  under  odd 
initial  conditions,  and  {x)  =  [1  -  (1  -  e^)  sec/i^  (x/xq)]^^^  under  even.  Here  xo  is 
associated  with  the  intensity  FWHM  of  this  beam  and  <  1.  In  this  example,  the 
intensity  FWHM  of  these  dark  beams  is  taken  to  be  10  /xm  and  their  normalized 
backgroimd  intensity  is  p  =  3.  Let  again,  be  0.55°  with  a  source  coherence  length 
Ic  8.8  fitn.  It  is  important  to  note,  that,  during  linear  diffraction  the  coherence 
length  of  an  incoherent  dark  beam  becomes  a  fimction  of  x  which  implies  that  eld 
uctuations  are  no  longer  stationary.  Moreover,  unlike  their  coherent  counterparts, 
incoherent  dark  beams  tend  to  quickly  disappear  during  diffraction,  without  leaving 
any  trace  of  their  previous  existence.  For  example,  after  1  cm  of  propagation,  an 
odd  dark  beam  whose  initial  FWHM  is  10  /xm  will  diffract  to  138  pm,  and  its 
grayness  will  be  ~  90%,  and  the  dark  notch  almost  disappears.  Our  numerical 
simulations  show  that  photorefractive  nonlinearity  plays  again  an  important  role. 
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With  the  application  of  ~  -300  V,  after  a  certain  distance  of  propagation,  an 
initially  odd  beam  becomes  a  gray-like  incoherent  quasi-soliton  whose  FWHM  and 
grayness  oscillate  around  constant  values.  The  intensity  and  the  coherence  length 
of  this  incoherent  beam  is  depicted  in  Figures  7.3  and  7.4.  Fig.7.3  shows  that  even 
though  the  beam  itself  exhibits  a  somewhat  oscillatory  behavior,  its  spatial  FWHM 
varies  between  9  /xm  and  5.6  fim  and  its  grayness  between  56  %  and  31  %  over 
a  distance  of  2  cm.  Furthermore,  its  coherence  length,  Fig.7.4,  becomes  higher 
(by  r\j  20  %)  within  the  dark  notch  with  a  mild  depression  at  the  center.  Note 
that  a  gray-like  fundamental  incoherent  dark  quasi-soliton  was  also  suggested  by 
Hasegawa  [85], [86], [87]  by  averaging  the  phase  of  all  the  quasi-particles  involved  in 
the  Vlasov  transport  equation.  Our  analysis  on  the  other  hand,  shows  that  the  phase 
is  extremely  important  (even  more  than  the  normalized  intensity  p),  i.e.  a  single 
dark  incoherent  soliton  is  only  possible  provided  that  initially  a  7r-phase  shift  is 
imposed  on  the  wavefront.  Fig.7.5  shows,  what  happens  to  this  beam  if  we  increase 
the  bias  voltage  to  —550  V.  In  this  case,  the  odd  incoherent  dark  beam  leads  to 
a  triplet.  This  differs  from  coherent  triplets  [9], [34], [108]  since  the  central  soliton 
is  gray-hke.  Fig.7.6  shows  that  the  coherence  length  increases  again  by  ~  20  % 
around  the  soliton  regions.  If  the  incoherent  beam  at  the  input  is  even,  a  gray 
soliton  doublet  occurs.  Our  simulation  shown  in  Fig.7.7,  demonstrates  that  this 
is  true  when  the  bias  voltage  is  —550  V,  the  intensity  FWHM  is  initially  10  pm 
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and  <C  1.  Fig.  7.8  shows  the  correlation  length  associated  with  these  incoherent 
doublets. 

Depending  on  the  initial  conditions,  an  even  or  odd  number  of  incoherent  dark¬ 
like  structures  can  appear  at  higher  voltage  values.  From  these  1-D  results  one  may 
speculate  that  2-D  incoherent  vortex  solitons  may  also  be  possible  in  photorefractives 
and  in  other  non-instantaneous  nonlinear  media.  Furthermore,  we  emphasize  that 
the  behavior  of  these  bright/dark  incoherent  quasi-soHtons  is  fundamentally  different 
from  that  of  their  coherent  counterparts.  For  example,  unlike  a  coherent  gray  sohton, 
the  fundamental  incoherent  dark  soliton  does  not  exhibit  a  transverse  velocity  in 
spite  of  its  grayness. 

7.4  Conclusions 

In  conclusion,  we  have  shown  that  under  appropriate  initial  conditions,  bright  as 
well  as  dark-like  incoherent  quasi-solitons  can  form  in  biased  photorefractives.  Using 
numerical  simulations  we  have  found  that  incoherent  gray  sohton  pairs  as  well  as 
triplets  are  possible  under  even  or  odd  initial  conditions.  The  coherence  properties 
of  these  beams  were  found  to  be  signi  cantly  affected  during  propagation  by  the 
photorefractive  nonhnearity.  Relevant  examples  have  been  provided. 
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Figure  7.5:  Formation  of  an  incoherent  soliton  triplet  when  p=3  and 
the  bias  is  -550  V. 
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Figure  7.6:  The  coherence  length  4  of  the  triplet  shown  in  Fig.7.5  as  a  function  of 
distance. 


120 


2 


Chapter  8 

Dark  incoherent  soliton  splitting 
2ind  phcuse-memory  effects 

In  this  chapter,  the  effects  of  incoherence  on  the  evolution  of  incoherent  dark  soliton 
doublets  are  investigated  both  theoretically  and  experimentally.  We  show  that  the 
dynamics  of  these  incoherent  self-trapped  entities  are  associated  with  strong  phase- 
memory  effects  which  are  otherwise  absent  in  the  linear  regime. 
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8.1  Introduction 


In  the  last  chapter,  we  have  shown  that  incoherent  dark  solitons  may  be  possible  in 
biased  photorefractives  [109] .  Based  on  this  prediction,  incoherent  dark  planar  and 
2-D  dark  solitons  (  vortices  )  were  observed  in  a  SBN:60  crystal  [llOj.  This  was 
achieved  by  employing  the  photorefractive  self-defocusing  nonlinearity  associated 
with  screening  sohtons  [23], [24].  As  predicted  in  [109],  the  incoherent  dark  solitons 
were  found  to  be  gray.  Moreover,  these  dark  incoherent  solitons  were  efficiently  ex¬ 
cited  provided  that  an  initial  7r-phase  ip  was  imposed  on  the  incoherent  wavefront. 
Incoherent  Y-soliton  splitting  was  also  predicted  in  the  last  chapter.  Nevertheless, 
thus  far  we  did  not  consider  the  behavior  of  these  incoherent  soliton  doublets,  for 
different  degrees  of  coherence. 

In  this  chapter,  we  report  the  rst  experimental  observation  of  dark  incoher¬ 
ent  sohton  Y-splitting  in  a  noninstantaneous  self-defocusing  nonlinear  medium  as 
predicted  in  Chapter  7.  The  evolution  of  these  incoherent  soUton  doublets  is  then 
systematically  investigated  as  a  function  of  their  coherence,  both  theoretically  and 
experimentally.  Surprisingly,  we  nd  that  over  a  wide  range  of  parameters,  the  Y- 
splitting  is  approximately  the  same  irrespective  of  coherence.  Moreover,  we  show 
that  the  dynamical  behavior  of  this  incoherent  Y-splitting  process  is  associated  with 
strong  phase-memory  effects  which  are  otherwise  absent  in  the  linear  regime.  In 
other  words,  we  show  that  dark  incoherent  self-trapped  entities  (dark  incoherent 
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solitons)  are  characterized  by  a  strong  memory  effect  that  lasts  throughout  propa¬ 
gation  and  governs  their  propagation  behavior  (single  soliton  versus  Y-sohton  spht- 
ting,  etc.).  This  is  in  sharp  contrast  to  all  known  so  far  about  linear  propagation  of 
incoherent  beams,  in  which  all  phase  information  is  fully  washed  out  after  a  nite 
distance  [111], [112]. 


8.2  Coherent  and  incoherent  Y-juntion  solitons 


8.2.1  Diffraction  of  coherent  and  incoherent  beams 


Our  experiments  were  carried  out  in  SBN;60  crystals.  For  this  reason,  here  we  use 
the  (1+1)  D  saturable  nonlinearity  of  the  form  1/(1  +  I)  [23],  [24]  so  as  to  make 
direct  comparisons  with  experiment.  We  emphasize,  however,  that  our  results  hold 
for  any  noninstantaneous  nonlinearity  that  can  give  rise  to  dark  soUtons.  In  this 
material  system  (photorefractives),  the  normalized  intensity  1^  =  l/Id  (where  Id 
is  the  dark  irradiance)  of  the  incoherent  dark  beam  evolves  according  to  following 
normalized  nonlinear  integro-differential  equation  [79], [82], [109]: 


+  +  ^ _ I _ 


=  0 


(8.1) 


where 
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lNisX)=  I  \f{s,c,0)f 
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and  at  C  =  0,  the  coherent  density  /  is  given  by 

f{C  =  o,s,e)  =  p^^^Gfie)<f>,{s)  .  (8.3) 

In  the  above  equations,  we  have  used  the  following  normalized  coordinates;  C  = 
z/{kxl)  and  s  =  x/xq  where  Xo  is  an  arbitrary  spatial  scale  associated  with  the 
intensity  FWHM  of  the  beam.  Moreover,  a  =  kxo9,  =  {k‘^XQ/2)n\T^^  |£o|  (1  +  p) 
where  6  represents  an  angle  (in  radians)  with  respect  to  the  z  axis,  k  =  koUe  is  the 
wavenumber,  ko  =  27r/Ao,  Ug  is  the  extraordinary  refractive  index  of  the  material, 
and  r33  is  the  electrooptic  coefficient  involved.  Eq  =  —V/W  is  the  value  of  the  space 
charge  eld  at  a:  ±oo  where  V  is  the  reverse  applied  bias  and  W  the  a;- width 

of  the  crystal.  Gn{9)  is  the  normalized  angular  power  spectrum  of  the  incoherent 
source  and  ^o(^)  i®  input  complex  spatial  modulation  fimction.  In  this  study, 
Gi\f{9)  is  assumed  to  be  Gaussian,  i.e.  Gn{9)  =  exp(— where  9o  is 

associated  with  the  width  of  angular  power  spectrum.  Finally,  p  is  the  normalized 
intensity  of  the  dark  beam  at  x  — >  ±oo.  Here,  as  usual,  we  assume  that,  the  beam  at 
the  input  obeys  a  stationary  random  process.  In  general,  the  coherence  properties 
of  these  beams  can  be  followed  using  a  version  of  the  Van  Cittert-Zernike  theorem 
as  in  the  last  chapter.  The  coherence  length  of  the  beam  at  z  =  0  can  be  readily 
obtained  from  Gn{9)  and  it  is  given  by  Ig  =  ik9o).  When  the  beam  is  fully 

coherent  (^o  =  0);  the  coherence  length  of  the  beam  becomes  in  nite,  i.e.  Ig  oo. 
In  this  case,  the  system  of  Eqs.  8.1  and  8.2  collapse  to  a  standard  single  differential 
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equation  given  in  Chapter  2  [23], [24]. 

Before  we  present  our  experimental  results,  it  may  prove  bene  cial  to  rst  discuss 
the  behavior  of  such  incoherent  dark  beams  from  a  theoretical  point  of  view.  As  in 
the  experiment,  let  us  consider  a  biased  SBN:60  crystal  with  rie  =  2.3,  =  250 

pm/V,  Ao  =  514  nm,  W  =  5.3  mm  [110].  We  let  the  spatial  modulation  fimction 
at  the  input  be  (f)o{x)  =  tanh(x/xo)  under  odd  initial  conditions,  and  <f>o{x)  = 
[1  —  e^sech^  {x/xo)f^^  under  even.  The  quantity  de  nes  the  beams  grayness. 
Throughout  this  work,  we  assume  that  at  z  =  0,  «  1  (almost  black  even  dark 

beams).  The  input  intensity  FWHM  of  the  even  and  odd  dark  beams  is  taken 
here  to  be  25  pm  as  shown  in  Fig.8.1(a).  Moreover,  the  normalized  background 
intensity  is  p  =  3.  First,  we  consider  linear  diffraction  of  coherent  and  incoherent 
dark-beams  imder  odd  and  even  initial  conditions.  Figs.8.1(b)  and  (c)  show  the 
diffracted  intensity  pro  les  of  coherent  odd  and  even  dark  beams  respectively  after 
~  12  mm  of  propagation.  In  this  case,  the  intensity  FWHM  of  the  odd  dark  beam 
at  the  output  is  ~  42  pm  whereas  that  of  the  even  is  ~  76  pm.  It  is  important 
to  note  that  after  diffraction,  the  intensity  of  the  odd  coherent  beam  is  always  zero 
at  the  center  whereas  that  of  the  even  is  gray-like.  Fig.  8.1(d),  on  the  other  hand, 
demonstrates  how  an  odd  or  even  incoherent  dark  beam  will  diffract  after  12  mm  of 
propagation  when  at  the  input  0o  —  5  mrads  or  =  17  pm.  This  latter  gure  shows 
that,  the  intensity  pro  les  of  the  odd  and  even  incoherent  dark  beams  are  almost 
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identical  with  an  output  FWHM  of  ^  100  ^im.  Simulations  suggest  that  the  same 
also  applies  for  the  lc{x)  curves  corresponding  to  these  two  cases  as  shown  in  Fig. 
8.1(e).  Thus,  from  diffraction  data  alone,  it  is  extremely  difficult  to  distinguish  an 
odd  dark  beam  from  an  even  one.  In  other  words,  the  randomly  changing  speckled 
structure  of  an  incoherent  beam  leads  to  a  loss  of  phase  memory.  Therefore,  as 
a  result  of  this  phase  washing  effect,  a  sufficiently  incoherent  dark  beam  diffracts 
approximately  the  same  way  regardless  of  the  phase  information  initially  imposed  on 
it.  An  important  distinction  between  diffraction  of  a  coherent  and  an  incoherent  dark 
beam  comes  from  the  structure  of  their  background.  Figure  8.1  clearly  demonstrates 
that  a  diffracted  coherent  dark  beam  involves  intensity  ripples  in  its  background. 
These  oscillations  tend  to  disappear  in  the  case  of  an  incoherent  beam  as  a  result 
of  its  speckled  structure. 

When  on  the  other  hand  the  nonlinearity  is  activated,  the  dynamics  of  these  in¬ 
coherent  beams  depend  on  initial  phase  information.  As  previously  predicted,  gen¬ 
eration  of  a  single  incoherent  dark  (which  is  in  reality  gray)  beam  or  a  higher-order 
triplet  requires  a  7r-phase  shift.  Conversely,  starting  from  even  initial  conditions,  an 
incoherent  gray  soliton  pair  or  Y-soliton  sphtting  can  be  obtained  [109].  In  other 
words,  in  the  presence  of  nonlinearity,  an  incoherent  dark  beam  tends  to  remember 
its  origins  and  identity,  i.e.  a  phase-memory  effect  is  established.  Thus,  the  beam 
starts  to  behave  in  a  quasi-coherent  fashion  [34]. 
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Figure  8. 1 :  (a)  Intensity  profile  of  a  25  /tm  odd  or  even  dark  beam  at  the  input. 
Diffraction  of  an  (b)  odd  coherent  dark  beam,  (c)  even  coherent  dark  beam,  (d) 
incoherent  odd  or  even  dark  beam  with  4  =  7.3  fim  after  12  mm  of  propagation,  (e)  4 
in  jMn  as  a  function  of  x  for  the  odd  (dashed  curve)  and  even  (solid  curve)  diffracted 
incoherent  dark  beam  shown  in  (d). 
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8.2.2  Experiment  results 


Experiments  with  an  amplitude  notch  (even  initial  conditions)  are  performed  by 
using  both  coherent  and  spatially  incoherent  light  sources  for  comparison.  Details 
regarding  the  coherent  dark  soliton  experiments  can  be  found  in  Ref.  [34].  The  laser 
used  is  a  Ar  ion  laser  (Aq  =  514  nm).  An  arti  cial  background  dark  irradiance  is 
provided  by  uniformly  illiuninating  the  entrance  face  of  the  crystal  (SBN:60)  along 
the  ordinary  axis.  The  maximum  intensity  ratio  (at  the  tails)  of  the  dark  beam  with 
respect  to  dark  irradiance  is  approximately  1.5.  The  dark  beam  is  also  broad  enough 
to  cover  entire  input  face  of  the  crystal.  For  the  incoherent  case,  a  rotating  dilfuser 
is  employed  to  provide  random  phase  uctuations  across  the  beam  [35]  ,[36]  ,[110] .  In 
this  case,  we  generate  a  dark  notch  on  a  broad  partially  spatially  incoherent  beam 
with  controllable  degree  of  coherence.  The  experimental  arrangement  is  the  same 
as  that  in  Ref.  [110],  except  that  the  phase  mask  is  now  replaced  by  an  amplitude 
mask.  Incoherent  Y-jimction  solitons  are  generated  and  then  compared  with  the 
coherent  ones.  Fig.8.2  shows  typical  experimental  results.  When  the  dark  beam  is 
coherent,  it  diflEracts  from  a  FWHM  of  25  nm  (left)  to  about  58  after  ~  12  mm  of 
propagation  (middle)  when  no  nonlinearity  is  present.  Note  that,  with  the  exception 
of  the  dark  notch  FWHM  (which  from  simulations  is  expected  to  be  ~  76  yum),  its 
intensity  structure  is  in  agreement  with  Fig.8.1(c).  The  discrepancy  in  FWHM 
is  attributed  to  the  fact  that  the  re  ection  from  the  metallic  wire  introduces  a 
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quadratic  phase  which  is  not  accounted  in  our  simulations.  After  applying  a  voltage 
of  -350y  (negative  relative  to  the  c-axis),  the  dark  amplitude  notch  evolves  into 
a  pair  of  gray  solitons  (right).  The  second  and  third  rows  of  this  gure  depict  the 
same  data  when  the  dark  beam  is  incoherent.  The  right  column  of  the  gure  was 
obtained  ai  V  —  —350  V  and  with  an  input  FWHM  of  25  fim.  As  seen  in  Fig.8.2, 
the  grayness  of  the  soliton  pair  increases  as  the  incoherence  of  the  beam  increases. 
Nevertheless,  the  spacing  of  these  two  sohtons  at  the  crystal  output  face  is  about 
the  same  for  a  varying  degree  of  coherence. 

8.2.3  Simulation  results 

These  experimental  results  are  now  compared  with  numerical  simulations.  By  keep¬ 
ing  in  mind  that  in  the  experiment,  the  input  speckle  size  of  the  incoherent  beams 
is  30  fim  for  Fig.8.2(b)  and  15  /xm  for  Fig.8.2(c)  and  by  considering  their  diffraction 
behavior,  we  estimate  that  the  width  of  the  angular  power  spectrum  in  these  two 
cases  is  ~  3.5  mrads  and  5.2  mrads  respectively.  The  simulation  shown  in  Fig.8.3(a) 
demonstrates  how  a  coherent  soliton  doublet  forms  from  a  25  fim  even  dark  beam 
after  12  mm  of  propagation  when  V  =  —450  V.  For  the  same  bias  voltage  and 
initial  beam  width,  the  intensity  pro  le  of  an  incoherent  doublet  after  12  mm  of 
propagation  is  shown  in  Fig.8.3(b)  when  9o  =  3.5  mrads  or  4  =  25  /xm.  Fig.  8.3(c) 
depicts  similar  data  when  Oq  =  5.2  mrads  (Ic  =  17  /xm)  and  again  V  =  —450  V. 
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Figure  8.2:  Experimental  observation  of  coherent  and  incoherent  Y-splitting;  (a) 
coherent  dark  beam;  (b)  (c)  incoherent  dark  beam  with  an  average  speckle  size  of  30 
im  and  15  fim  respectively.  The  first  column  depicts  the  input  intensity,  the  second 
one  diffraction  data,  and  the  third  one  Y-splitting  at  -350  V.  In  all  the  cases  the 
intensity  FWHM  of  the  beam  at  the  input  is  25  fim. 
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Both  gures,  8.3(b)  and  (c),  were  obtained  by  numerically  solving  Eqs.  8.1-8.3  as 
done  in  previous  chapters.  In  agreement  with  the  experiment,  Fig. 8. 3  demonstrates 
that  the  doublet  becomes  grayer  as  the  incoherency  increases.  Surprisingly,  for  this 
range  of  parameters,  both  theory  (Fig.8.3)  and  experiment  (Fig.8.2)  suggest  that  the 
Y-splitting  angle  or  the  doublet  separation  does  not  depend  strongly  on  the  degree 
of  coherence.  To  further  understand  this  Y-splitting  process,  we  carried  out  another 
set  of  simulations.  In  this  latter  set,  the  intensity  FWHM  of  the  even  dark  beam 
was  chosen  to  be  10  /xm  (in  order  to  accelerate  splitting  process)  and  varied  from 
oo  down  to  3.4  /xm.  Figs.8.4(a)-(c)  were  obtained  for  the  same  initial  conditions 
and  bias  voltage  {V  =  -2400  V)  after  12  mm  of  propagation  for  different  degrees 
of  coherence.  Even  in  this  case,  the  splitting  is  relatively  insensitive  to  6q.  This  is 
by  itself  very  interesting  considering  the  range  in  which  Ic  varies.  This  is  another 
manifestation  of  the  phase  memory  effect  discussed  earlier.  As  the  incoherency 
of  the  dark  beam  increases,  a  higher  bias  voltage  is  required  to  establish  a  doublet. 
Fig.8.4(d)  shows  Y-splitting  of  a  10  /xm  even  incoherent  dark  beam  after  12  mm, 
when  V  =  —4000  V  and  Ic  =  3.4  /xm.  Finally,  at  lower  1^  s,  the  doublet  practically 
disappears  (because  of  its  gra3mess)  and  the  splitting  angle  is  further  reduced. 
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Figure  8.3:  Intensity  profile  of  a  soliton  doublet  a.tz  =  12  mm  when  the  external  bias  is 
-450  V  and  the  beam  is  (a)  coherent  or  incoherent  with  (b)  4  =  jum,  (c)  Ic-  17  fJfn. 
In  all  cases  the  initial  intensity  FWHM  of  the  beam  is  25  jum. 
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Normalized  Intensity 


X  (|im) 


Figure  8.4;  Intensity  profile  of  a  soliton  doublet  at  z  =  72  mm  when  the  external  bias  is 
-2400  V  and  the  beam  is  (a)  coherent  or  incoherent  with  (b)  4  =  9.3  jxm,  (c)  4  =  7.3 
fm.  Same  information  when  the  external  bias  is  -4000  V  and  l^  =  3.4  ^im.  In  all  cases 
the  initial  intensity  FWHM  of  the  beam  is  70  jum. 
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8.3  Conclusions 


In  conclusion,  incoherent  dark  soliton  Y-splitting  has  been  demonstrated  experi¬ 
mentally.  Using  the  coherent  density  approach  we  have  shown  that  the  evolution 
of  incoherent  dark  solitons  in  non-instantaneous  nonlinear  media  is  associated  with 
strong  phase-memory  effects  which  are  otherwise  absent  in  the  linear  regime.  The 
higher-order  behavior  of  these  dark  beams  has  been  compared  under  the  same  ini¬ 
tial  conditions  but  for  different  degrees  of  coherence.  It  was  found  that  over  a  wide 
range  of  parameters,  the  Y-splitting  is  approximately  the  same  irrespective  of  spa¬ 
tial  coherence.  Experimental  observations  are  in  good  agreement  with  theoretical 
predictions. 
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Chapter  9 


Theory  of  Incoherent  Dark 
Solitons 


In  this  chapter,  we  formulate  the  theory  of  incoherent  dark  spatial  solitons  in  non- 
instantaneous  self-defocusing  nonlinear  media.  We  nd  that  the  basic  modal  con¬ 
stituents  of  these  incoherent  dark  soliton  entities  are  radiation  modes  as  well  as 
bound  states.  Our  results  explain  for  the  rst  time  why  incoherent  dark  solitons  are 
in  fact  gray  and  why  a  transverse  7r-phase  ip  can  facilitate  their  observation. 
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9.1  Introduction 


The  numerical  results  of  Chapter  7  suggested  that  incoherent  dark  quasi-solitons 
can  be  effectively  excited  in  self-defocusing  (reverse  biased)  photorefractive  crystals, 
provided  that  at  the  origin,  a  tt  phase  jump  is  imposed  on  the  incoherent  wavefront 
[109].  Even  more  importantly,  unlike  their  coherent  counterparts  [113],  these  dark 
incoherent  solitons  were  always  found  to  be  gray!  The  gray  character  of  these 
solutions  is  in  qualitative  agreement  with  some  earlier  predictions  of  random-phase 
envelope  solitons  made  by  Hasegawa  two  decades  ago  within  the  context  of  plasma 
physics  [85], [86], [87].  In  that  early  work,  the  average  dynamics  of  all  the  random 
quasi-particles  involved  were  treated  using  a  Vlasov  transport  equation. 

Subsequently,  incoherent  dark-stripe  and  dark-hole  (vortex)  solitons  were  exper¬ 
imentally  demonstrated  in  a  biased  photorefractive  crystal  [110].  In  all  cases,  these 
self- trapped  incoherent  beams  carried  the  characteristic  signature  of  dark  incoherent 
soliton  structures  [85], [86], [87], [109],  i.e.,  they  were  found  to  be  gray.  Moreover,  in 
agreement  with  predictions  of  Chapter  7  [109],  the  incoherent  dark  solitons  were 
experimentally  observed  when  an  appropriate  phase  pro  le  was  imposed  on  the 
wavefront.  Yet,  at  this  point,  several  important  questions  remain  unanswered.  First 
of  all,  are  there  truly  stationary  incoherent  dark  solitons  and  why  are  they  gray? 
Furthermore,  why  is  the  7r-phase  jump  necessary  for  their  excitation  and  how  is 
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it  possible  for  this  initial  phase  imprint  to  survive  in  the  midst  of  random  phase- 
uctuations?  The  answers  to  the  above  questions  can  not  be  obtained  from  the 
coherent  density  method  (because  of  its  inherent  complexity)  nor  from  the  approx¬ 
imate  Vlasov  approach.  These  issues  can  only  be  resolved  by  identifying  the  modal 
composition  of  these  dark  incoherent  soliton  states,  as  was  done  in  the  case  of  their 
bright  soliton  counterparts  [83], [114]- [116]. 

In  this  chapter,  by  means  of  an  exact  solution,  we  demonstrate  that  stationary 
incoherent  dark  sohtons  can  exist  in  nonUnear  self-defocusing  media.  These  sohtons, 
involve  in  general,  a  belt  of  radiation  modes  (both  odd  and  even)  as  well  as  bound 
states.  The  presence  of  even  radiation  and  bound  modes  explains  why  these  struc¬ 
tures  are  in  fact  gray.  Moreover,  we  nd  that  the  odd  radiation  modes  dominate 
within  the  dark  region  of  the  beam,  which  justi  es  the  7r-phase  shift  required  to  ex¬ 
cite  these  dark  incoherent  soliton  states.  The  coherence  properties  of  these  sohtons 
are  also  considered  and  they  are  found  to  be  in  good  agreement  with  the  results  of 
Chapter  7. 

9.2  Theoretical  formulation 

Let  us  consider  a  self-defocusing  nonlinear  medium  of  the  Kerr  type,  i.e.  = 
n\  —  n2/,  where  no  is  the  linear  refractive  index  of  the  material,  is  the  Kerr 
coefficient  and  /  the  optical  intensity.  We  also  make  the  important  assumption 
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that  the  nonlinearity  responds  much  slower  than  the  characteristic  phase  uctua- 
tion  time  across  the  beam  so  as  to  avoid  beam  breakup  due  to  speckle  instabilities 
[74], [83].  Thus  in  this  regime,  the  material  will  experience  only  the  time-averaged 
beam  intensity.  Such  non-instantaneous  Kerr-hke  media  include  for  example,  biased 
photorefractives  at  low  intensity  ratios  and  materials  with  appreciable  thermal  non- 
linearities  [23], [24],  [92], [93].  For  example,  a  typical  phase  uctuation  time  is  1  fis, 
whereas  a  photorefractive  crystal  responds  within  0.1  s  [110].  Let  the  time-averaged 
intensity  pro  le  of  this  planar  dark  incoherent  soliton  be  of  the  form 

Is  =  /o  [l  -  e^sech^  (9.1) 

where  the  parameter  <  1  is  associated  with  its  grayness  and  Xq  is  related  to  its 
spatial  extent.  The  partially-spatially  incoherent  dark  beam  is  quasi-monochromatic 
and  it  propagates  along  z.  Furthermore,  let  the  electric  eld  of  all  the  modes  com¬ 
prising  this  beam  be  written  as  E  =  U{x)  exp(i/3z)  where  0  is  the  mode  propagation 
constant.  Using  Eq.  9.1,  the  modal  function  U  is  then  found  to  obey  the  following 
Helmholtz  equation: 

^  +  (5  +  /sech2(a))c/  =  0  (9.2) 

where  s  =  x/xq,  g  =  [A:o(no  —  n^Io)  —  and  /  =  klxQe'^n2lo-  In  the  spirit 
of  Ref.  [83],  the  next  task  will  be  to  identify  an  appropriate  modal  composition 
such  that  the  time  averaged  intensity  fy  gives  rise  to  a  nonlinear  index  change 
which  is  self-consistent  [98]  with  the  composition  assumed  in  the  very  beginning.  In 
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general,  Eq.  9.2  exhibits  two  types  of  eigenfunctions:  radiation  modes  and  bound 
modes.  As  shown  schematically  in  Fig.9.1,  bound  states  are  possible  whenever 
g  =  -(^  ox  0^  >  kl{nl  -n-ilo),  whereas  radiation  modes  require  that  5  =  or 
<  kl{nl  -  na/o). 

At  this  point,  let  us  rst  assume  that  the  waveguide  induced  by  this  dark  beam 
can  only  support  one  bound  mode.  This  latter  requirement  can  be  met  provided 
that  the  coefficient  of  the  sech^(s)  potential  is  set  equal  to  two,  i.e.  /  =  2  or 
xg  =  2/ {kle^n2h)  [117], [118].  In  this  case,  all  possible  modes  allowed  by  Ekj.  9.2 
are  given  by  (see  Appendix  F)  [119]: 

Ub  =  sech(s)  (9.3) 

Ur,e  =  Qcos{Qs)  —  tanh(s)  sin((5s)  (9.4) 

Ur,o  =  Q  sin((5s)  +  tanh(s)  cos{Qs)  (9.5) 

For  f  =  2,  Ub  is  the  only  allowed  bound  state  (at  =  1)  and  the  two  degenerate 
eigenfrmctions  Ur,e  and  Ur,o  are  part  of  the  radiation  mode  continuum.  It  is  impor¬ 
tant  to  note  that  Ur,e  is  an  even  radiation  mode  whereas  Ur,o  is  odd.  Following  these 
results,  the  total  electric  eld  is  given  by  [117], [118] 

E  =  C6[4(s)  exp(i^tz)  (9.6) 

00 

+  j  dQ  [ce{Q)Ur,e{s,  Q)  -|-  Co{Q)Ur,o{s,  Q)]  exp  [ip0Q)z] 

0 
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Figure  9.1:  Eigenvalue  diagram  associated  with  a  first-order  incoherent  dark  soliton. 
The  bound  state  intensity  as  well  as  the  intensities  of  the  even  and  odd  radiation 
modes  (at  Q  =  0.1)  are  also  depicted.  The  dark  stripe  on  the  right  shows  the  spatial 
extent  of  the  soliton  induced  waveguide  when  e?  =  0.5. 
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where  Cb  and  Ce,o  are  modal  eld  coefficients  that  in  general  vary  randomly  in  time 
[83],[114]-[116].  In  Eq.  9.6,  the  upper  limit  of  the  integral  is  taken  at  in  nity.  As  we 
will  see,  Q  is  typically  in  the  neighborhood  of  Q  «  0,  far  away  from  Q^ax  =  - 

which  in  turn  justi  es  the  upper  limit  of  this  integral.  Under  incoherent 
excitation,  the  following  relationships  hold  true:  (cb,m<t,n)  ^  ^mn,  —  0) 

(coc^)  =  0  and  (ce(Q)c^(Q'))  =  (coiQKiQ'))  D{Q)6{Q-Q').  In  other  words,  the 
statistical  time  expectation  value  of  the  c-  eld  coefficients  is  zero  between  different 
bound  modes  and  the  same  applies  between  bound  modes  and  radiation  modes  (odd 
or  even).  Furthermore,  the  odd  and  even  radiation  modes  are  always  imcorrelated. 
Among  the  even  radiation  modes,  the  Cg  coefficients  correlate  only  for  the  same 
value  of  Q  and  this  is  also  true  for  the  odd  radiation  elds.  The  last  relationship 
also  implies  that  the  odd  and  even  radiation  modes  are  equally  excited  at  the  same 
Q  (with  strength  D(Q)).  This  is  because  the  random  source  shows  no  preference 
to  either  odd  or  even  radiation  modes.  The  positive  function  D{Q)  represents  a 
radiation  mode  distribution. 

By  utilizing  these  latter  relationships,  the  intensity  I  oc  {E{s,z)E*{s,z))  can 
then  be  obtained  from  Eq.  9.6,  i.e., 

OO 

I  =  Ahech\s)  +  j  D{Q)  +  tanh2(s)]  dQ  (9.7) 

0 

where  in  Eq.  9.7  we  made  use  of  the  fact  that  |f/r,e|^  +  lUr,<,|^  =  +  tanh^(s) 

and  a  A^.  The  rst  term  in  Eq.  9.7  arises  from  the  bound  mode  whereas 
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the  second  one  from  the  combined  intensity  of  odd  and  even  radiation  modes.  For 
self-consistency,  it  is  required  that  the  intensity  given  by  Eq.  9.7  is  identical  to  Is 
of  Eq.  9.1.  This  is  satis  ed  provided  that, 


OO 


/o  =  1  D{Q){Q^  +  l)dQ 

0 

(9.8) 

OO 

j  D{Q)[l-e'^{Q^  +  l)\dQ  . 

(9.9) 

0 


The  analytical  solution  given  by  Eqs.  9.8  and  9.9  clearly  demonstrates  that  station¬ 
ary  incoherent  dark  solitons  indeed  exist.  This  is  the  rst  time  we  know  of,  that  a 
soliton  was  found  to  involve  a  continuum  of  radiation  modes  as  well  as  bound  states. 
Even  more  importantly,  this  new  class  of  solitons  is  gray  because  of  the  presence 
of  even  bound  and  radiation  modes.  It  is  evident  from  Eqs.  9.8  and  9.9  that  the 
radiation  mode  distribution  function  D{Q)  is  by  no  means  unique.  In  fact,  in  nitely 
many  self-consistent  solutions  can  be  obtained,  depending  on  the  particular  choice 
of  D{Q). 

To  further  illustrate  our  results,  let  D{Q)  be  Boltzmann-like,  i.e.  D{Q)  = 
Do  exp{—Q/Qo)  where  Qo  represents  the  Q-width  of  this  distribution.  The  expo¬ 
nentially  decreasing  character  of  D{Q)  can  be  justi  ed  whenever  the  angular  power 
spectrum  of  the  incoherent  source  decreases  with  the  launch  angle  [79], [82], [109].  As 
a  result,  more  power  is  expected  to  be  coupled  into  small-angle  {Q  ^  0)  radiation 
modes  than  in  those  at  higher  Q  s.  For  this  speci  c  choice  of  D{Q)^  one  quickly 
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nds  that  Dq  =  {IqIQq){2QI  +  1)  ^  and  that  =  /q  [1  -  e‘^{2Ql  +  1)]  {2Ql  +  1)  ^ 
Thus,  Is  =  h  +  Ir  where 

h  =  5^1^  [l  -  sH2Ql  +  1)]  sech*(s)  (9.10) 

U  =  ,  l^Ql  +  1  -  sech^(s)]  (9.11) 

In  Eqs.  9.10  and  9.11,  /6  is  the  bound-mode  intensity  component  of  the  dark  inco¬ 
herent  soliton  and  Ir  is  the  intensity  pro  le  of  the  radiation-mode  belt.  It  is  also 
clear  from  Eq.  9.10  that  this  soliton  exists  provided  that  e'^{2Ql  4- 1)  <  1-  The  com¬ 
plex  coherence  factor  fiy  2(^1)  ^2)  [69]  of  this  incoherent  soliton  can  then  be  obtained 
from  Eq.  9.6,  by  evaluating  the  quantity 

{E{si,z)E*{s2,z))  oc  A^sech(si)sech(s2)  (9-12) 

00 

-H  I  DiQ)[Ur,e{Sl)UrAs2)  +  Ur,o{Sl)Ur,o{^2)]  dQ 
0 

00  1  l2 

In  turn,  its  correlation  length  can  be  found  from  Ids)  =  xq  f  /ii,2(S)  5  4- 5)  dS 

~0O  '  ‘ 

[109]. 


9.3  Results  and  discussion 

Let  us  now  physically  interpret  these  results.  Prom  the  /  =  2  condition,  one  can 
deduce  that  for  a  given  712/01  fhe  width  xq  of  the  dark  soliton  increases  with  its 
grayness.  Moreover,  it  is  important  to  note  that  Qo  de  nes  the  correlation  length 
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at  the  tails  (s  — +  ±cxd)  of  this  dark  incoherent  soliton.  In  these  regions,  the  bound 
states  disappear  and  the  soliton  correlation  length  is  determined  by  the  width  Qo 
of  the  radiation  mode  belt.  At  the  tails,  Ic  decreases  as  Qo  increases  and  vice  versa. 
In  fact,  Ic  in  these  regions  coincides  with  the  correlation  length  of  the  source.  In 
the  limit  Qo  =  0,  =  1,  Eqs.  9.10  and  9.11  reduce  to  the  well  known  coherent 

dark  spatial  soliton  solution  [113].  In  this  case,  h  —  0  and  the  soliton  consists  of 
an  odd  tanh{s)  mode  at  cut-off  with  Ic  =  oo  everywhere.  On  the  other  hand,  for 
Qo  =  0  and  7^  1,  we  obtain  an  incoherently-coupled  dark-bright  soliton  pair, 
identical  in  nature  to  that  previously  considered  in  photorefractive  crystals  [50], [53]. 
Prom  the  condition  e^{2Qo  -I- 1)  <  1,  it  is  also  clear  that  the  dark  soliton  becomes 
more  gray  as  its  incoherence  increases.  Another  interesting  possibility  arises  in  the 
limit  £^(2Qq  -I- 1)  =  1.  In  this  case  the  bound  state  is  empty  {h  =  0)  and  thus  the 
dark  incoherent  soliton  consists  of  only  radiation  modes.  As  previously  noted,  the 
dark  incoherent  soHton  is  actually  gray  because  of  the  presence  of  even  radiation 
and  even  bound  modes.  To  further  illustrate  these  issues,  let  us  consider  a  practical 
example.  Let  uq  =  2,  Aq  =  0.5  fim  and  712/0  =  10“^.  Let  the  soliton  grayness  be  50% 
or  =  0.5.  These  parameters  are  in  fact  close  to  the  those  previously  considered 
in  photorefractives  in  Chapters  7  and  8  [109], [110].  In  this  case,  Xq  c::;  5  fjLm  and 
this  sohton  exists  for  Qo  <  l/v^-  Figures  9.2(a)  and  (b)  show  the  soliton  intensity 
pro  le  and  correlation  length  when  Qo  =  0.4.  The  correlation  length  of  the  source 
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w  13.5  fim.  The  depression  in  4  at  s  w  0,  is  due  to  the  presence  of  the  bound  mode. 
Figs.9.3(a)  and  (b)  provide  the  same  information  when  Qo  =  0.7.  This  corresponds 
to  a  source  correlation  length  of  ~  5.3  fim  and  in  this  case  the  bound  mode  is  almost 
absent.  For  this  reason,  4  increases  around  the  dark  notch.  Overall,  the  behavior  of 
the  I3  and  Ic  curves  is  in  qualitative  agreement  with  the  ndings  of  previous  chapters 
[109].  From  the  above  results,  it  becomes  apparent,  that  for  our  choice  of  D{Q),  the 
radiation  modes  are  mostly  con  ned  within  a  narrow  belt  around  Q  ~  0.  Because  of 
this,  the  odd  radiation  modes  dominate  in  the  soliton-induced  waveguide  as  shown 
schematically  in  Fig.9.1.  This  behavior  can  be  easily  understood  by  considering  Eqs. 
9.4  and  9.5  in  the  neighborhood  of  s  w  0  when  Q  Ri  0.  Thus,  in  order  to  effectively 
launch  this  dark  incoherent  soliton  the  phase  must  be  properly  manipulated  so  as  at 
the  center  (s  =  0),  the  eld  distribution  is  mostly  odd.  This  explains  why  a  7r-phase 
shift  can  greatly  feicilitate  their  observation  [109], [110].  The  even  bound  mode  will 
subsequently  appear  as  a  result  of  evolution. 

Similarly,  higher-order  dark  incoherent  solitons  can  be  obtained  for  /  =  6, 12, 20, ... 
[117], [118],  [119].  For  example,  if  /  =  6,  the  soliton-induced  potential  can  support 

two  bound  states,  i.e.  Uh  =sech^(s)  and  UI2  =sech(s)  tanh(s)  at  gj  =  4  and  ql  =  1 

respectively.  In  this  case  (/  =  6),  the  radiation  modes  are  given  by 

C/^g  =  1^1  -h  —  3tanh^(s)j  cos(Qs)  —  3(5tanh(s)  sin(Qs)  (9.13) 

Ulg  =  \i  +  Q'^  —  3tanh^(s)j  sin(Qs)  4-  3(5  tanh(s)  cos(Qs)  .  (9-14) 
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Thus  I'h  =  J?^sech‘*(s)  +  C^sech2(s)  tanh^(s).  As  in  the  case  of  rst-order  incoherent 
dark  solitons,  D{Q)  is  not  unique.  If  we  choose  however,  D{Q)  =  DQexp{—Q/Qo), 
we  nd  that  =  (/o/2Qo)(12(5o  +  5Qo  +  2)“^ 

=  Qo£>o  [(3  +  6Q0)  -  2e"  (l2Qo  +  5Qo  +  2)] 

and  CP  =  8“^  +  9QoD'q.  The  intensity  component  of  the  radiation  belt  is  given  by 

j;  =  DoQo[2(l2(5^  +  5Q^  +  2)+9sech^(s)-(l2  +  6(5^)sech2(s)]  .  (9.15) 

Again,  1^  +  where  4  is  given  by  Eq.9.1.  Prom  the  above,  C  >  B  and  this 

soliton  is  possible  provided  that  >  0.  In  the  limit  S  =  0,  the  rst  (even)  bound 
mode  is  empty  and  the  soliton  involves  only  radiation  modes  and  the  next  odd  bound 
state.  Furthermore,  for  a  given  degree  of  grayness  and  nonlinear  index  change,  the 
/  —  6  dark  soliton  is  broader  than  the  rst-order  one  by  a  factor  of  \/3.  Figs. 9.4(a) 
and  (b)  depict  the  intensity  pro  le  and  correlation  length  of  such  a  second-order 
incoherent  dark  soliton,  under  the  same  parameters  used  before  when  Qo  =  0.5, 
^2  _  Q  49  2:0  8.8  /xm.  The  source  correlation  length  is  11  fjim.  As  shown  in 

Fig.9.4(b),  the  correlation  length  curve  now  exhibits  a  richer  sub-structure  within 
the  local  Ic  minimum  around  s  ~  0.  This  is  due  to  the  presence  of  the  additional 
odd  bound  mode.  Figs.9.2(b),  9.3(b),  and  9.4(b)  also  suggest  that  what  was  found 
in  the  computational  study  of  Chapter  7  was  actually  a  dark  incoherent  soliton 
of  the  rst-order  type.  More  speci  cally,  the  Ic  curve  of  the  dark  quasi-soliton  of 
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this  chapter  was  very  similar  to  that  of  Fig.9.3(b)  and  had  no  sub-structure  around 
s  «  0  which  is  characteristic  of  a  higher-order  dark  soliton  state.  It  is  also  very 
important  to  note,  that  in  this  latter  case  (/  =  6),  the  narrow  {Q  «  0)  radiation 
mode  belt  is  dominated  by  even  radiation  modes  within  the  waveguide  region.  On 
the  other  hand,  our  analysis  shows  that  this  second-order  dark  incoherent  soliton 
has  a  strong  contribution  from  the  second  odd  bound  state  which  by  the  way  is  never 
empty.  Such  higher  order  dark  solitons  can  be  launched  by  properly  engineering 
the  input  beam  in  a  manner  similar  to  that  of  Ref.  [120].  The  modal  composition 
of  even  higher-order  dark  incoherent  solitons  (/  =  12, 20, ...)  can  be  obtained  in  a 
similar  fashion.  These  results  are  however  more  involved. 


9.4  Conclusions 

In  conclusion  we  have  theoretically  demonstrated  the  existence  of  dark  incoherent 
spatial  solitons  in  non-instantaneous  self-defocusing  nonlinear  media.  These  new 
dark  soliton  entities  were  found  to  involve  radiation  modes  as  well  as  bound  states. 
Our  results  explain  for  the  rst  time  why  these  solitons  are  gray  and  why  a  7r-phase 
shift  tends  to  facilitate  their  experimental  observation.  Finally,  we  emphasize  that 
our  results  can  be  in  principle  extended  in  two  dimensions,  i.e.  in  the  description  of 
two-dimensional  dark  incoherent  solitons. 
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Chapter  10 


Conclusions 


In  this  thesis,  we  have  investigated  several  issues  regarding  incoherent  spatial  soh- 
tons  and  their  coherence  properties.  In  particular,  we  have  developed  two  com¬ 
plimentary  approaches  in  order  to  explain  the  observed  incoherent  self-focusing 
behavior  in  biased  photorefractive  media.  The  rst  one  is  the  so-called  coherent 
density  approach  where  the  underlying  evolution  model  takes  the  form  of  a  non¬ 
linear  Schrbdinger-like  integro-differential  equation  provided  that  at  the  origin  the 
coherent  density  is  appropriately  scaled  with  respect  to  the  angular  power  spectrum 
of  the  incoherent  source.  The  second  one  is  essentially  a  self-consistent  incoher¬ 
ent  multimode  approach.  In  brief,  in  this  procedure,  an  incoherent  spatial  soliton 
is  sought,  which  intensity- wise  is  a  superposition  of  all  the  modes  self-consistently 
guided  in  its  nonlinearly  induced  waveguide.  The  numerical  simulations  resulting 
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from  these  theoretical  approaches  were  found  to  be  in  good  agreement  with  the 
experimental  observations.  We  wish  to  emphasize  that  our  theoretical  approaches 
apply  not  only  to  photorefractive  media  but  also  to  any  other  nonlinear  system 
whose  temporal  response  time  is  much  bigger  than  the  phase  actuation  time  of  the 
incoherent  beam. 

We  have  shown  that  under  appropriate  initial  conditions,  bright  as  well  as  dark¬ 
like  incoherent  solitons  are  possible  in  biased  photorefractive  crystals.  Our  numerical 
simulations  have  demonstrated  that  the  coherence  properties  of  these  beams  are  sig- 
ni  cantly  affected  by  the  self-trapping  process.  More  speci  cally,  in  the  case  of  a 
bright  incoherent  beam,  we  have  found  that  coherence  length  remains  approximately 
constant  around  the  center  of  the  beam,  whereas  it  increases  at  the  margins.  The  in¬ 
crease  of  the  coherence  length  at  the  margins  of  the  beam  can  be  explained  from  the 
fact,  that  away  from  the  center,  only  the  higher  order  modes  contribute  to  the  coher¬ 
ence.  Our  numerical  simulations  also  show  that  further  increase  of  the  bias  voltage 
sets  the  bright  beam  into  cycles  of  compression  and  expansion.  In  this  latter  case, 
the  coherence  length  decreases  at  the  center  of  the  beam  during  compression  and  it 
increases  when  the  beam  expands.  Our  analysis  has  also  demonstrated  that  a  single 
dark  incoherent  soliton  can  be  excited  only  if  initially  a  7r-phase  shift  is  imposed 
on  the  incoherent  wave  front.  Moreover,  our  simulations  show  that  incoherent  dark 
solitons  are  in  fact  gray.  In  this  case,  the  coherence  length  was  foimd  to  be  higher 
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within  the  dark  notch,  with  a  depression  at  the  center.  Th  dynamics  of  incoherent 
dark  solitons  have  been  studied  under  different  initial  conditions.  Depending  on  the 
initial  conditions,  an  even  or  an  odd  number  of  incoherent  dark-like  structures  can 
be  obtained  in  self-defocusing  media.  It  is  important  to  note  that  the  behavior  of 
bright /dark  incoherent  solitons  is  fundamentally  different  from  that  of  their  coher¬ 
ent  counterparts.  For  example,  unlike  a  coherent  gray  soliton,  an  incoherent  dark 
sohton  does  not  exhibit  a  transverse  velocity  in  spite  of  its  grayness.  Furthermore, 
we  have  shown  that  the  evolution  of  incoherent  dark  solitons  in  non-instantaneous 
nonlinear  media  is  associated  with  strong  phase-memory  effects  which  are  other¬ 
wise  absent  in  the  linear  regime.  The  higher-order  behavior  of  these  dark  beams 
have  been  compared  under  the  same  initial  conditions  but  for  different  degrees  of 
coherence.  It  was  found  that  over  a  wide  range  of  parameters,  the  Y-splitting  is 
approximately  the  same  irrespective  of  the  initial  spatial  coherence.  Experimental 
observations  are  in  good  agreement  with  our  theoretical  predictions. 

In  addition  to  bright  solitons  in  photorefractive  media,  we  have  shown  that  par¬ 
tially  incoherent  spatial  solitons  are  also  possible  in  other  saturable  nonlinear  media. 
An  exact  Gaussian  solution  was  obtained  in  the  case  where  the  nonlinearity  is  of 
the  logarithmic  type.  We  found  that  these  incoherent  Gaussian  solitons  can  exist 
as  long  as  their  spatial  width  is  appropriately  interrelated  with  the  strength  of  the 
nonUnearity  and  the  width  of  the  incoherent  angular  power  spectrum  of  the  source. 
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The  behavior  of  these  beams  above  and  below  this  limit  was  also  investigated  us¬ 
ing  computer  simulations.  By  using  the  self-consistent  multimode  approach,  we 
have  shown  that  incoherent  spatial  solitons  in  logarithmically  saturable  nonUnear 
media  are  multimoded.  These  solitons  can  exist  as  long  as  their  mode-occupancy 
function  obeys  a  Poisson  distribution.  We  have  found  that  two  approaches,  that  is 
the  dynamic  coherent  density  description  as  well  as  the  self-consistent  multimode 
method  lead  in  this  case  to  exactly  the  same  results.  Even  though  the  saturable 
logarithmic  nonlinearity  differs  from  the  photorefractive,  it  provides  nevertheless  a 
platform  (perhaps  the  only  platform)  upon  which  the  equivalence  of  the  two  previ¬ 
ously  mentioned  approaches  can  be  established  in  closed  form.  In  the  same  vein,  we 
have  shown  that  two  dimensional  incoherent  Gaussian  solitons  can  also  exist  in  this 
nonlinear  system.  Moreover,  these  incoherent  solitons  can  be  elliptical  or  circular 
depending  on  whether  the  angular  power  spectrum  of  the  incoherent  source  is  sym¬ 
metric  or  not.  The  possibility  of  generating  two  dimensional  elliptical  solitons  in 
isotropic  nonlinear  media  seems  to  be  unique  to  incoherent  multimode  solitons,  since 
it  has  been  shown  that  their  elliptical  coherent  counterparts  change  their  widths  pe¬ 
riodically  during  propagation. 

We  have  also  shown  that  multimode  incoherent  spatial  solitons  are  possible  in 
non-instantaneous  Kerr-like  media.  Closed  form  solutions  were  obtained  using  the 
self-consistency  method  provided  that  their  intensity  pro  le  is  of  the  sech^  {xj xq) 
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type.  The  coherence  properties  of  these  incoherent  soliton  states  were  also  investi¬ 
gated  in  detail  and  explained  by  means  of  their  modal  composition.  Note  that  the 
cross-sections  of  the  complex  coherence  factor  ^12  were  very  different  in  character 
from  those  found  in  the  case  of  logarithmically  saturable  nonlinear  media.  In  Log 
systems,  the  statistical  process  is  everywhere  stationary,  i.e.  /ii2  depends  only  on 
<5  and  it  is  Gaussian  in  nature.  For  the  Kerr  case,  the  coherence  curves  depend 
also  on  the  position  r]  and  do  not  always  go  to  zero  as  6  — >  ±00.  This  feature  was 
also  encountered  in  saturable  Kerr  media  of  the  type  I /  (1  -I-  /).  This  important 
difference  is  due  to  the  nite  number  of  modes  associated  with  the  Kerr  incoherent 
solitons.  In  particular,  the  induced  waveguides  in  Kerr  and  Kerr-saturable  media 
exhibit  cutoffs  (  nite  number  of  modes)  whereas  the  logarithmic  one  does  not. 

Finally,  we  have  theoretically  demonstrated  the  existence  of  dark  incoherent 
spatial  solitons  in  non-instantaneous  self-defocusing  nonlinear  media  of  the  Kerr 
type.  These  dark  soliton  entities  were  found  to  involve  radiation  modes  as  well 
as  bound  states.  The  presence  of  even  radiation  and  bound  modes  explains  why 
these  structures  are  in  fact  gray.  Moreover,  we  have  foimd  that  the  odd  radiation 
modes  dominate  within  the  dark  region  of  the  beam  (for  rst-order  incoherent  dark 
beams),  which  justi  es  the  vr-phase  shift  required  to  excite  these  dark  incoherent 
soliton  states.  Thus,  in  order  to  effectively  launch  dark  incoherent  solitons  the 
phase  must  be  properly  manipulated  so  as  at  the  center,  the  eld  distribution  is 
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mostly  odd.  The  even  bound  mode  subsequently  appears  as  a  result  of  evolution 
The  coherence  properties  of  these  solitons  were  also  considered  and  they  were  found 
to  be  in  good  agreement  with  the  results  of  numerical  simulations  in  the  case  of 
photorefractive  nonlinearity.  The  existence  of  such  incoherent  dark  solitons  was 
later  veri  ed  experimentally. 

In  closing,  we  would  Uke  to  note  that  other  issues  may  also  merit  further  investi¬ 
gation.  For  example,  as  in  the  coherent  case,  the  dynamics  of  two  dimensional  Log 
type  incoherent  bright  solitons  can  be  described  in  closed  form.  Furthermore,  our 
results  for  one  dimensional  incoherent  dark  solitons  can  be  in  principle  extended 
in  two  dimensions  to  analyze  two-dimensional  dark  incoherent  solitons.  Detailed 
investigations  of  coherent  and  incoherent  soliton  collisions  might  be  another  inter¬ 
esting  avenue  for  future  studies.  The  outcome  of  collisions  from  incoherent  beams 
is  expected  to  be  more  involved  since  incoherent  beams  are  multimoded.  It  also 
remains  to  be  seen  how  one  can  control  the  coherence  of  an  optical  beam  by  using 
incoherent  beams. 
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Appendix  A 

Normalized  mode  functions 

In  this  Appendix  we  provide  the  normalized  mode  functions  (77).  Here  we  de  ne 
T  =  tanh  (77)  and  S  =  sech{ri).  The  (77)  functions  appear  in  accord  to  their,  modal 
order  (lowest  order  mode  rst). 

For  71  =  1 : 


u\  (77)  =  S  (A.l) 

For  n  =  2: 

u\  (77)  =  5^  (A.2) 


u\{r})  =  ST 


(A.3) 
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For  n  =  3: 


ul  (77)  =  ^Jl5/16S^ 

ul  iv)  = 

(77)  =  (1/4)  5  (4 -552) 

For  77  =  4: 

as  (»)  =  /vss- 

nUri)  =  {i/i)V7S^T 
ul  (v)  =  (v^/4)  S’  (6  -  7S’) 
5SW  =  {l/4)Sr(4-7S’) 

For  n  =  5: 

2^  (rj)  =  (x/^/16)  5® 


(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.IO) 

(All) 
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ul  iv)  =  (1/2) 

(A.12) 

ul  (tj)  =  (^/16)  (8  -  952) 

(A.13) 

(^)  =  (V7/2)  (2  -  352) 

(A.14) 

u\  {ri)  =  (1/8)  5  (8  -  285^  +  215^*) 

(A.15) 

56  {ri)  -  (3/16)  \/^5® 

(A.  16) 

u|(77)  =  (5/16)  \/^5®T 

(A.17) 

«6  (^)  =  (^/3/4)  54  (10  -  1152) 

(A.18) 

ul  (77)  =  (3/16)  \/io53r  (8  - 1152) 

(A.19) 

(t,)  =  (^/Io/16)  52  (16  -  485^  +  335“) 

(A.20) 

[t])  =  (1/8)  ST  (8  -  365^  +  335'*) 


(A.21) 
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Appendix  B 


Complex  coherence  factors 


In  this  Appendix,  the  functional  form  of  the  complex  coherence  factors  is  pro¬ 
vided  up  to  n  =  5.  Again  we  have  de  ned  T  (77)  =  tanh  (77)  and  S  (77)  =  sech(77). 

For  n=\: 


l‘n(v,V  +  ^)  =  1 


(B.l) 


For  71  =  2: 


fJ-n  {V,V  +  S)=T  (77)  T{rf  +  S)  +  S{r])S{r]  +  S)  (B.2) 


For  71  =  3; 


f^i2{v,V  +  S)  =  1  + (5/2)  {5  (77)  5  (77 +  5)  (B.3) 
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xlT(fl)T(7)  +  ii)+SMS(i7  +  «)|  (B.4) 

-(l/2)[5"W  +  Sn,  +  fl]} 


For  n  =  4: 


Ati2  (r?, 7?  +  5)  =  T(r))T{r]  +  6)  i 


+S{v)S{v  +  S){ 


1  -  (7/4)  [52  (7?)  +  52  (7?  +  6)] 

+752(77)52(77  +  6) 

9/2 -(21/4)  [52  (77) +  52  (77 +  6)] 
+752  (77)52  (77  +  6) 


(B.5) 


For  77  =  5: 


f^u{v,V  +  S)  =  1  +  7{T(77)T(77  +  6)  5(77)5(77  +  6)  (B.6) 

X  [1  -  (3/2)  (52  (77)  +  52  (77  +  6))  +  352  (77)  52  (77  +  5)] 

-  (1/2)  (52  in)  +  52  (77  +  6))  +  (3/8)  (5"  (77)  +  5^*  (77  +  6)) 
13/4-3(52  (77) +  52  (77 +  6)) 

+352  (77)52(77  +  6) 


+52  (77)  52  (77  +  6) 


} 
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Appendix  C 


Diffraction  equations  of 
Gauss-Hermite  beams 


This  Appendix  details  the  steps  that  lead  to  the  diffraction  equation  of  Gauss- 
Hermite  beams. 

An  optical  beam  propagates  in  a  linear  material  by  obeying  the  following  differ¬ 


ential  equation: 


.dU  1  d'^U  ^ 


(C.1) 


where  the  normalized  coordinates  are  ^  =  z/kxo,  s  =  x/xq,  and  k  —  koUo  is  the 
wavenumber  in  this  material.  Let  as  assume  that  at  ^  =  0  the  optical  eld  is  given 
by  Um  =  Aoexp{-s^/2}  Hm{s)  where  are  Hermite  functions.  Now,  let  the 
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diffraction  solution  be 


Ur, 


=  A  (C)  exp  {-s^l2u?  (O}  Hm  {s/u  (0)  exp  \i  (0  +  F  (0  s^] }  (C.2) 


where  A,  w,  9,  F  are  real  functions.  Substituting  this  solution  into  Equation  C.l 
gives 

(i) 


+ 


''->2  2  AFs 


-  AF^s^  -  i- 


^  +  i2F\H„ 


2 

H — 

OJ 


s 

- -  > 

(jT 


where  A  =  dA/d^,  etc.,  H'^{v)  =  dHm/dri,  =  (FHm/drf  where  ri  =  s/w. 

Prom  the  imaginary  part  of  this  equation  we  get 


2- -4 

L  u  u)^ 


AH's 


2—  -  ^ 
L  u>^ 


[A  +  FA]Hm  + 

By  equating  the  coefficients  of  the  s-polynomial  to  zero  we  get 

A  + FA  =  0 


AHms"^  =  0  (C.4) 


(C.5) 


and 


OJ 


U) 


=  2F 


From  the  real  part  of  Equation  C.3  we  obtain 


-AHm{9  +  Fs‘^)+^ 


2uj‘^ 


-  2AF^s^  - 


A 

2a;2 


A 

r  «  25  1 

rjff  Tjf 

(C.6) 


=  0  (C.7) 
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This  relation  can  be  further  simpli  ed  by  using  the  following  de  nition 


^  dHm  „  „ 

2rj— —  =  —2mHr, 


drf 


dr] 


and  we  end  up  with 


-[e  +  Fs^)  AHm  + 


1  2_2FV-^ 

2uj^  2w^ 


AHm  -  \mAHm  =  0 


(C.8) 


(C.9) 


Thus 


^  =  “2:? 


F  +  2F2  -  —  =  0 

2w4 


^(jj  =  constant  = 


(C.IO) 

(C.ll) 

(C.12) 


where  we  assume  that  the  initial  spot  size  luq  is  equal  to  1,  i.e.  w  =  0)  =  1.  Eq.C.6 
can  be  written  in  the  following  form 


F  =  U)  . 

4a;2  ^  V  ; 


(C.13) 


Substituting  this  result  into  Equation  C.ll,  we  obtain  the  following  relationship 


UJ  = 


(C.14) 


Prom  this  equation  one  can  solve  for  u.  This  is  given  by 


1/2 


(C.15) 
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Furthermore,  using  this  solution  in  Eq.C.lO,  one  may  calculate  the  phase  accumu¬ 
lated  ^  (0- 


=  (C.16) 


where  9o  is  the  initial  value  of  the  phase  associated  with  this  beam.  Again  substi¬ 
tuting  Eq.C.15  into  Eq.C.12  one  directly  nds  the  amplitude  of  the  diffracting  mode 
which  is  given  by 


-^0 


{i+ey'" 


(C.17) 


One  can  also  nd  the  function  F  (C  by  substituting  Eq.  C.15  into  Eq.  F.ll  to  get 


F(i)  = 


2(l  +  ?") 


(C.18) 


Substituting  A,  u,  9,  F  into  Eq.  C.2  and  taking  the  absolute  value  square,  the 
intensity  of  a  diffracting  Gauss-Hermite  beam  is  given  by: 


In  general,  if  the  initial  value  of  the  spot  size  is  not  unity,  i.e. 
Elqs.C.lS  takes  the  following  form 


\ 

(C.19) 

/ 

u;  Wo  at  ^  =  0, 


w  (^)  =  Wo  (l  +  C^/a^o)  (C.20) 


and  the  diffracting  beam  is  given  by 
Ao 


\Un.{0\  = 


(l  +  1  (l  +  f /‘^o) 


Hi 


180 


Appendix  D 


Intensity  and  coherence  function 
of  incoherent  Gaussian  becuns 


In  this  Appendix,  by  using  modal  description,  we  compute  the  time-averaged  in¬ 
tensity  and  complex  coherence  function  of  an  incoherent  Gaussian  beam  during 
diffraction. 

The  optical  eld  of  an  incoherent  Gaussian  beam  can  be  expressed  through 

OO 

superposition,  i.e.,  U  —  Yl  Here  the  Gauss-Hermite  modes  at  ^  =  0  are 

m=0 

given  by 

Um  =  Hm  exp  (-as^/2)  (D.l) 
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at  ^  =  0  whereas  they  are 


— 


1 

{l  +  a^ey^^ 

xHm 

f 

a 


[l  + 

172  I  exp  |i 


(D.2) 


0(0+  r 


2  [1  +  ' 

at  ^  7^  0  from  Appendix  C.  Furthermore,  Cm  are  the  mode-occupancy  coefficients 
that  vary  randomly  in  time  and  let  the  mode-occupancy  be  described  by  a  Poisson 
distribution,  that  is  (|cm|^)  =  rexp  (-p/2)  (p/2)"'  /m!  where  p  is  the  Poisson  para¬ 
meter  and  f  =  r  (1  —  p^)  exp  —p/2  is  a  constant.  The  time-average  intensity  of  this 
beam  can  then  be  obtained  from 


00 

Id  =  (PI')  =  E  (D-3) 

m=0 

where  we  made  use  of  the  fact  that  the  time  average  of  the  cross- interference  terms  is 
zero  {(^CiCj^^ij)  under  incoherent  excitation.  Therefore,  the  time-average  intensity 
of  this  multimoded  beam  during  linear  diffraction  is  given  by 


rexp 


X  exp  <  — 


This  expression  can  be  simpli  ed,  by  using  Mehler  s  formula: 


^  (i/2)’"  ..2  ,  X  1  f2ta:2-2tV] 


(D.4) 

(D.5) 


(D.6) 
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where  t  <  1.  Thus  as  a  result  one  nds 


{l  +  a'^ey 


l  +  L  1  +  P 


Knowing  that  a  (1  -  p)  /  (1  +  p)  =  1,  we  further  simplify  Eq.D.7  which  gives 


(l  +  a^e) 

The  coherence  function  is  de  ned  as  follows: 


'l  +  a2^2]J  ■ 


Mi2  (si,S2,0  = 


{Uisu0U*{s2,0) 

[Id{Su0Id{S2,0]'^" 


Therefore,  in  our  case 


{U{Si,0U*{S2,0)  =  ( 


oo 

E  (|CmO“m(Sl,0«m  («2,0 


(D.IO) 


since  the  time  average  of  the  cross-interference  terms  is  zero.  Substituting  the 
functional  form  of  Um  (si,  and  {s2, 0  above  equations,  and  by  using  again 
Mehler  s  Formula,  and  after  some  simpli  cation,  we  obtain  the  following  complex 


coherence  function 


=  (D.ll) 


where  s  =  si  and  5  =  S2  ~  51. 
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Appendix  E 


Diffraction  of  an  incoherent 
Gaussian  beam 


In  this  Appendix,  by  using  coherent  density  approach,  we  analyze  the  diffraction  of 
an  incoherent  Gaussian  beam. 

The  diffraction  equation  in  normalized  coordinates  and  quantities  is  given  by: 


i 


2  5s2 


(E.1) 


where  s  =  x/xq,  C  =  z/kxo,  k  =  fcono,  and  «  =  kxoO.  Now  let  ^  =  C  and  r]  =  s-a^. 
After  this  transformation,  the  above  equation  takes  the  following  form 


Let  at  the  input,  the  coherent  density  function  be  given  by 


/  (C  =  0.  ^)  =  exp  ^-77^/2)  ^  exp  . 


Therefore  at  C  7^  0,  the  intensity  of  this  function  evolves  according  to  following 
expression 


or 


\f{Cv)f  = 


(E.3) 


(E.4) 


The  total  intensity  of  the  incoherent  beam  can  be  calculated  from  Id  =  \Uf  — 
7  \f  (s,  d9.  Substituting  Eq.E.4  into  this  latter  equation  gives 


Id  = 


r  exp 


where  V  =  kxoOo  and  A  =  6/60.  The  integral  of  the  above  expression  can  be 
simpli  ed  by  using  /  exp(— +  bx)  =  vr^^^/oexp  (6^/4a^).  Therefore,  the  total 

—  00 

intensity  of  a  incoherent  Gaussian  beam  during  diffraction  is  given  by 
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Appendix  F 
Radiation  modes 


In  this  Appendix,  we  compute  the  radiation  modes  of  the  Helmholtz  equation  when 
the  potential  is  of  the  sech"^  (s)  type. 

In  this  case,  the  Helmholtz  equation  is  given  by 

^+[5  +  /sech2(s)]l7  =  0  (F.l) 

where  g  is  required  to  be  a  positive  number,  i.e.  g  =  Q^,  for  radiation  modes.  We 
assmne  that  the  induced  waveguide  can  support  only  one  bound  mode,  i.e.  /  =  2. 
In  that  case  this  equation  takes  the  following  form 

^  +  [q2  +  2  sec  h^  (s)]  U  =  0.  (F.2) 

To  nd  the  odd  radiation  modes,  let  us  assume  that 

U  =  Uo  =  F{s)  cos  (Qs)  +  G  (s)  sin  (Qs)  (F.3) 
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After  substituting  this  solution  into  Equation  F.2  and  by  selecting  the  cos((5s) 
terms  we  nd 

F  +  2(5G  +  2sech^  (s)  F  =  0  (F.4) 

whereas  by  selecting  sin  {Qs)  terms  we  get 

G  -  2QF  +  2  sec  h^  (s)  G  =  0  (F.5) 


Let  G  =  Go  is  a  constant.  Therefore  Eqs.F.5  and  F.4  can  be  simph  ed  to  give 


F  +  2  sec  h^  (s)  F  =  0 
G 

®  sec  h^  (s) 
or 

sec  h  [s)  =  -p:;- 
Go 

Substituting  Eq.F.8  into  Eq.F.6  results  with  the  following  relation 

F+^F  =  0 
Go 


(F.6) 

(F.7) 


(F.8) 


(F.9) 


Integrating  the  last  equation  once  yields 

F+^F2  =  G  (F.IO) 

Go 

Moreover,  from  its  asymptotic  behavior  we  know  that  F  (s  ^  ±oo)  =  ±1  and 
F{s  ->■  ±oo)  =  0.  The  constant  G  can  be  found  by  using  these  boundary  con¬ 
ditions  and  letting  Go  =  Q  and  it  is  G  =  1.  Now,  Eq.F.lO  can  be  integrated  one 
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more  time  to  nd  F  (s)  which  is  given  by 


F  (s)  =  tanh  (s)  . 


(F.ll) 


Finally,  substituting  the  functional  form  of  F  (s)  and  constant  value  oi  G  =  Q  into 
Eq.F.3  yields 

Uo  =  Q  sin  (Qs)  +  tanh  (s)  cos  (Qs)  .  (F- 12) 

The  even  radiation  modes  can  be  generated  from  the  odd  one  by  letting  t4  =  uUo 
and  by  using  variation  of  parameters.  By  doing  so,  we  nd 

=  =  (F.13) 


and 


U,  =  [Q  sm  {Qs)  +  tanh  {»)  cos  (Q»)|  /  |q  W  ,<,3  (F- 1^) 


Therefore  the  even  radiation  modes  are  given  by 


Ue  =  Q  cos  (Qs)  —  tanh  (s)  sin  (Qs)  (F-15) 

Since  UoUe  —  UoUg  =  —Q^  —  Q  ^  0  these  two  degenerate  modes  are  Uneaxly  inde¬ 
pendent.  Moreover  U^  +  Ug  —Q'^  +  1,  i.e.  constant,  at  s  ±oo. 
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